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1.  INTRODUCTION 


A  laboratory  Tilt  Rotor  Flying  Machine  (TROFM)  was  designed  and  built  for 
the  purpose  of  conducting  an  experimenttj  study  of  the  effects  of  comm2uid  systems, 
displays  auid  tasks  on  human  pilot  behavior.  The  aircraft  is  mounted  on  a  two  axis 
pivot  stand  which  gives  the  fuselage  two  degrees  of  freedom;  pitch  and  yaw.  Figure  1.1 
shows  a  diagram  of  the  TROFM.  Rotation  of  the  mast  is  referred  to  as  yawing  motion 
and  is  denoted  Angular  displacement  of  the  pipe  about  the  pivot  point  is  called 
pitch  which  is  denoted  0.  The  main  rotor  is  tilted  by  a  Futaba  9102  servo.  Tilt  or 
roll  angle  is  denoted  $.  Main  rotor  speed  is  the  primary  pitch  axis  controller  while 
the  rotor  tilt  angle  is  used  for  yaw  axis  control. 

The  angles  vp,  0  and  $  are  measured  by  potentiometers  (1  Kfli  %1  linearity) 
configured  as  voltage  dividers  powered  by  a  10  V,  200  mA  power  supply  with  %0.5 
ripple.  A  potentiometric  joystick  is  also  configured  as  a  voltage  divider  and  shares  the 
sensor  power  supply.  The  rotor  is  driven  by  a  Maxon  18  VDC  motor  which  is  powered 
by  a  Copley  model  201  Pulse  Width  Modulated  (PWM)  amplifier.  The  motor  has  a 
DC  tachometer  mounted  on  its  base  for  rotor  velocity  measurements.  The  rotor  is  a 
12  inch  fixed  (#6)  pitch  model  airplane  propeller. 

Sensors  and  actuators  are  interfaced  to  a  Mac  Ilsi  digital  computer  via  a  National 
Instruments  NB-MIO-16  A/D,  D/A  board.  The  board  has  12  bit  accuracy,  9/i  sec 
A/D  conversion  time  and  a  maximum  sampling  rate  of  100,000  samples/sec.  Two 
analog  output  channels  are  used  to  send  commands  to  the  actuators  and  six  differ¬ 
ential  analog  inputs  are  used  to  sample  sensor  and  joystick  signals.  The  pilot  in  the 
loop  experiments  make  use  of  a  custom  built  X-Y  plotter  which  uses  four  digital  I/O 
lines  on  the  NB-MIO-16  board  to  drive  two  stepper  motors 
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Figure  1,1  Schematic  of  Tilt  Rotor  Flying  Machine  (TROFM). 


Analytical  models  are  derived  using  Lagrangian  and  Newtonian  techniques.  Iden¬ 
tification  of  the  aircraft  characteristics  is  based  on  a  non-p£U'ametric  Discrete  Fourier 
Transform  technique  and  a  maximum  likelihood  parameter  estimation  scheme.  Clas¬ 
sical  control  design  techniques  are  used  to  design  continuous  compensators  for  rate 
and  attitude  command  systems.  Digital  filters  are  designed  to  emulate  the  continuous 
compensators  and  are  implemented  in  real  time  on  the  digital  computer.  Rate  com¬ 
mand,  attitude  command  and  open  loop  (unaugmented)  systems  are  developed.  The 
attitude  command  system  use  feedback  from  the  pitch  and  yaw  sensors  to  maintain 
a  desired  attitude  specified  by  the  position  of  the  pilot’s  joystick.  The  rate  command 
system  uses  rate  estimates  for  feedback  to  maintain  angular  rates  specified  by  the 
position  of  the  joystick.  The  rate  estimates  are  generated  by  a  digital  Kalman  filter 
which  is  based  on  a  complete  linear  dynamic  model  of  the  TROFM  because  rate 
sensors  are  unavailable,  attitude  sensors  are  too  noisy  to  be  differentiated,  and  the 
model  does  not  account  for  non-linear  phenomena  and  unknown  disturbances. 

A  study  of  the  effects  of  command  systems,  displays  and  tasks  on  human  pilot 
behavior  is  conducted.  Loren  P.  Dean,  general  aviation  pilot  used  the  Cooper-Harper 
rating  scale  [2]  to  evaluate  the  flying  qualities  of  the  TROFM  using  three  command 
systems  for  two  types  of  tasks  and  two  different  display  types.  Pilot  and  pilot-vehicle 
describing  functions  are  obtained  from  time  histories  of  dynamic  trficking  tasks  where 
the  pilot  attempts  to  track  a  random  appearing  input.  Compensatory  instrument  dis¬ 
plays  and  a  visual  pursuit  scene  are  used  for  the  dynamic  tracking  tasks.  Correlations 
between  the  Cooper-Harper  ratings  and  quantitative  descriptions  of  pilot  behavior  are 
made.  The  experimental  results  are  compared  to  the  classical  theory  of  man-machine 
system  dynamics  which  is  based  on  the  crossover  model  [9]. 


3 


2.  EQUATIONS  OF  MOTION 


In  order  to  develop  control  laws  for  a  physical  system,  a  mathematical  model 
of  the  motion  of  the  system  must  be  developed.  For  simple  systems  a  Newtonian 
approach  is  advisable,  but  for  more  complex  systems  a  Lagrangian  approach  is  more 
suitable,  since  it  allows  one  to  directly  obtain  the  differential  equations  of  motion  of 
a  physical  system  from  its  kinetic  energy  expression. 

The  physical  system  of  interest  is  a  Tilt  Rotor  Flying  Machine  (TROFM)  whose 
angular  position  is  controlled  by  an  aerodynamic  actuator  or  rotor.  The  TROFM  is 
free  to  rotate  in  pitch  and  yaw  about  a  fixed  point.  The  pitch  angle  is  controlled  by 
the  lift  generated  by  the  rotor  and  the  yaw  angle  is  controlled  by  tilting  the  rotor  in 
the  direction  of  desired  rotation.  A  schematic  of  the  machine  is  shown  in  Figure  1.1. 

An  idealization  of  the  TROFM  consists  of  four  rigid  bodies; 

1  M.  Mast  assembly 

2  P.  Pipe,  ballast,  tilt  servo  ,  tilt  angle  transducer,  idler  wheels  and  mount¬ 
ing  bracket. 

3  T.  Tilt  mechanism,  i.e.  friction  wheel,  motor  mount,  motor  stator. 

4  R.  Rotor,  motor  armature. 

One  inertially  fixed  and  four  body  fixed  coordinate  axis  systems  are  defined.  In 
order  to  use  Lagrange’s  equations,  the  kinetic  energy  of  the  entire  system  must  be 
expressed  in  inertial  coordinates;  therefore,  relations  between  the  body  axis  systems 
and  an  inertial  frame  were  developed. 
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2.1  Coordinate  Systems 


An  inertial  coordinate  system  (Figure  2.1)  i  with  origin  at  the  top  of  the  mast  is 
defined.  The  ix,iy  and  i,  unit  vectors  are  earth  fixed  and  defined  positive  North,  Elast 
and  toward  the  center  of  the  earth  respectively. 

A  Cartesian  coordinate  system  is  defined  which  rotates  with  the  mast  and  whose 
origin  is  at  the  top  of  the  mast.  The  mast  is  free  to  rotate  an  angle  ^  about  the  iz 
axis.  The  maist  body  axis  system  is  denoted  m  and  is  shown  in  Figure  2.2. 

A  vector  quantity  v  may  be  written 

V  =  mx  rhx  +  m„  m„  +  m,  m,  [2.1] 

—  "b  *y  d" 

where  rhi,  rhy,  and  are  unit  vectors  coincident  with  the  body  fixed  axes  of  the 
mast.  The  scalars  m^,  rriy,  and  iy,  are  the  x,  y  and  z  components  of  the 
vector  V  measured  with  respect  to  the  m  and  i  coordinate  systems  respectively.  It  is 
convenient  to  make  use  of  linear  algebra  when  showing  relationships  between  vectors 
and  deriving  expressions  for  the  kinetic  energy  of  a  rigid  body;  therefore,  matrix 
notation  is  introduced.  The  relationship  between  m*,  rriy,  m*  and  i*,  »y,  is  given 
by 


mi 

cos(^) 

—  sin(^) 

0 

»x 

ruy 

= 

sin(^) 

cos(^) 

0 

*y 

[2.2] 

rriz 

m 

0 

0 

1 

ix 

i 

A  pipe  or  p  coordinate  axis  system  is  defined  which  rotates  with  the  second  rigid 
body  about  an  origin  fixed  at  the  top  of  the  mast.  This  axis  system  is  free  to  rotate 
an  angle  ^  about  iz  and  an  angle  0  about  rhy.  The  relationship  between  the  p  and 
i  coordinate  system  is  illustrated  in  Figure  2.3. 

A  vector  quantity  v  may  be  written 

V  =  Px  Pi  +  p»  Py  +  Px  Px  [2.3] 

—  lx  “1“  ly  ly  “I"  Iz  Iz 
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The  following  equation  shows  the  mathematical  relationship  between  the  x,  y,  and  z 
components  of  v  measured  with  respect  to  the  p  body-axis  system  and  the  x,  y  and 
z  components  of  v  measured  with  respect  to  the  inertial  frame. 


• 

• 

cos(0)cos(’i') 

cos(0)8in(^?)  —  8in(0) 

tx 

Py 

= 

-sin(^') 

cos(4')  0 

Pz 

•  9 

P 

sin(0)co8(4') 

sin(0)sin(4’)  cos(0) 

iz 

»  m 

(2.4] 


Next  a  coordinate  system  t  (Figure  2.4)  is  defined  which  rotates  with  the  tilt 
mechanism  and  has  an  origin  located  at  the  top  of  the  mast.  The  axis  system  is  free 
to  rotate  an  angle  4  about  t,,  an  angle  0  about  ihy  and  an  angle  $  about  the  pg 
axis. 

The  following  matrix  transforms  components  of  vectors  reckoned  with  respect  to 

the  inertial  frame  to  components  reckoned  with  respect  to  the  tilt  body-axis  system: 
<pi  to  t  _ 


eoi(6)on(<l')  oai(6)Mii(<l)  -iia(0) 

«n(4)sia(d}ea((#) -coa(4}«m(<l)  aaf#)tia(9)na(*l)  iia(4)cat(6) 

c(M(4)siii(9)oai(<f)’fsin(4)siii(<l)  o(»<#>sin(8)sia(<l) -sia(4)oai<<t)  om(4)cw(6) 


Finally  a  rotor  body-axis  system  r  is  defined  whose  origin  is  located  at  the  center 
of  mass  of  the  main  rotor  and  motor  armature  and  with  unit  vectors  aligned  with  the 
tilt  body-2kxis  system.  Figure  2.5  shows  the  relationship  between  the  rotor  body-axis 
system  and  the  inertial  frame. 

Consider  a  vector  quantity  v  such  as  position  or  angular  velocity.  Let  Vr  be  the 
vector  quantity  measured  with  respect  to  the  rotor  body-axis  system  and  Vt  be  the 
vector  quantity  measured  with  respect  to  the  tilt  body-axis  system. 

Vr  =  r.  f,  -I-  r,  f„  +  r,  f,  (2.6) 

Vt  =  fr  +  ty  iy  +  it  U 
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Figure  2.4  Orientation  of  the  tilt  body-axis  system  with  respect  to  the  inertial 

frame. 
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The  vector  v  as  observed  from  an  inertial  frame  is  given  by 

V  =  V,  +  Vt  [2.7] 

But  since  t  and  r  are  parallel, 

V  =  (rx  +  tx)  tx  +  (ty  +  ty)  ty  +  (r,  t* )  t,  [2-8] 

The  relationship  between  the  components  of  v  reckoned  with  respect  to  the  tilt 
axis  system  and  the  components  of  v  reckoned  with  respect  to  the  inertial  frame  are 
given  by 

“h  ^x  ®x 

r^^ty  [2.9] 

r,  +  i. 

Consider  the  position  of  the  origin  of  the  rotor  body-axis  system  relative  to  the  tilt 
body-axis  system.  Note  that  this  position  is  constant.  Let  lx,ly  and  /*  denote  the  x, 
y  and  z  distances  measured  from  the  origin  of  the  tilt  body  system  to  the  origin  of 
the  rotor  body  system  (Note  ly  =  0).  These  constants  are  useful  in  formulating  an 
inertia  tensor  for  the  rotor  center  of  mass. 

A  relationship  between  the  angular  velocities  of  a  rotating  rigid  body  axis  system 
and  a  set  of  generalized  inertial  angular  velocities  (i.e  Euler  angles)  is  also  needed. 
Figures  2.2,  2.3,  and  2.4  show  the  orientations  of  the  Euler  angular  velocities  with 
respect  to  the  inertial  frame.  The  Euler  angles  are  considered  inertial  because  they 
uniquely  define  the  orientation  of  a  body.  Body  fixed  angular  velocities  are  non- 
inertial  since  they  do  not  uniquely  define  the  orientation  of  the  body.  Consider  a 
body  which  has  three  rotational  degrees  of  freedom  and  no  translational  degrees  of 
freedom  (i.e.  tilt  body).  Assume  that  the  origin  of  the  body-axis  system  and  the 
origin  of  the  inertial  frame  are  coincident.  The  angular  velocity  of  the  body  reckoned 
with  respect  to  the  inertial  frame  ‘u;*’  may  be  written  as 

hx  -f  by  •¥  bi  [2-10] 
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[2.111 


=  **  +  *»  +  ** 

=  i,  +  0  +  $  Pr 

The  mathematical  relationship  between  the  body-axis  angular  velocities  and  the 
Euler  angular  velocities  is  given  below  [6]. 


wj,  1  0  —  sin(0) 

=  0  cos($)  cos(©)sin($)  |  |  0  [  [2-12] 

0  —  sin($)  cos(0)cos($) 

j  b  L 

Note  that  the  tilt,  pipe  and  mast  body  angular  velocities  can  be  transformed  to 
Euler  angultir  velocities  by  direct  application  of  Elquation  2.12. 


$ 

0 

4/ 

2.2  Kinetic  Energy 


In  general  the  total  kinetic  energy  of  a  system  of  n  rigid  bodies  is  given  by 

^  i  ^  II  {w,}^  (4  {wj  [2.13] 

^  i=l  ^  i=l 

where  u,  is  the  velocity  of  the  center  of  mass  of  the  t  th  body,  [7]^  is  the  inertia  tensor 
of  the  i  th  body,  and  {a;,}  is  the  angular  velocity  of  the  i  th  body.  The  total  kinetic 
energy  of  the  TROFM  is  equal  to  the  sum  of  its  four  components. 


T  =  Tm  +  Tj,  +  TT  +  TR  [2.14] 

In  order  to  apply  Lagrange’s  equation,  the  kinetic  energy  of  the  system  must  be 
expressed  in  inertial  coordinates. 

Consider  the  kinetic  energy  of  the  m«ist  Tm-  The  inertia  tensor  of  the  mast  is 
constant  when  observed  from  the  mast  coordinate  system  and  the  rrixy  and  nigz 
planes  are  planes  of  symmetry. 


I 


m 


H  ^ yym  7 yzm 
0  fyzm  Izzm 


|2.15| 
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Since  the  m  frame  does  not  translate,  the  kinetic  energy  of  the  mast  is 


r„  =  i  {>«“}’■  [Il„  [2.16! 

where  [I]^  is  the  inertia  tensor  of  the  mast  in  m  coordinates.  Since  m,  and  »,  are 
coincident  and  the  rotation  is  about  ii  only,  Equation  2.12  reduces  to 


1,  [2.17] 

Substituting  Equation  2.17  into  the  Equation  2.13  one  obtains  an  expression  for  the 
kinetic  energy  of  the  mast  in  generalized  inertial  coordinates. 

=  [2.18] 

Consider  the  kinetic  energy  of  the  pipe  body  Tp.  Since  the  pipe  is  fixed  at  the 
origin  of  the  pipe  body  axis  system,  the  kinetic  energy  of  the  pipe  body  is  given  by: 

r,  =  i  1%  [2.19] 

where  [I]^  is  the  inertia  tensor  of  the  pipe  body  as  observed  from  the  p  coordinate 
axis  system.  The  inertia  tensor  is  constant  when  observed  from  the  p  coordinate 
axis  system.  Since  the  pxy  and  Py^  planes  are  planes  of  symmetry,  the  inertia  tensor 
becomes 


[2.20] 


The  angular  velocities  are  transformed  to  inertial  coordinates.  Recall  that  the  pipe 


body  has  two  degrees  of  freedom,  0  and  thus,  from  Equation  2.12  one  obtains 


—  sin(0)'j 

{■u;P}  =  0  [2.21] 

cos(0)'P 
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Substituting  Equation  2.21  into  Ekjuation  2.19,  one  obtains  an  expression  for  the 
kinetic  energy  of  the  pipe  in  inertial  coordinates. 

Tp  =  ^/xxp sin2(0)4'^  -  /„p sin(0)4-" cos(0)  +  i/^p©^  +  co8*(0)4'2  [2.22] 

The  kinetic  energy  of  the  tilt  body  system  is  given  by 

T.  =  \  [I],  C*..'}  +  (2.23) 

where  Itcmcr  is  the  perpendicular  distance  between  the  center  of  maiss  of  the  tilt  body 
and  the  shaft  of  the  tilt  servo.  The  inertia  tensor  is  given  by 

fzxt  0  Iizt 

It  =  0  lyyt  0  [2.24] 

Izzt  0  lilt 

The  angulax  velocity  vector  {‘w*}  is  expressed  in  terms  of  its  Euler  angle  compo¬ 
nents.  Recall  that  the  tilt  body  has  three  degrees  of  freedom,  0,  and  thus, 
Equation  2.12  may  be  substituted  directly  into  Equation  2.23. 

Tt  =  ^/xxt  -  Izzt  sm{e)m  +  sin^(0)4'2  -  sin($)0  $  [2.25] 

Izzt  sin($)  sin(0)4'0  -|-  Int  cos(0)cos($)^  $  —  I^it  cos(0)  cos(4>)4'^  sin(0)  -t- 
]:Iyyt  cos^($)0^  -t-  lyyt  cos($)  cos(0)sin{$)4'0  -|-  ^lyyt  cos^(0)sin^($)4'^  + 

^Izzt  sin^($)0^  —  lilt  sin($)0  cos(0) cos($)^  +  ^hzt  cos^(0) cos^($)^* 

'  2”*  ^tcmcT^ 

Finally  the  kinetic  energy  of  the  rotor  is  found.  Recall  that  the  rotor  coordinate 
axis  system  is  fixed  to  the  center  of  mass  of  the  rotor  body  and  oriented  parallel  to 
the  tilt  body  axis  system.  This  means  that  the  rotor  body  axis  system  translates  as 
well  as  rotates  when  observed  from  the  inertial  frame.  The  rotor  body  is  assumed  to 
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be  axisymmetric  about  the  f,  axis  since  the  rotor  is  spinning  at  high  speed  relative 
to  the  other  rigid  bodies  in  the  system.  This  assumption  is  baised  on  actuator  disk 
theory  as  it  applies  to  propellers.  The  kinetic  energy  of  the  rotor  is  the  sum  of  its 
rotation<il  and  translational  components. 

r.  =  5  III,  +  im,  |{>v'}f  12.26) 


The  rotor  has  four  degrees  of  freedom,  ^,0,  and  /?.  The  angle  j3  is  the  angulair 
position  of  the  rotor  blade  with  respect  to  the  r  coordinate  system.  Since  the  rotor  is 
assumed  axisymmetric,  the  inertia  tensor  as  observed  from  the  r  coordinate  system 
is  constant  and  diagonal. 


Ir 


/xxr  0  0 

0  /„r  0 

0  0 


[2.27] 


The  rotor  spins  about  the  axis  with  an  angular  velocity  of  0.  The  rotor  axis  system 
also  rotates  with  respect  to  the  inertial  frame.  The  total  angular  velocity  of  the  rotor 
is  the  vector  sum  of  the  angular  velocity  of  the  tilt  body  components  and  the  rotor 
body  components  (see  page  12) 


1  0  —  sin(0) 

$ 

0 

0  cos($)  cos(0)sin($) 

0 

+ 

0 

0  —  sin($)  cos(0)cos($) 

Px 

0 

[2.281 


Substituting  Equation  2.28  into  the  first  term  of  Equation  2.26  one  obtains  an  ex¬ 
pression  for  the  rotational  kinetic  energy  of  the  rotor  body  in  inertial  coordinates. 


T,,«  =  s\n{Q)U  -h  sin2(0)'f2  -h  cos\^)Q^  +  [2.29] 

Ixxr  cos($)  cos(0)sin($)'i'0  +  cos^(0)  sin^($)'i'^  -f-  ^hzT  sin^($)0^  - 

IzzT  COs($)  cos(0)sin($)4’0  —  /„r  sin($)0/3  -f-  ^/zzr  cos^(0)cos^($)4'^  -h 
^zzr  cos(0)  COs($)'i'  0  +  ^/zzr  0^ 
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The  translational  kinetic  energy  of  the  rotor  body  is  given  by 


m,  (/,)^  0  m,  (/,)(/,) 

0  TTlr  (Ixf  0 

m,  (/,)(/,)  0  mrilrf 


-  Ircmcr)^ 


[2.30] 

where  Ircmcr  is  the  perpendicular  distance  between  the  center  of  mass  of  the  rotor 
and  the  shaft  of  the  tilt  angle  position  servo.  Substituting  Exjuation  2.12  into  Ekiua- 
tion  2.30  one  obtains  an  expression  for  the  translational  kinetic  energy  of  the  rotor 
body  in  inertial  coordinates. 


Tr,r.„  =  \rnr  (/.)'  -  m,  (/,)'  sin(0)4'$  +  (/,)*  sin*(0)4»2  +  [2.31] 

^mr  (Ir)^  cos^($)0^  +  im,  (Ix)^  cos^(0)  sin^($)4'*  + 

^rrir  (Ixf  sin^($)0^  +  (/x)^  cos^(0)  cos^($)4'^  + 

nir  (lx)  (Iz)  sin($)  sin(0)0^)  +  {lx)  (Iz)  cos(0)  cos($)'i'$  — 

(fx)  (Iz)  cos(0)  cos($)  sin(0)4’^  -  TTlr  (lx)  (/z)  sin($)0$ 


The  kinetic  energy  of  each  of  the  rigid  bodies  is  now  expressed  in  inertial  coordinates. 
Substituting  Equations  2.18,  2.22,  2.25,  2.29  and  2.31  into  Equation  2.14  one  obtains 
an  expression  for  the  kinetic  energy  of  the  entire  system  in  inertial  coordinates. 

T  =  +  ^/ripSin^(0)4'^  -  /xxpSin(0)’i'*COs(0)  +  ;^4yp0*  + 

^Ixxt  -  Ixxt  sin(0)4'$  +  ^/xxt  sin^(0)^^  -  Ixzt  sin($)0$  + 

Ixzt  cos(0)cos($)'j  $  —  Ixzt  cos(0)cos($)>p*  sin(0)  +  \:lyyt  cos*($)0^  + 
lyyt  cos($)  cos(0)  sin($)4'0  +  ^lyyt  cos^(0)  sin^(<&)^^  +  \:Izzt  sin^($)0^  — 
Izzt  sin($)0  cos(0)cos($)4'  +  cos*(0)  cos^($)^*  + 

\lxxr  -  Ixxr  sin(0)'j$  +  ^Ixxr  sin2(0)4'2  +  ^I^xr  COS^($)0'  +  [2.32] 

IxxT  cos($)  cos(0)sin(^>)^0  +  \:Ixxt  cos^(0)  sin^($)'P^  +  \Izzt  sin^($)0*  - 
Izzr  COs($)  cos(0)sin($)'i'0  -  Ixzt  sin($)0/3+  COS^(0)  COS*($)^^  + 
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Ixtr  cos(0)  cos($)4’  /j  +  (/,)*  -  m,  (/,)*  sin(©)4’$  + 

^TTir  Hxf  cos*($)0*  +  im,  (/*)*  cos^(0)  sin^($)\i'^  + 

^mr  (/rf  sin*($)0*  +  im,  (/,)’cos*(0)co8^($)«’  +  ^mW’cmcr^*  + 

»T»r  (/i)  (/»)  sin($)  sin(0)04'  +  m,  (/*)  (/,)  cos(0)  cos($)4'$  — 
rrir  (lx)  (/* )  cos(0)  cos($)  sin(0)4'*  —  m,  (/,)  (/*)  sin(<^)0$  + 
im,  (/,)*sin*(0)4'*  +  /x*<  sin($)  sin(0)4'©  +  cos*(0)®’ 

2.3  Lagrange’s  Equation 


Lagrange’s  equation  is  used  to  determine  the  equations  of  motion  from  the  kinetic 
energy  expression  of  Equation  2.32. 


dt  \dqj )  dqi  ‘ 


[2.33] 


where  q,-  is  the  i  th  generalized  coordinate  and  Q,  is  a  generalized  force.  In  this  case 
four  generalized  coordinates  are  used,  0,  ^  and  /?.  The  application  of  Equation  2.33 
results  in  four  equations  of  motion. 

US}-S-'- 

Evaluation  of  these  equations  yields  the  full  nonlinear  equations  of  motion  of  the 
system. 


2.4  Linearization  of  Equations  of  Motion 


The  kinematic  equations  which  result  from  the  evaluation  of  the  left  hand  side  of 
Lagrange’s  equation  are  linearized  in  two  steps.  First  the  following  assumptions  are 
made: 
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•9 


•  Perturbations  about  operating  point  are  small,  i.e.  A0,A®,A$  <  15®  =>■ 
sin(x)  =  x;  cos(x)  =  1. 

•  Products  of  small  angular  velocities  are  negligible,  i.e.  0'!'  w  0,  «  0,  etc.. 


•  Products  of  Inertia  such  as  and  4*  are  negligible. 

The  resulting  equations  of  motion  are 

F*  =  /a$  +  4«-/j(0  +  +  40^  [2.38] 

Fe  =  40  +  4^^  -  Izzrk^  -  04')  -  [2.39] 

F*  =  4^0  +  44'  +  I„rP  -  4zry3(00  -  <!»$)  [2.40] 

Fff  =  +  ^)  [2.41] 

where 


4  ~  (4rr  4“  4it  "i"  ^t^tcmcr  ^^(4ctncr  4)  "I"  ^r4) 

4  ”  (4yp  "h  4vt  "h  4xr  "f"  ^r4) 
4  “  (  4*t  "h  4xr  "I"  4y*  4zr) 
4  =  (24z»n  +  4zt  +  4zp  +  4zr  + 

4  TTlyl  ^  4if  4xt’ 


12.42] 


2.4.1  Generalized  Forces 

A  generalized  force  is  not  always  a  true  force.  If  qr,  is  a  translational  displacement, 
then  Qi  is  a  true  force;  however,  if  q,  is  an  angular  displacement  then  Qi  is  a  torque. 
In  this  case  all  generalized  forces  are  torques.  To  evaluate  each  generalized  force, 
one  imagines  that  one  of  the  coordinates,  g,  is  increased  by  a  small  amount  8qi  while 
all  other  coordinates  which  appear  in  the  kinetic  energy  expression  remain  fixed. 
Next,  an  expression  for  the  work,  done  by  all  driving  and  dissipative  forces  is 
formulated.  In  general  the  work  may  be  expressed  as 


5W„  =  F„  Sqi 
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[2.43] 


Equation  2.43  is  then  solved  for  the  generalized  force  Work  is  taken  positive  or 
negative  depending  on  whether  the  forces  tend  to  increase  or  decrease  qi.  The  method 
described  above  is  used  to  derive  expressions  for  F#,Fe,F*  and  Fp.  First  consider 
the  forces  acting  in  the  ^  direction.  There  is  a  lateral  force  exerted  by  the  tilted  rotor 
and  some  frictional  forces  in  the  bearing  at  the  base  of  the  mast  opposing  the  rotor 
force.  Equation  2.43  is  used  to  determine  the  generalized  force  F*  .  One  may  write 
an  expression  for  the  total  work  done  by  the  driving  and  dissipative  forces  through 
an  incremental  displacement  as 

SW’^  =  F*  [2-44]  < 

or 

SW^  =  (f,™  +  (2.451 

The  torque  produced  by  the  lateral  thrust  of  the  tilted  rotor  is  given  by 

F'irotcr  =  Fh  sin  $  cos  0  /x  [2.46] 

The  mast  bearing  exerts  three  types  of  frictional  torques,  static,  coulomb  and  viscous. 

The  viscous  bearing  torque  is  modeled  by  Equation  2.47. 

Ft,...™.  =  -('♦  |2-47| 

The  goal  of  this  development  is  to  derive  a  linear  model  of  the  system;  therefore, 
static  and  coulomb  friction  are  not  considered  since  these  phenomenon  are  highly 
nonlinear.  The  effects  of  such  nonlinear  phenomena  are  assessed  in  Chapter  3.  The  , 

generalized  force  F.^  immediately  follows  from  Equations  2.45,  2.46  and  2.47. 

F*  =  Fr  sin  $  cos  0  /x  —  4*  [2-48]  ^ 

where  p  is  the  viscous  friction  coefficient  of  the  yaw  bearing. 

Now  consider  the  generalized  force  Fq.  From  Equation  2.43  one  obtains 

6Wq  =  Fe  6Q  [2.49] 


20 


Note  that  Fq  will  be  similar'  to  F*  since  similar  driving  and  dissipative  forces  are  at 
work.  There  exists  a  driving  force  exerted  by  the  rotor  which  tends  to  increase  0. 
This  force  is  described  by  the  following  expression. 

Fe  =  FrIx  cos  $  [2.50] 

There  also  exists  static,  coulomb  and  viscous  frictional  torques  in  the  pipe  bearing 
similar  to  those  in  the  mast  beairing.  The  viscous  frictional  torques  are  modeled  by 
Equation  2.51. 

12.51) 

At  the  nominal  flight  condition  the  main  rotor  generates  thrust  to  overcome  gravi¬ 
tational  effects.  There  is  also  a  platform  located  approximately  1.5  rotor  diameters 
below  the  main  rotor  level  flight  condition.  The  platform  is  perpendicular  to  the  down- 
wcish  of  the  main  rotor.  This  configuration  results  in  a  ground  effect  phenomenon 
which  causes  the  fuselage  to  behave  as  though  there  is  a  torsional  spring  located  at 
the  pivot  point.  In  other  words,  for  a  given  main  rotor  angular  velocity,  the  thrust 
decreases  with  an  increase  in  altitude  (pitch  angle)  and  increaises  with  a  decrease  in 
altitude.  A  ground  effect  spring  constant  is  defined  such  that 

Fees  =  ©.  [2.52] 

The  gravitational  force  on  the  center  of  mass  of  the  entire  system  acting  over  a  length 
leg  Px  exerts  a  torque  which  tends  to  decrease  0. 

Fe,  =  —m  g  leg  cos  0  [2.53] 

Equation  2.43  can  now  be  written  as 

SWe  =  (FrIx  cos  ^  —  ly  Q  —  m  g  leg  cos  0  -  kesQ)  Se  [2.54] 

It  follows  that 

Fe  =  FrIx  cos  ^  —  u  Q  —  m  g  leg  cos  0  —  [2.55] 

The  generalized  forces  F^  and  Fp  are  evaluated  using  the  results  of  the  system  iden¬ 
tification  experiments  of  Chapter  3.  The  identification  experiments  have  shown  that 
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the  main  rotor  can  be  modeled  by  a  second  order  differential  equation  in  angular 
velocity  when  the  pitch  and  yaw  axes  are  restrained  in  their  nominal  positions. 

P=  -Tifi  -  np  +  TiPeom  [2.56] 

Identification  experiments  have  also  shown  that  a  second  order  differential  equation 
adequately  describes  the  tilt  servo  or  roll  axis  dynamics. 

$  =  +  ks^com  [2.57] 


Substituting  these  relations  into  Equations  2.38  and  2.41  for  the  test  conditions,  one 
may  obtain  differential  equations  for  the  generalized  forces  (torques)  Fp  and  F*. 

ni0c<»n  -/?)-/? 


Fp  =  Izzrfi  =  /* 


T2 


[2.58] 

[2.59] 


Ft  =  =  lai-ke^  -  ks^  +  ks^com) 

The  propeller  force  Fr  and  drag  torque  To  are  nonlinear  functions  of  main  rotor 
angular  velocity  which  may  be  linearized  about  a  nominal  angular  velocity  to  obtain 

[2.60] 

dTo, 


0(3 


To 


d(3 


[2.61] 


The  experiments  which  determine  the  values  of  the  gradients  and  constants  of  these 
equations  are  described  in  detail  in  Chapter  3. 

Substituting  the  generalized  force  expressions  into  Equation  2.38,  2.39,  2.40  and 
2.41  and  applying  the  aissumptions  of  t  Section  2.3  one  obtains  a  set  of  nonlinear 
coupled  differential  equations 


QFr  a  , 
j-(3  +  k2 


[d(3 


+  k,^,^)  =  lA  +  /«r/j(0  +  +  /.0'i'  [2.62] 

mgl,,  -  kcsQ  =  hQ  +  -  0^)  -  [2.63] 


QFr  o  ,  I. 


dp 


dP 

niPcom  -  h-  (3 


=  +  +  [2.64] 

=  /„,(^-$0  +  ^)  [2.65] 
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2.4.2  Linearization  using  Taylor  Series 


The  second  step  of  the  linearization  process  is  to  solve  Equations  2.62,  2.63,  2.64 
and  2.65  for  the  highest  order  derivatives  and  approximate  the  resulting  equations 
using  the  first  two  terms  of  a  Taylor  series  expansion.  Let  f(p,  <)  be  a  column  vector 
representing  the  highest  order  derivatives.  The  Taylor  series  expauision  becomes 


f(p,/)  =  f(p,t)  + 


(p-p)-h  HOT 


[2.66] 


where  p  =  [x^  u^j  is  the  augmented  state  and  control  vectors  and  p  contains  the 
nominal  values  of  p.  A  state  space  realization  follows  directly  from  this  step. 


where 


X  =  Ax  +  Bu 
y  =  Cx 

x= 

^  (/?com  ^com] 


[2.67] 


[2.68] 

[2.69] 


X  =  [0  0  0  0  —300  rad  sec  0  0  0]^ 
u  =  [—300  radfsec  0]^ 


di 

a* 

ai 

ai 

a* 

ai 

a<t 

a* 

a* 

ae 

a* 

ap 

ap 

a« 

ae 

a* 

ae 

ae 

ae 

ae 

ae 

ae 

ae 

ae 

a* 

ap 

dP 

a<t 

ae 

a* 

a* 

a* 

a4 

a* 

a» 

a* 

ai 

a* 

ae 

a# 

ap 

dP 

a* 

ae 

a* 

m 

ajl 

m 

m 

M 

M 

d* 

ae 

a* 

dP 

dP 

a^ 

ae 

a* 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

[2.70] 


[2.71] 


[2.72] 
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• 

94 .  ae  a*  0  0  0  0 

oT  _  a/Jcom  a/Jcom 

^  a*  ae  a*  -M_  0  0  0  0 

.  d^eotn  a*£<lfn  a^com  a*»n» 

The  non-zero  elements  of  the  state  space  model  are  given  below; 

[2.73] 

?i=-k. 

di 

(2.74! 

• 

i-<35l 

W 

1 

It 

«9'1-<D 

<Ci  1<SJ 

[2.75] 

• 

[2.76] 

dQ  '  Ud 

(2.771 

. 

nA  ” 

com 

[2.78] 

ae  /«ri3 
a$  h 

[2.79] 

a0  _  — 
h 

[2.80] 

ae 

a/3  fk 

(2.81] 

' 

a0  [^^  +  ^4]  (/<!  + ^c) 

a$  ~  hh 

[2.82] 

v| 

dQ  kcE 

^  "  /fc 

[2.83] 

_  ii 
a'f  ^d 

[2.84] 

24 

d0  Id 

_  ~Ic  +  *^2)  —  fnglcg^  +  /fc/r(^/?  +  kj) 

—  =  _ 


[2.85] 


[2.86] 


?i=-^ 

dft  Id 


[2.87] 


d  _T2H 

~  Id 


[2.88] 


d  (3  ^ 

dfi  ~  Id 

{!>>-  Ic-  Id)^T  +  {I c  +  Id)Tn gleg 
_  =  ^r2  — 


[2.89] 


[2.90] 


d  ji  _  T2(-/j-f  /„r) 

a/J  “  Id 

d^com 

2.5  Newtonian  Equations  of  Motion 


12.91) 


[2.92| 


A  simplified  Newtonian  approach  is  used  to  derive  the  TROFM  equations  of  mo¬ 
tion  to  verify  the  Lagrangian  results.  The  analysis  2issumed  that  pitch  and  yaw  axes 
are  uncoupled.  This  assumption  is  valid  for  small  deviations,  A0  <  15°,  A'P  <  15° 
about  the  operating  point,  0  =  <5  =  0. 

The  following  torques  and  forces  acted  on  the  system:  rotor  thrust,  rotor  drag 
torque,  viscous  friction,  ground  effect  ,  gyroscopic  couples  ,  and  gravitational  force. 
Nonlinear  friction  is  not  considered  here. 

Figure  2.6  shows  a  free  body  diagram  of  the  pitch  axis.  Applying  Newton’s  law 
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to  the  pitch  axis  one  obtains 


lyyQ  —  Fr  Ir  cos  $  —  mglcg  cos  0  -  uQ  —  kcE^  —  To  sin  $  +  (2.93) 

Where  Fr  and  To  are  experimentally  determined  to  be  linear  functions  of  the  main 
rotor  angular  velocity  about  the  operating  point  /?  in  Section  3.3.3. 


Fr  = 


it'**- 


(2.94) 


Td  = 


dTo 

[d^ 


+  ^4 


[2.95] 


Substituting  Equations  2.94  and  2.95  into  2.93  and  solving  for  0  one  obtains 


dFi 


0  =  - 


sin  $  + 


COS  $  — 


cos  0  — : - 


lyy  lyy  lyy  lyy  I  yy 


+  ■ 


lyy 

[2.96] 


Since  $  and  0  are  assumed  to  be  small  (<  15°),  one  may  write 


0 


a/j koE^  IzzT0^ 

*  r  r  r  T 


lyy  lyy  lyy  lyy  lyy  lyy 


[2.97] 


Figure  2.7  shows  a  free  body  diagram  of  the  yaw  axis.  Applying  Newton’s  laws, 
one  may  write 


/„'!>  =  FrIt  sin  cos  Q  +  To  cos  $  —  fi'i  —  hzrP  cos  $ 


[2.98] 


Solving  for  'I'  and  assuming  0  and  ^  <  15°  one  obtains 


't 


dFp 

aa 


/?  +  /C2 


Iz 


-lr^  + 


dTp 

90 


Izzr0 


Iz 


[2.99] 


Figure  2.8  shows  a  free  body  diagram  of  the  roll  axis.  Applying  Newton’s  laws 
and  the  results  of  the  tilt  servo  identification  of  Chapter  3  one  may  write 


ir  cos(0) 


X 


♦ 

z 


Figure  2.7  Free  body  diagram  of  the  yaw  axis. 


Figure  2.8  Free  body  diagram  of  the  roll  (tilt)  axis. 
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Solving  for  $  one  obtains 


4  =  +  (-*6$  -  +  iks^com)  -  7^/j(0  +  [2.101] 

ixxt 


Figure  2.9  shows  a  free  body  diagram  of  the  main  rotor  and  motor  armature. 
Applying  Newton’s  law  and  using  the  results  of  the  main  motor  identification  of 
Chapter  3  one  may  write 

-  $0)  =  -  /?)-  /?)  [2.102] 


Equations  2.97,  2.99,  2.101  and  2.102  are  solved  simultaneously  for  the  highest 
order  derivatives.  Substituting  the  expressions  for  the  highest  order  derivatives  into 
Equation  2.66  and  evaluating  the  first  two  terms  of  the  Taylor  series,  one  obtains  a 
linear  approximation  of  the  system  which  may  be  written  in  the  state  space  form  of 
Equation  2.67.  The  non-zero  elements  of  the  state  space  matrix  are  given  below: 


[2.103] 


[2.104] 


[2.105] 


_  -I„J 

dQ  Ixxt 

dQ  I„r0 

5$  lyy 

~  lyy 


[2.106] 

[2.107] 

[2.108] 
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de  _ 

dfi  In 

[2.109] 

de  _  -1/ 

de  lyy 

[2.110] 

de  —koE 
de  lyy 

[2.111] 

dfi  /« 

[2.112] 

d^  +  ^2]  It 

~  /„ 

[2.113] 

d^  —fx 

M  ~  hz 

[2.114] 

__  -hzr 

d(3  hz 

[2.115] 

dp  T2H 

d^  hz 

[2.116] 

dP  +  n 

dp  hz 

[2.117] 

—  hz)h  [^§fP  +  ^2]  +  hz  ^ghg 

I zz  lyy 

[2.118] 

d  P  _  Tiihzr  ~  hz) 

dP  ~  hz 

[2.119] 

dp 

dPcom 

[2.120] 
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2.6  Comparison  of  Linear  Lagrangian  and  Newtonian  Models 

Table  3.1  contains  the  numerical  values  of  the  physical  parameters  needed  to 
compute  the  elements  of  the  state  space  matrices.  Substituting  these  values  into  the 
results  obtzuned  in  Sections  2.4  and  2.5  one  obtains  the  theoretical  A  and  B  state 
space  matrices. 


^Lagrange  ““ 


-14.18  52.63 

-.6383  -1.062 


0  -153.42  1.5136  0 

-.07066  -2.138  -2.722  0 


-.7153  -.00202  -.01465  7.162 

160.627  -17.63  -654.94  -168.54 


0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

[2.121 


'Lagrange 


0  0  0  655.2  0  0  0  0 

153.2  0  0  0  0  0  0  0 


[2.122] 


-14.18  52.67 

-.6383  -1.062 


0  -153.42  2.032  0 

-.07066  -2.143  -2.722  0 


-.7153  -.00202  -.01465  7.168 

160.627  -17.63  -654.94  -168.64 


0  0 

0  0 

0  0 

0  0 

0  0 

[2.123] 
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^Newton  — 


T 


0  0  0  655.2  0  0  0  0 

153.2  0  0  0  0  0  0  0 


[2.124] 


One  can  see  that  the  two  methods  of  zmalysis  generate  nearly  identical  linear 


models  with  the  exception  of  the  An  term. 
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3.  SYSTEM  IDENTIFICATION 


3.1  System  Identification  Methods 

Accurate  knowledge  of  the  physical  parameters  which  comprise  the  elements  of 
the  state  space  model  of  Chapter  2  is  essential  to  the  development  of  linear  con¬ 
trol  systems.  Numerical  values  of  most  physical  parameters  could  not  be  obtained 
cinalytically;  therefore,  experimental  procedures  are  used  to  obtain  estimates  of  the 
parameters.  Parametric  (maximum  likelihood)  and  non-parametric  (frequency  re¬ 
sponse)  system  identihcation  methods  are  used  to  obtain  those  estimates. 

3.1.1  Continuous  Frequency  Response 

The  frequency  response  of  a  system  is  defined  as  the  steady  state  response  of 
a  system  to  a  sinusoidal  input  signal.  For  BIBO  (bounded  input  bounded  output) 
stable  linear  systems,  the  output  signal  will  be  sinusoidal  in  steady  state  but  different 
from  the  input  in  amplitude  and  phase.  The  frequency  response  function  describing 
the  sinusoidal  steady  state  can  be  obtained  by  replacing  s  with  ju  in  the  system 
transfer  function  T{s). 

=  TUu,)  [3.1) 

Recall  the  definition  of  the  Laplace  transform  of  a  continuous  signal  /(<)>  where  f{t) 
is  defined  for  /  >  0; 

£{/(()}  =  FW  =  /“/(‘K"*  [3.2] 

Jo 

thus; 

F(»)U  =  FO")  =  Hm)  =  J"  [3.31 
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which  is  defined  as  the  Fourier  transform  of  f{t).  Consider  ein  input  signal  whose 
Fourier  transform  is 

=  jT  (3.4) 

Applying  the  input  signal  to  a  linear  system  T[jui),  the  output  signal  is  given  by 

0{juj)  =  T{jui)I{ju).  [3.5] 

it  follows  that 

md  •  \  fo  4jl 

=  7^:0 

One  may  now  deduce  a  method  to  experimentally  determine  the  transfer  function 
T{s)  of  a  linear  system, 

•  Observe  the  input-output  time  histories; 

•  Evaluate  the  ratio  of  the  output  and  input  Fourier  Transforms; 

•  Plot  the  magnitude  and  phase  characteristics  and  determine  the  pole-zero  loca¬ 
tions  of  T{s)  using  Bode  analysis. 

3.1.2  Method  of  One  Frequency  at  a  Time  (Continuous) 

One  method  of  experimentally  determining  a  transfer  function  is  as  follows:  apply 
a  sinusoidal  input  of  amplitude  A  and  frequency  a;  to  the  system  of  interest; 

i{t)  =  A  sin  u)t  [3.7] 

wait  until  the  transients  have  died  out  and  observe  the  steady  state  system  response 
on  an  oscilloscope  to  determine  the  magnitude  and  phase  of  the  output. 

o{t)  =  Bsin(u;t  -f  <f>)  [3.8] 

Evaluating  the  ratio  of  the  output  to  input  one  can  deduce  the  frequency  response 
function  of  the  system  for  the  test  frequency  lj. 

TQu,)  =  |e>*  (3,91 
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Figure  3.1  Determining  amplitude  ratio  and  phase  angle  from  an  X-Y  plot  of 
output  sinusoid  versus  input  sinusoid. 

The  amplitude  ratio  BjA  and  the  phase  angle  4>  may  be  determined  by  plotting  o(t) 
versus  *(t)  using  the  X-Y  mode  of  an  oscilloscope.  The  equation  of  the  output  sinusoid 
may  be  written 

o(t)  =  B  sin{u}t)  cos{<l>) B  cos{u;t)  s\n{<i>)  [3.10] 

Note  that  when  i{t)  =  0,  sin(u;t)  =  0;  thus,  cos(u;t)  =  ±1  and  So  =  o(f)|i(i)=o  = 
Bs\n{<i>).  The  phase  angle  may  then  be  computed  from  Ek}uation  3.11. 

<f>  =  iarcsin  —  [3-ll] 

D 

Figure  3.1  shows  how  these  parameters  are  obtained  from  the  oscilloscope  display. 
The  sign  of  <{>  may  be  determined  by  noting  that  when  i{t)  goes  through  zero  (point 
q)  from  negative  to  positive  then  cos(u;t)  =  -|-1  [17].  This  means  that  if  the  curve 
advances  clockwise  in  time,  the  sign  of  is  positive  (phase  leawl).  Similarly  if  the 
curve  advances  counter-clockwise  in  time  the  sign  of  <f>  is  negative  (phase  lag).  The 
major  disadvantage  of  this  method  of  continuous  time  frequency  response  system 
identification  lies  its  inability  to  evaluate  the  frequency  response  function  for  more 
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than  one  frequency  at  a  time.  This  makes  the  continuous  time  method  slow  and 
tedious. 


3.1.3  Discrete  Frequency  Response 

Digital  computers  specially  equipped  for  real  time  data  acquisition  are  well  suited 
to  experimentally  determine  the  frequency  response  function  of  a  sampled  data  sys¬ 
tem.  Since  digital  computers  cannot  evaluate  a  continuous  time  Fourier  transform, 
one  must  define  a  Discrete  Fourier  Transform  (DFT)  which  can  be  calculated  from  a 
finite  record  of  sampled  data. 

Recall  the  definition  of  a  2-transform  of  a  discrete  sequence  of  sampled  data  f{kT) 
where  A:  €  0, 1, 2 ...  oo  and  T  is  the  sampling  period. 

£{/(»:))  2  =  (3.12) 

*=0 

It  can  be  shown  that  the  relationship  between  the  s  and  z  domain  is  given  by 

2  =  [3.13] 

Thus,  the  frequency  response  of  a  discrete  system  with  transfer  function  T{z)  = 
0{z)l I{z)  can  be  obtained  by 


[3.14] 


If  the  sequence  is  periodic,  i.e.  f{kT)  =  f{kT  +  NT)  then  the  Discrete  Fourier 
Transform  can  be  defined  as  the  finite  sum 


E/* 


=  DFTiM 


13.15) 


DFT(M  =  F,  =  E  A'" 

k=0 

The  DFT  has  the  following  properties; 


[3.16] 


•  It  is  discrete  in  frequency  since  it  is  defined  only  at  tj„  =  {2nn)/{NT) 


•  It  is  periodic  in  frequency  since 
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13.171 


which  follows  from  Euler’s  formula 

=  cos(^)  +  j  sin(5)  =  cos{S  +  2nl)  +  j  sin(d  +  27r/)  [3.18] 

where  /  is  an  integer. 

3.1.4  Method  of  One  Frequency  at  a  Time  (Discrete) 

One  method  of  evaluating  the  frequency  response  of  a  sampled  data  system  is 
analogous  to  the  continuous  time  method  of  one  frequency  at  a  time.  Apply  a  discrete 
time  sinusoid  with  zero  order  hold  to  the  system  to  be  identified. 

i{kT)  =  AsiniiJnkT)  [3.19] 

Record  the  steady  state  response  of  the  system  in  discrete  time  using  a  data  acquisi¬ 
tion  system. 

o{kT)  =  B  sin{u}nkT  +  4>)  [3.20] 

Evaluate  the  discrete  frequency  response  function  of  the  test  system  using 

the  following  facts: 

B/A  =  |T(e'"-’')| ,  I3.21J 

=  arg(r(e-'"-7)) 

The  amplitude  ratio  and  phase  angle  may  be  determined  by  plotting  o[kT)  versus 
i{kT)  on  the  computer  display  and  using  the  technique  described  in  Section  3.1.2.  The 
accuracy  of  the  estimated  transfer  function  is  affected  by:  the  number  of  samples  {N), 
noise,  transients  and  nonlinearities  in  the  test  system.  The  solution  to  these  problems 
is  to  compute  the  best  fit  sinusoid  to  the  output  by  estimating  B  and  which  best 
fit  o{kT)  [4].  Define  the  estimate  of  the  output: 

6{kT)  =  B  s\n{u}nkT  +  <i>n) 


B  sin(<^„)  cos{u;„kT)  -f 
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[3.22] 

[3.23] 


B  co8(0n)  8in(u;Mib7') 

BcCos(uj^kT)  +  B,nn{wnkT) 


Form  a  quadratic  cost  function; 


^  k=0 

=  jjE(o(kT)-B.co>(w.kT)- 


k=0 


B,  sin(u;„fc7'))^ 

Find  Be  and  B,  such  that  J  is  minimized, 


8  7  1 

=  jy  E2(»(<:r)-B,cosK*r) 


*=0 


BtS\n{uJnkT)){— COs{LiJnkT)  =  0 


87  1 

=  -5:2(<,(tr)-B.cos(u.,tr) 


k-o 


Bf  sin(uJnkT)){-  sin(w„fcr)  =  0 


For  Un  =  27rn/N 


2 

fie  =  o(f^'^)cos(2wnk/N) 

^  k=0 
2 

fi,  =  —  ^  o{kT)sin{2nnk/N) 
^  fe=0 


(3.24] 

[3.25] 


[3.26] 


[3.27] 


[3.28] 

[3.29] 


It  follows  that; 


thus, 


Be-jB,  = 

^  fc=0 

=  ^DF^oikT)) 

DFT{o(kT))=!^lB,-jB.] 


[3.30] 

[3.31] 


[3.32] 
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Figure  3.2  Typical  chirp  waveform 

Therefore  the  transfer  function  is  given  by 

^  DFTiojkT)) 

^  ^  DFT{i{kT)) 

One  may  conclude  that  the  DFT  automatically  computes  the  best  fit  sinusoid  to  the 
output  in  a  least  squares  sense. 

3.1.5  FVequency  Sweep  Method 

An  alternative  to  the  method  of  one  frequency  at  a  time  is  to  apply  a  sinusoidal 
input  whose  frequency  varies  from  low  to  high  over  the  range  of  interest.  Such  a 
signal  is  called  a  chirp.  Figure  3.2  shows  a  typical  chirp  waveform.  The  equation  of 
a  discrete  time  chirp  is  given  below: 

i{kT)  =  Ao  +  A  sin(u;tA:7')  [3.34] 

where 

k 

+  ~^{^end  ~  <^itaTt)  [3.35] 

The  frequency  sweep  system  identification  method  is  outlined  below. 

•  Determine  frequency  range  of  interest. 
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•  Select  sampling  rate  such  that  l/T  >  (l(kjen<<)/(2n').  This  follows  from  Shan¬ 
non’s  sampling  theorem. 

•  Choose  nominal  operating  point  Ao- 

•  Apply  the  chirp  input  to  the  system  T(z)  and  observe  the  output. 

•  Compute  the  DFT  of  the  input  and  output  sequences  using  Equation  3.16. 

•  Plot  magnitude  and  phase  of 

=  DFT{o{kT)) 

^  ’  DFT{i{kT)) 

versus 

2nn 

=  - . 

"  NT 

•  Estimate  the  pole- zero  locations  of  the  continuous  time  transfer  function  using 
Bode  analysis. 

3.1.6  Maximum  Likelihood  Parameter  Estimation 
Consider  a  system  of  linear  dynamic  equations: 

x(0)  =  xo  [3.38] 

x(t)  =  Ax(t)  +  Bu(t)  -I-  n(t)  [3.39] 

z{t,)  =  Cx{ti)  -h  Du(t,)  -I-  r]{ti)  [3.40] 

Equation  3.39  is  a  continuous-time  differential  equation  and  Equation  3.40  is  defined 
only  at  the  instants  t,.  This  is  representative  of  a  dynamic  system  with  data  sam¬ 
pled  by  a  digital  computer.  The  measurement  noise  rj  is  assumed  to  be  a  sequence 
of  independent,  zero  mean,  diagonal  covariance  Gaussian  vectors.  The  A,B,C,D 
matrices  are  functions  of  a  vector  of  unknown  parameters  An  estimate  of  the  true 
value  of  ^  based  on  the  measured  response  z(-)  to  an  input  u(')  is  computed  using 
an  implementation  of  a  maximum  likelihood  estimation  program  MMLE  described 
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in  Reference  [12].  The  general  form  of  a  maximum  likelihood  estimator  is  given  in 
Equation  3.41. 

f[z(-)]  =  argmaxL(z(-),^]  [3.41] 

where  L  is  the  likelihood  function.  For  the  system  described  above,  the  likelihood 
function  takes  the  form  of  a  multivariate  probability  distribution  function  [8]. 

=  l(2>-rdrtlGG-ll-'''’exp|i^(i<(i,)  -  2(ii)l-(GG-)-'[Z((<i)  -  *(<.)! 

13.42] 

where  z'{(t,)  is  the  Kalman  filter  estimate  of  z{ti)  based  on  measurements  up  to  f,_i. 
The  discrete  Kalman  filter  takes  the  form  of  a  current  state  estimator. 

x(f.)  =  +  r^[u{f.)  +  u(/,_i)]  [3.43] 

z^(f,)  =  Cx((<.)  +  Du(f,)  [3.44] 

X4(f.)  =  X((ti)  +  K[z(f.)  -  Z((<,)]  [3.45] 

where  $  =  F  =  e-^’’Bdr,  K  is  the  Kalman  gain  matrix  and  GG*  is  the 

covariance  matrix  of  the  residuals  z(ti)  —  One  can  see  that  the  estimate  Z4(ti) 

is  a  function  of  ^  because  the  Kalman  filter  depends  on  the  system  matrices. 

Initial  guesses  for  (  are  given  to  the  MMLE  program.  An  iterative  optimization 
routine  varies  the  values  of  ^  until  the  likelihood  function  is  maximized.  The  program 
actually  finds  ^  which  minimizes  the  less  complex  negative  log  likelihood  function  of 
Equation  3.46  which  is  equivalent  to  maximizing  Equation  3.42  since  In  Z.[z(-),^]  is  a 
monotonic. 

■>U)  =  (l^flM  -  2(fi)r(GG-)->|zf('.)  -  2('i)]]  +  5Af|n(det(GG-))  [3.46] 

The  MMLE  program  returns  the  maximum  likelihood  estimate  of  the  parameter 
vector 

3.1.7  Pseudo  Random  Binary  Sequences  (PRBS) 

Pseudo  Random  Binairy  Sequences  (PRBS)  or  maximum  length  sequences  (m- 
sequences)  generated  by  feedback  shift  registers  have  properties  of  a  good  input  test 


sign2d  for  system  identification  work  [3].  PRBS  test  signals  aure  generated  by  a  com¬ 
puter  implementation  of  Figure  3.3.  PRBS  are  periodic  with  period  P  =  (2"  —  1)7’, 
where  n  is  the  number  of  shift  registers  T  is  the  sampling  period.  The  implementation 
used  in  the  TROFM  identification  test  uses  13  shift  registers  with  a  sampling  period 
of  T  =  0.02245  sec.  The  minimum  time  the  signal  remained  in  a  positive  or  negative 
state  is  set  to  5T  =  .1123  sec  and  the  initial  state  of  eaich  shift  register  is  randomly 
set  before  each  test. 

3.2  Actuators 

The  main  rotor  and  tilt  servo  are  the  two  TROFM  actuators.  Their  principals  of 
operation  are  discussed  below. 

3.2.1  Main  Rotor 

The  the  thrust  produced  by  the  main  rotor  Fr  is  a  function  of  the  angular  velocity 
of  the  rotor  0.  The  angular  velocity  of  the  rotor  is  described  by  a  differential  equation. 
The  main  rotor  is  driven  by  a  DC  motor  and  its  angular  velocity  is  measured  by  a 
DC  tachometer.  The  commanded  rotor  angular  velocity  calculated  by  the  control 
computer  is  converted  to  a  analog  reference  voltage  by  a  D/A  converter  .  This 
reference  voltage  is  sent  to  a  pulse  width  modulated  amplifier  (PWM)  which  controlled 
the  motor  speed.  The  amplifier  compared  the  reference  voltage  generated  by  the 
computer  to  the  feedback  voltage  generated  by  the  tachometer  and  adjusted  the 
current  supplied  to  the  motor  to  decrease  the  angular  velocity  difference.  Figure  3.4 
shows  a  block  diagram  of  the  velocity  tracking  control  system  implemented  by  the 
PWM. 

3.2.2  Tilt  Servo 

The  actuator  which  tilted  the  main  rotor  is  a  rotary  position  servo  commonly 
found  in  remote  control  model  airplanes.  The  actual  servo  is  a  Futaba  9102  which 
produces  69.5  oz-in  of  torque  and  slews  at  273‘’/sec.  The  servo  position  is  controlled 
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Clock  Pulse 


Figure  3.3  Typical  PRBS  and  PRBS  signal  generator  consisting  of  digital 

components. 
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Pulse  Width  Modulated 

Position  Command  Position  Command  Angle 


Figure  3.5  Block  diagram  of  tilt  servo  angular  position  control  system. 

by  the  width  of  a  pulse  applied  to  the  servo's  control  terminal.  A  proportional 
control  system  is  built  into  the  servo  by  the  manufacturer  for  the  purpose  of  tracking 
position  commands.  An  external  potentiometer  configured  as  a  voltage  divider  is  used 
to  determine  the  actual  angular  position  of  the  output  shaft.  The  control  computer 
produces  analog  outputs  for  actuator  control,  but  the  servo  requires  that  position 
commands  be  in  the  form  of  pulses  of  varying  widths.  A  pulse  width  modulator  is 
built  to  convert  the  analog  output  of  the  computer  to  a  periodic  signal  of  square  waves 
whose  duty  cycle  is  proportional  to  the  voltage  produced  by  the  analog  output,  A 
schematic  of  this  circuit  is  included  in  Appendix  A.  Figure  3.5  shows  a  block  diagram 
of  the  servo  control  system  and  the  apparatus  used  to  measure  and  command  the 
angular  position  of  the  output  shaft. 

3.3  Identification  of  State  Space  Model 

The  form  of  the  state  space  model  is  determined  analytically  in  Chapter  2.  Several 
important  parameters  could  not  be  determined  analytically;  therefore  ,  experiments 
are  conducted  to  estimate  these  unknowns  as  well  as  correct  deficiencies  in  the  ana¬ 
lytical  model.  Figure  3.6  shows  the  system  identification  process  used  to  determine 
the  state  space  model  of  the  TROFM. 
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te  Limit  Noted 


Figure  3.6  Flow  chart  of  the  system  identification  process  used  to  determine  the 

state  space  model  of  the  TROFM. 


of  Unknown  Pi 


3.3.1  Main  Motor  (Rotor)  Identification 

Both  frequency  response  (non-parametric)  and  maximum  likelihood  (parametric) 
methods  are  used  to  identify  the  characteristics  of  the  main  motor.  A  chirp  input 
signal  is  used  to  excite  the  main  motor.  The  chirp  consists  of  2048  data  points  and 
sweeps  a  frequency  range  of  6.5  rad/sec  to  46  rad/sec.  The  nominal  rotor  speed  is  300 
rad/sec  and  the  chirp  amplitude  is  10.47  rad/sec.  The  roll  (tilt),  pitch  and  yaw  aoces 
are  restrained  at  $  =  0,'t  =  0  and  0  =  0  rad  during  this  test  to  insure  that  coupling 
between  the  axes  do  not  affect  the  results.  The  Bode  plot  generated  from  the  DFT 
suggests  that  the  system  was  second  order  because  of  the  -40  dB/dec  magnitude  slope 
starting  at  u;„  24  rad/sec  and  the  characteristics  of  the  phase  curve  (going  from 

0°  at  low  frequencies  to  -90"  at  24  rad/sec  and  to  -180°  at  40  rad/sec).  The  2  dB 
resonant  peak  indicates  that  the  system  has  a  damping  ratio  of  (^^oin  ^  0.35.  The 
transfer  function  ${$) / /3com{s)  is  assumed  to  be  of  the  form: 


/3com('S)  d”  ^Cmain^n  moin  d"  i 


13.47) 


A  maximum  likelihood  analysis  is  then  performed  to  determine  the  best  estimates 
of  Wnmoin  and  Cmoin-  The  estimates  from  the  DFT  results  are  used  as  initial  guesses 
to  start  the  numerical  search.  The  zero  order  hold  on  the  input  signal  introduces 
approximately  a  half  sampling  period  delay  [12].  This  was  is  corrected  in  the  analysis 
by  replacing  the  input  with  Equation  3.48. 


jW  —  20com{k  d-  1)  +  0com(k)) 


|3.48) 


The  maximum  likelihood  estimate  of  the  main  motor  transfer  function  is  given  in 
Eiquation  3.49. 

= _  13491 

17.62s +  655.2  '■  ' 

It  follows  that  ijJn  main  =  25.597  rad/sec  and  Cmain  =  0.3442.  This  compares  favorably 

to  the  initial  estimate  obtained  from  the  DFT  results.  Figure  3.7  shows  the  simulated 

and  actual  (rate  limit  induced  DC  offset  removed)  responses  of  the  output  0  to  the 
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Figure  3.8  Main  motor  frequency  response. 
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chirp  0com  input  signal.  Figure  3.8  shows  the  Bode  plots  of  the  maximum  likelihood 
and  DFT  transfer  function  estimates. 

An  unusual  trend  is  observed  in  the  actual  output  time  history.  At  frequencies 
greater  than  25  rad/sec,  the  positive  peak  values  of  /?(<)  decrease  as  the  frequency 
of  /3com(0  increases.  The  negative  peaks  of  the  output  increase  to  about  2/?com(0 
the  resonant  frequency  and  slowly  decrease  at  frequencies  higher  than  the  resonant 
frequency.  This  non-linear  phenomenon  has  the  effect  of  lowering  the  nominal  value  of 
/3(t)  at  frequencies  greater  than  25  rad/sec  by  6  to  7  rad/sec.  The  time  histories  used 
to  identify  the  linear  model  are  based  on  this  nonlinear  data  with  the  effective  high 
frequency  DC  offset  removed.  The  linear  model  follows  the  nonlinear  time  history  up 
to  the  resonant  peak  at  25  rad/sec  and  follows  amplitude  well  throughout  the  test 
frequency  range  with  a  DC  offset  of  6  to  7  rad/sec  for  test  frequencies  greater  than 
25  rad /sec. 

Close  examination  of  the  time  history  allows  one  to  ciscertain  the  type  of  nonlinew 
phenomenon  at  work.  The  slope  of  /3(f)  going  up  is  less  than  the  slope  of  /3(f) 
going  down  which  indicates  that  a  rate  limit  is  present.  The  /3(f)  time  history  is 
differentiated  according  to  Equation  3.50. 

m  =  |3,50| 

It  was  found  that  P{t)  is  limited  to  about  180  rad/sec^  in  the  positive  direction 
(acceleration).  Rate  limiting  is  not  factor  in  the  negative  direction  (deceleration). 
In  order  to  simulate  the  nonlinear  main  motor,  the  structure  shown  in  Figure  3.9  is 
assumed.  The  structure  of  the  linear  model  is  shown  for  comparison  in  Figure  3.10. 
The  linear  model  values  a  =  17.67  and  K  =  30.08  are  used  to  start  an  iterative  search 
for  a  and  K  which  best  fit  the  nonlinear  time  history  data  in  a  visual  sense.  Values 
of  a  and  K  which  fit  the  nonlinear  time  history  well  are  found  to  be  a  =  23  and 
K  =  100.  Figure  3.11  shows  the  actual  and  nonlinear  simulated  response  of  the  main 
rotor  speed  to  a  chirp  input.  One  can  see  that  the  nonlinear  simulation  models  the 
actual  response  quite  well.  Chapter  3  will  show  that  the  linear  model  is  suitable  for 
pitch  and  yaw  axis  control  design  since  the  closed  loop  bandwidths  of  the  command 
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PWM  Controller 


Figure  3.9  Block  diagram  of  rate  limited  main  motor  model. 

systems  axe  at  least  a  factor  of  2  below  the  frequency  where  the  nonlinear  rate  limit 
effects  appear. 

The  rate  limit  described  above  is  actually  an  acceleration  limit  because  main  motor 
speed  is  rate  limited.  The  motor  is  acceleration  limited  in  the  positive  direction 
for  two  reasons:  drag  torque  and  limited  power.  At  the  nominal  /?,  the  motor  is 
drawing  power  from  the  amplifier  to  overcome  drag  torque.  To  accelerate  in  the 
opposite  direction  of  the  nominal  drag  torque  vector  the  motor  requires  more  power 
from  the  amplifier  to  overcome  an  increase  in  drag  torque  as  well  as  inertial  torques. 
Since  the  amplifier  is  already  using  significant  power  to  maintain  the  nominal  $,  it 
is  unable  to  increase  the  rotor  speed  as  quickly  in  the  direction  of  the  nominal  drag 
torque  vector.  Decreasing  rotor  speed  from  the  nominal  (3  can  be  accomplished  by 
simply  removing  the  power,  since  drag  torque  and  friction  paissively  cause  the  rotor  to 
decelerate.  The  amplifier  is  also  capable  of  driving  the  rotor  in  reverse  which  means 
that  at  the  nominal  /?,  it  has  much  more  power  available  to  decelerate  than  accelerate. 
Figure  3.12  graphically  illustrates  these  principles.  The  limited  power  also  introduces 
another  important  nonlinear  effect,  0  saturation.  The  amplifier  eventually  runs  out 
of  power  to  accelerate  the  rotor  since  all  of  the  power  is  used  to  overcome  drag  torque. 
This  critical  speed  was  experimentally  determined  to  be  /?mar  =  362  rad/sec.  The 
saturation  nonlinearity  is  considered  in  the  control  system  design. 


53 


PWM  Controller 

Gain  Motor 


Figure  3.10  Block  diagram  of  linear  main  motor  model. 


3.3.2  Tilt  Servo 


Frequency  response  and  maximum  likelihood  methods  are  used  to  identify  the 
characteristics  of  the  tilt  servo.  A  chirp  input  signal  is  applied  to  the  control  terminal 
of  the  servo.  The  chirp  consists  of  2048  data  points  with  a  frequency  range  of  2.75 
to  30.5  rad/sec.  The  output  shait  position  is  sampled  at  50.84  Hz  or  319.28  rad/sec. 
The  pitch  and  yaw  axes  are  restrained  at  0  =  0  and  =  0  rad  while  the  main 
rotor  maintains  its  nominal  speed  of  300  rad/sec.  The  Bode  plot  generated  from  the 
DFT  suggests  that  the  system  was  second  order  because  of  the  -40  dB/dec  slope  at 
frequencies  greater  than  a;„  tut  «  12  rad/sec  aP.d  the  charaicteristics  of  the  phase  curve 
(going  from  0°  at  steady  state  to  -90°  at  u;  =  12  rad/sec  to  180°  at  25  riid/sec).  The 
shape  of  the  Bode  magnitude  curve  in  the  region  of  u;„  ti/t  ^  12rad/sec  indicates  that 
Ctiit  *2  0.6.  The  transfer  function  $(s)/$<^(s)  is  assumed  to  be  of  the  form 


^ _ u;l  _ 

^com('S)  T  ^^tiltUJn  tilt  "b  tilt 


[3.51] 


A  miiximum  likelihood  analysis  is  performed  to  determine  the  best  estimates  of  u;„  tUt 


and  The  estimates  from  the  DFT  results  are  used  as  initial  guesses  to  start 
the  numerical  search.  A  zero  order  hold  is  used  on  the  input  signal;  therefore,  in  an 
approximate  sense,  a  half  sampling  period  delay  was  introduced.  To  compensate  for 
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Continued 


Figure  3.11  Actual  and  nonlinear  simulated  response  of  main  rotor  speed  to  a  chirp 

input. 
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this  time  delay,  the  input  is  repla^red  by  Equation  3.52. 

+  1)  -  *c^(k)  p 

The  maximum  likelihood  estimate  of  the  tilt  servo  transfer  function  is  given  by  Ex]ua- 
tion  3.53. 

*(»)  153.42 

s=+ 14. 18i+  153.42  '■  ' 

It  follows  that  Wntiit  =  12.39  rad/sec  and  Qiit  =  0.5724.  This  compares  favorably  to 

the  initial  estimate  obtained  from  the  DFT  results.  Figure  3.13  shows  the  simulated 
and  actual  time  history  responses  of  the  output  $  to  the  chirp  ^com  input  signal.  Fig¬ 
ure  3.14  shows  the  Bode  plots  of  the  maximum  likelihood  and  DFT  transfer  function 
estimates. 

3.3.3  Drag  Torque,  Thrust  and  Rotor  Speed 

The  relationship  between  rotor  thrust  and  rotor  velocity  is  determined  experimen¬ 
tally.  The  apparatus  shown  in  Figure  3.16  is  used  to  measure  the  rotor  force  as  the 
rotor  speed  is  slowly  varied  from  0  to  362  rad/sec.  The  rotor  speed  was  measured  by 
a  DC  tachometer  mounted  on  the  motor.  A  second  order  polynomial  is  fitted  to  the 
experimental  measurement  using  a  least  squtu-es  procedure.  The  curve  is  lineau'ized 
about  the  nominal  rotor  speed  j3 =300  rad/sec.  This  linearized  relation  is  given  in 
Equation  3.54  and  shown  in  Figure  3.15. 

SFr  ■  Ibf  • 

+  =  -0.003057 — -  0.59/6  f  [3.54] 

5/?  rad/ sec 

An  experiment  was  designed  for  the  purpose  of  determining  the  rotor  drag  torque  as 
a  function  of  main  rotor  angular  velocity  (/?).  The  pitch  axis  is  restrained  at  0  =  0 
during  these  experiments.  Figure  3.17  shows  a  free  body  diagram  of  the  system. 
Summing  torques  about  the  axis  one  obtains; 

^  A/ij  =  Full  sin  $  —  Fsch  ~  Fsf  —  Tq  =  0  [3.55] 

where  Fr,  Fsgi  Tsf,  and  To  denote  the  force  generated  by  the  main  rotor,  the  force 
measured  by  the  strain  gage  force  transducer,  yaw  bearing  static  friction  and  the 
rotor  drag  torque  respectively. 
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Figure  3.13  Response  of  tilt  servo  to  chirp  input. 
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Figure  3.15  Rotor  force  Fr  as  a  function  of  main  rotor  speed  l^j. 
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Figure  3.16  Sketch  of  apparatus  used  to  measure  forces. 
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Top  View 
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Side  View  at  Main  Rotor 


Figure  3.17  Free  body  diagram  of  system 


The  main  rotor  is  tilted  to  =  20"  and  0  is  varied  from  0  to  362  rad/sec.  The 
force  is  measured  by  the  computer  and  plotted  as  a  function  of  rotor  angular  velocity. 
The  main  rotor  is  then  tilted  to  $2  =  —20"  and  the  process  repeated.  If  the  magnitude 
of  Tsf  is  the  same  for  positive  and  negative  $  and  there  is  no  drag  torque  generated 
by  the  rotor  then  the  force  versus  velocity  curves  would  be  identical;  however,  drag 
torque  exists  and  the  static  friction  of  the  yaw  bearing  is  direction  dependent.  Two 
moment  equilibrium  equations  may  now  be  written. 

Afti,  =  Ffl/i  sin  $1  —  Fscxh  —  TsFi  —  Fd  cos($)  =  0  [3.56] 

52  A/2tz  =  FrIi  sin  ^2  -  FsGjh  -  TsFi  -  To  cos($)  =  0  [3.57] 


Since  $1  =  —  $2,  one  may  add  these  equations  and  solve  for  the  drag  torque. 

rp  _  (FsGa  -  FsGi  )l2  +  {TsFj  -  TsFx  ) 

^  2  cos 

The  static  friction  torques  are  determined  using  the  apparatus  shown  in  Figure  3.16. 

A  force  is  applied  to  the  end  of  the  pipe  which  causes  the  force  transducer  to  deflect. 
The  force  is  slowly  removed  until  the  spring  force  of  the  force  transducer  acting  over 
I2  and  static  friction  torque  balance. 

Tsf  =  Fsf  h  [3.59] 

The  balancing  force  Fsf  is  measured  in  both  the  positive  and  negative  *5  directions. 
It  has  been  determined  that  Fsf{+'^)  ^  1.30oz  and  Fsf{  —  '^)  ^  108  oz.  Figure  3.18 
shows  FsGt  ^  *  function  of  0  for  $1  =  20"  as  measured  by  the  computer.  A  second 
order  polynomial  least  squares  fit  is  applied  to  the  data  to  facilitate  the  implementa¬ 
tion  of  Equation  3.58.  Figure  3.19  shows  Fsg2  ^  *  function  of  main  rotor  angular 
velocity  for  $2  =  —20°.  By  applying  Equation  3.58  and  using  the  least  squares  fits 
in  Figures  3.18  and  3.19  one  can  obtain  To{0)-  Figure  3.20  shows  the  relationship 
between  To  and  0.  The  drag  torque  curve  of  Figure  3.20  is  linearized  about  the 
nominal  rotor  speed  of  /3=300  rad /sec.  This  linear  relationship  is  found  to  be 

Td  »  ^0  +  A:4  =  -0.000496;g  -  0.08 ft.lbf.  [3.60] 

d0  radfsec 
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Figure  3.18  Lateral  force  Fsai  =  20®  as  a  function  of  /?. 


Figure  3.19  Lateral  force  Fsa^  for  ^2  =  -20®  as  a  function  of  0. 
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Figure  3.20 


Main  Rotor  Drag  torque  as  a  function  of 
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3.3.4  Pitch  to  Rotor  Speed  TVansfer  f\inction 


The  pitch  angle  to  rotor  speed  transfer  function  0(s)//j(s)  is  to  be  determined. 
The  roll  and  yaw  axes  are  restrained  at  $  =  0  and  ^  =  0  during  this  test  to  insure 
inter-axis  coupling  does  not  affect  the  results.  A  chirp  ficom  is  issued  and  the  outputs 
$  and  0  are  measured.  The  input  chirp  consists  of  4096  points  sampled  at  49.54 
Hz  and  covers  the  frequency  range  between  .05  and  3.2  rad/sec.  Frequency  response 
and  MMLE  methods  axe  used  to  estimate  the  transfer  function  Q{s)/0{s).  One  At 
a  Time  (OAT)  frequency  response  tests  are  also  conducted  to  verify  the  chirp  test 
results.  The  pitch  axis  is  extremely  difficult  to  identify  because  numerous  nonlinear 
effects  are  present.  The  ground  effect  phenomenon  introduces  a  nonlinear  spring.  For 
a  fixed  rotor  speed  and  attitude  0  where  the  TROFM  is  in  ground  effect,  a  decrease 
in  0  results  in  a  ground  effect  force  which  tends  to  restore  0  to  its  nominal  value. 
Similarly  an  increase  in  0  results  in  a  loss  of  ground  effect  force  which  tends  to  restore 
0  to  its  nominal  value.  These  ground  effect  restoring  forces  Fqe  are  nonlinear  because 
dFcsIdQ  =  /(0)  [7].  For  large  nominal  values  of  0  (large  distances  between  platform 
and  rotor),  the  gradient  dFoE/^Q  =  0.  In  other  words  the  TROFM  is  out  of  the 
ground  effect.  Other  sources  of  nonlinearity  include  static  and  coulomb  friction  which 
are  the  probable  cause  of  the  hysteresis  and  beat  phenomena  which  are  observed  in 
some  of  the  OAT  frequency  response  tests. 

Despite  the  nonlinearities,  a  linear  model  is  found  which  adequately  describes  the 
system  for  linear  control  design  purposes.  The  transfer  function  0(s)//j(s)  is  assumed 
to  be  of  the  form  of  Equation  3.61  which  follows  from  the  analysis  of  Chapter  2. 

ew  ^  ^  ^  p 

/3{s)  5  +  ^Cpttch^n  pitch  pitch 

The  maximum  likelihood  estimate  of  this  transfer  function  is  given  by  Equation  3.62. 


0(s)  0.0359 

0{s)  ~  s2  +  1.064s +  2.7225 


[3.62] 


The  MMLE  estimate  nicely  matches  the  low  frequency  portion  of  the  DFT  generated 
Bode  plot  but  underestimates  the  resonant  peak  and  break  frequency.  Based  on  the 
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results  of  the  OAT  frequency  response  tests,  a  modification  to  the  MMLE  estimate 
of  dFiildfi  is  made.  The  OAT  tests  indicate  that  the  DPT  chirp  underestimates  the 
low  frequency  magnitude  curve  by  about  5  to  10  dB;  therefore,  the  MMLE  result  is 
modified  to  Equation  3.63. 

e(s)  _  0.05  .  , 

/3(s)  ~  +  1.064s  +  2.7225  ‘  ' 

OAT  results  indicate  that  the  MMLE  result  made  a  good  estimate  of  the  pole  lo¬ 
cations.  The  OAT  tests  indicate  a  2-4  dB  droop  in  the  magnitude  curve  between 
u;  =  0.8  2Lnd  1.2  rad/sec  which  suggests  the  existence  of  a  third  pole  around  u  =  0.7 
rad/sec.  Assuming  the  pitch  axis  is  a  3rd  order  system,  the  approximate  -40  dB/dec 
slope  at  frequencies  above  2  rad/sec  suggests  the  existence  of  a  first  order  zero  around 
the  2  rad/sec  region.  No  physical  explanation  as  to  why  this  would  be  a  third  order 
system  was  found. 

The  modified  MMLE  transfer  function  estimate  was  found  to  be  an  acceptable 
low  order  linear  model  of  a  higher  order  nonlinear  system.  DFT  and  OAT  frequency 
responses  as  well  as  the  modified  MMLE  frequency  response  are  shown  in  Figure  3.21. 
The  actual  and  simulated  chirp  response  is  shown  in  Figure  3.22.  One  can  see  that 
the  simulated  response  is  in  phase  with  the  actual  response  throughout  the  test.  The 
amplitude  follows  well  in  the  very  low  and  high  frequency  ranges.  In  the  middle 
of  the  response  (mid  frequency  range)  the  MMLE  response  over-estimates  the  chirp 
amplitude.  This  is  further  evidence  of  a  possible  third  pole  around  0.7  rad /sec  and  a 
zero  around  2  rad/sec. 

3.3.5  Yaw  Axis  Friction  Coefficient 

In  this  test,  the  pitch  axis  is  restrained  at  0  =  0  and  the  tilt  angle  is  set  such 
that  it  cancels  the  drag  torque  generated  by  the  rotor  spinning  at  300  rad/sec.  The 
tilt  angle  is  increased  (decreased)  by  0.07  rad  (-0.07  rad)  and  the  steady  state  yaw 
angular  velocity  is  determined  from  the  slope  of  a  yaw  attitude  versus  time  plot.  The 
moment  equilibrium  equation  below  is  used  to  compute  the  yaw  axis  viscous  friction 
coefficient  /i. 
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Figure  3.21  Frequecy  response  of  0(s)//3  determined  by  OAT,  DFT  and  maximum 

likelihood  methods. 
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[3.64] 

[3.65] 


5;;a/„  =  7h+^4»  =  o 


The  test  is  repeated  several  times  in  both  directions.  The  average  clockwise  and 
counterclockwise  yaw  axis  viscous  friction  coefficients  are  found  to  be  /icw  =  0.01983 
(ft-lbf)  /  (rad/sec)  and  ficcw  =  0.02858  (ft-lbf)  /  (rad/sec)  respectively.  This  means 
that  the  yaw  axis  equation  of  motion  is  nonlinear  since  the  viscous  damping  coefficient 
is  direction  dependant.  For  the  purpose  of  linear  control  design,  the  CW  and  CCW 
viscous  damping  coefficients  are  averaged.  An  estimate  of  the  ^(s)/$(s)  transfer 
function  based  on  the  analysis  of  Chapter  1  and  the  results  of  this  section  is  given  in 
Equation  3.66. 


<^(s)  s{s  +  s{s  +  0.71526) 


[3.66] 


3.3.6  Gyroscopic  Couples 


The  elements  of  the  state  space  model  associated  with  gyroscopic  couples  were 
dQ/d^  and  d^/dQ.  In  this  experiment,  two  independently  generated  pseudo-random 
binary  sequences  with  amplitudes  of  15  rad/«ec  and  0.1  rad  are  simultaneously  applied 
as  inputs  to  the  main  rotor  and  tilt  servo  respectively.  The  input-output  time  histories 
are  sampled  at  44.54  Hz.  The  MMLE  program  is  used  to  estimate  the  values  of  the 
gyrocopic  stablity  derivatives.  The  results  obtained  from  the  MMLE  program  are 
given  below. 

d@ 

ss  -0.7030  sec  [3.67] 

— as  6.679  sec 

ae 


Actual  and  simulated  responses  of  the  system  to  a  set  of  PRBS  inputs  can  be  found 
in  Section  3.4. 
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3.3.7  Back  Substitution  of  Results  into  Analytical  Model 


The  following  stability  derivatives  are  computed  by  back-substituting  the  results 
of  the  system  identification  into  the  analytical  model  developed  in  Chapter  1. 


a/3  Tj/i 


=  160.6269  sec' 


a/3  , 

=  -168.5407  sec-2 

a4> 


^  ^  =  -0.0146s«- 

a/3  Id 


[3.68] 


[3.69] 

[3.70] 


Maximum  likelihood  estimation  of  the  stability  derivatives  d^/dQ  and  dQ/d^  in¬ 
dicate  that  the  model  is  insensitive  to  these  parameters.  Since  the  estimates  were 
small,  they  were  set  to  zero. 

3.4  Results 

Table  3.1  lists  the  measured  and  identified  physical  parameters  necessary  to  obtain 
the  state  space  model  derived  in  Chapter  2.  The  moments  of  inertia  are  determined 
by  measuring  and  weighing  the  components  of  the  TROFM  and  performing  weight 
and  balance  calculations. 

The  state  space  model  used  for  control  design  is  given  below.  Some  parameter 
estimates  in  this  model  are  based  on  earlier  system  identification  results.  A  refined 
linear  model  was  developed  as  time  progressed. 


-14.18  6.679 
-.638  1.064 


-153.42 


-2.7225  0 


0  -0.715260  -.002  -.0146  7.16 

0  160.63  -17.63  -655.2  -168.54 


|3.71| 


Table  3.1  List  of  important  physical  parameters. 


B  = 


0  0  0  655.2  0  0  0  0 

153.42  0  0  0  0  0  0  0 


T 


C  = 


0  0  0  0  1  0  0  0 
0  0  0  0  0  1  0  0 
0  0  0  0  0  0  1  0 
0  0  0  0  0  0  0  1 


Figure  3.23  shows  the  actual  and  simulated  response  of  this  system  to  two  inde¬ 
pendently  generated  PRBS  ^com  and  input  signals  where  T  =  0.05  sec  and 
Pcom  =  ±10  rad /sec  and  =  ±0.1  rad.  One  can  see  that  the  linear  model  follows 

the  actual  responses  reasonably  well  except  for  the  yaw  axis  response.  The  linear 
model  describes  the  short  term  response  of  the  yaw  axis  quite  well  but  diverges  from 
the  actual  response  as  time  progresses.  The  linear  yaw  axis  model  has  a  pure  integral 
term  which  explains  the  divergence  of  the  simulated  yaw  response.  In  reality  static 
friction  is  present;  therefore,  each  time  the  pitch  rate  equals  zero,  the  static  friction 
must  be  overcome  by  a  control  input.  This  explains  why  the  actual  response  does  not 
diverge.  Empirically,  the  effect  of  static  friction  in  the  yaw  axis  can  be  described  by 
replacing  the  pure  integral  term  with  a  low  frequecy  lag  with  a  pole  location  around 
s  w  —  .7153.  A  refined  system  model  is  given  in  Equation  3.72.  This  model  is  the  re¬ 
sult  of  futher  system  identification  work  inspired  by  physical  insights  obtained  during 
the  implementation  of  the  control  designs  of  Chapter  4.  Figure  3.24  shows  the  actual 
and  simulated  response  of  this  system  to  two  independently  generated  PRBS  0com  and 
^com  input  signals  where  T  =  0.05  sec  and  ^com  =  ±10  reid/sec  and  ^com  =  ±0.1  rad. 
One  can  see  that  the  refined  linear  model  follows  all  of  actual  responses  reasonably 
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Figure  3.23 


Response  of  control  design  model  to  PRBS  inputs. 
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well.  It  is  recommended  that  future  researchers  start  their  work  using  this  model. 
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[3.72] 


0  0  0  0  1  0  0  0 
0  0  0  0  0  1  0  0 
0  0  0  0  0  0  1  0 
0  0  0  0  0  0  0  1 
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Figure  3.24  Response  of  refined  model  and  actual  system  to  PRBS  inputs. 
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4.  CONTROL  SYSTEM  DESIGN 


Three  types  of  control  systems  are  designed  for  the  purpose  of  determining  their 
effects  upon  pilot  performance.  An  open  loop  uncompensated  system  is  designed  such 
that  a  joystick  commands  main  rotor  speed  for  pitch  axis  control  and  rotor  tilt  angle 
for  yaw  axis  control.  A  rate  command  system  is  designed  where  the  joystick  com¬ 
mands  pitch  and  yaw  rate.  This  is  accomplished  by  using  a  feedback  control  system 
with  lead-lag  compensators.  Pitch  and  yaw  angle  sensor  signals  can  not  be  differenti¬ 
ated  to  obtain  reliable  rate  feedback  signals  because  of  sensor  noise.  A  Kalman  filter 
is  designed  obtain  reliable  estimates  of  all  unmeasurable  states,  particularly  pitch 
and  yaw  rate  (0  and  'P),  in  spite  of  noisy  measurements,  unknown  disturbances  and 
model  uncertainty.  Attitude  command  systems  are  designed  such  that  the  joystick 
commands  pitch  and  yaw  attitude.  This  is  accomplished  by  feeding  back  attitude 
measurements  from  the  pitch  and  yaw  angle  sensors  to  the  main  rotor  and  tilt  servo 
respectively. 

Compensators  for  the  rate  and  attitude  command  systems  are  designed  in  con¬ 
tinuous  time  assuming  measurements  are  available  for  feedback.  Discrete  equivalents 
of  the  continuous  compensators  are  then  designed  by  emulation  using  a  frequency 
prewarped  bilinear  transformation.  The  command  systems  are  implemented  digitally 
on  a  Mac  Il-si  digital  computer  equipped  with  an  A/D,  D/A  board. 

4.1  Discrete  Models  of  Sampled  Data  Systems 

The  systems  developed  in  this  chapter  contain  both  continuous  and  discrete  com¬ 
ponents.  Compensators  are  designed  assuming  that  continuous  signals  are  available 
for  feedback.  Digital  compensators  (filters)  are  designed  to  emulate  their  continuous 


77 


is  equivalent  to  : 


ZOH 


G{s) 


]y_o(kT) 


“V 

G(z) 


Figure  4.1  General  sampled  data  system. 


counterparts  because  a  digital  computer  with  A/D  and  D/.A  converters  is  used  to  im¬ 
plement  control  laws.  The  A/D  converter  samples  continuous  time  signals  at  regular 
time  intervals  T.  The  computer  then  processes  the  samples  through  the  control  law 
(digital  filter)  and  sends  the  command  to  a  D/A  converter  which  holds  the  command 
constant  until  the  next  sample  is  received.  The  D/A  converter  acts  as  a  zero  order 
hold  (ZOH)  which  derives  its  name  from  the  order  of  the  polynomial  function  which 
describes  the  shape  of  the  output  between  samples.  Since  the  output  is  constant 
between  samples,  the  polynomial  is  of  order  zero.  Figure  4.1  shows  the  general  A/D, 
D/A  processes  with  a  continuous  time  system.  An  exact  discrete  representation  of 
Figure  4.1  is  given  by  Equation  4.1  [4]. 

G(z)  =  (1-z-')z|5^| 

where  2  denotes  the  z-transform  of  G(s)/s  which  is  given  by 

where  denotes  the  inverse  Laplace  transform.  Sampled  data  systems  of  this  type 
may  therefore  be  analyzed  using  z-transform  techniques. 


[4.1] 

[4.2] 
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4.2  Digital  Compensator  Design  by  Emulation 


Continuous  compensators  are  designed  using  classical  continuous  control  tech¬ 
niques.  It  is  assumed  that  continuous  measurements  are  available  for  feedback.  The 
pitch  and  roll  axes  are  assumed  to  be  uncoupled  and  the  yaw  axis  is  assumed  to  be 
unaffected  by  main  rotor  speed.  Digital  filters  or  discrete  compensators  are  designed 
to  emulate  the  continuous  compensator  designs  which  are  based  on  the  above  as¬ 
sumptions.  The  bilinear  transformation  (Equation  4.3)  maps  the  stable  region  of  the 
s-plane  (3l?e(s)  <  0)  into  the  stable  region  of  the  z-plane  (|z(  <  1). 

2  z-1 
Tz  + I 

A  considerable  amount  of  distortion  takes  place  even  though  the  stability  regions  are 
coincident.  This  causes  the  frequency  response  characteristics  of  the  digital  filter  to 
differ  from  those  of  the  continuous  compensator.  To  make  the  frequency  responses 
of  the  digital  and  continuous  filters  equal  at  a  particular  frequency  a  prewarping 
factor  is  used  [4].  The  prewarped  discrete  equivalent  compensator  may  be  realized  by 
implementing  Equation  4.4. 


<on(wir/2)  »+l 


4.3  Digital  Filter  Programming 


Consider  a  digital  filter  of  the  form  given  in  Equation  4.5. 


K(z)  = 


0(z)  _  bp  +  biz  *  4- ...  -f  bjnz'”' 
I{z)  1 -f- ai2~*  4- 02^'^  +  . . .  a„z“" 


Since  2  '  represents  a  unit  delay,  a  difference  equation  for  the  filter  output  may  be 
written. 


o{k)  =  -aio{k-l)  —  a2o{k  —  2)  —  ...—ano{k  —  n)-\-boi{k)  +  bii{k  —  l)  +  ...  +  bmi{k  —  Tn) 

14.6) 

It  can  be  seen  that  n  -f-  m  delays  or  storage  registers  must  be  used  to  realize  this 
difference  equation.  Implementation  of  Equation  4.6  is  called  the  direct  programming 
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approach  [15]  .  The  number  of  storage  registers  may  be  reduced  from  n  +  m  to  n 
where  n  >  m  by  using  the  st2indard  programming  method  [15].  The  pulse  transfer 
function  0{z)H{z)  may  be  written 


0(z)  _  0{z)  Hjz) 
I(z)  H{z)I[z) 

where 

r  =  bo  +  b^z  *  +  ...  +  2  "* 

H{z) 

and 

H{z)  _  _ 1 _ 

I{z)  1  +  a,2"*  +  022-2  +  . .  .anZ-" 
The  following  difference  equations  may  now  be  written 


(4.7] 


o(A:)  =  boh{k)  +  bih{k  —  1)  +  b2h{k  —  2)  +  . . .  +  6„/i(fc  —  m)  [4.8] 

/i(k)  =  i(k)  —  aih{k  —  1)  —  a2h{k  —  2)  —  ...  —  amh{k  —  m)  [4.9] 

Note  that  only  n  storage  registers  are  required  to  obtain  ©(it).  This  saves  computer 
memory  and  processing  time.  The  difference  equations  are  implemented  as  dot  prod¬ 
ucts. 


/i(A:)  =  [-Oo  -oi  ...  -  a„]  •  [i(A:) /i(A:  -  1)  . . .  h{k  -  n)]  [4.10] 

o(fc)  =  [bobi  ...  6„]  •  [h{k)  h{k  -  1)  ...  h{k  -  m)]  [4.11] 


4.4  Open  Loop  System 

The  joystick  is  used  to  control  main  rotor  speed  and  rotor  tilt  angle  for  pitch  and 
yaw  control  respectively.  No  stability  augmentation  systems  are  implemented  in  this 
case.  Figure  4.2  shows  a  block  diagram  of  the  open  loop  systems.  Joystick  commands 
are  sampled  and  sent  to  the  actuator  every  0.05  sec.  For  the  purpose  of  analysis,  the 
pitch  0  and  roll  $  axes  are  assumed  to  be  uncoupled  since  the  primary  pitch  axis 
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Figure  4.2  Open  loop  command  system  block  diagram. 

controller  is  main  rotor  speed  0.  Figure  4.3  shows  the  Bode  plot  of  the  open  loop 
pitch  axis  system  of  Figure  4.2.  One  can  see  that  it  has  a  flat  frequency  response  out 
to  about  1.65  rad/sec  and  very  low  gain  (-35  dB).  The  Bode  plot  of  the  open  loop 
yaw  system  is  shown  in  Figure  4.4.  One  can  see  that  the  system  crossover  frequency 
is  at  2.7  rad/sec  and  the  phase  margin  is  15".  The  low  frequency  -20  dB/dec  slope  is 
indicative  of  the  presence  of  a  pure  integral  or  1/s  term  in  the  system  dynamics. 

Actuators  are  saturated  digitally  before  they  hit  their  physical  limits  for  safety 
reasons.  The  main  rotor  is  limited  to  speeds  between  240  and  360  rad/sec  while  the 
travel  of  the  tilt  servo  is  limited  to  the  region  —0.3  <  $  <  0.3  rad.  The  rate  limit 
on  the  main  rotor  discussed  in  Chapter  3  is  not  considered  since  its  effects  are  not 
significant  in  the  frequency  range  of  interest  (see  Chapter  5).  • 

4.5  Attitude  Command  Systems 
4.5.1  Yaw  Axis 

The  following  specifications  are  set  for  the  yaw  attitude  command  system 

•  Closed  loop  stability. 

•  Gain  Margin  >  6  dB 

•  Phase  Margin  (pm  >  40". 
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open  Loop  Yaw  System 


CD  (rad/sec) 


Figure  4.4  Bode  plot  of  open  loop  yaw  system 
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Phase  (deg) 


Figure  4.5  Closed  loop  continuous  yaw  attitude  command  system. 


•  Closed  Loop  -3  dB  cutoff  frequency  >  3  rad/sec. 

•  Steady  state  error  to  a  step  <  5%  (A'p  >19)  where  Kp  is  the  position  error 
constant  [16]. 


Figure  4.5  shows  a  block  diagram  of  the  closed  loop  yaw  axis,  tilt  servo  dynamics  and 
cascade  compensator.  The  steady  state  requirement  is  met  without  compensation 
because  of  the  pure  integral  term  in  the  yaw  axis  dynamics. 

153.42  7.16  7.16 


Kp  =  lim 


*-0  52^.14.185  + 153  42  5(5 +  .7153)  0 


=  cx> 


[4.12] 


From  Figure  4.6,  one  can  see  that  the  uncompensated  system  has  the  following  prop¬ 
erties: 


•  Gain  crossover  frequency  u}gc  =  2.7  rad/sec 

•  Phase  margin  =  15° 


•  Gain  margin  =  5.33  dB 


The  target  crossover  frequency  is  selected  to  be  Wi  =  3  rad/sec.  The  phase  deficiency 
at  u;t  =  3  rad/sec  is  43°.  A  lead  compensator  is  designed  to  contribute  60°  of  phase 
lead  at  3  rad/sec.  The  form  of  the  lead  compensator  is  given  in  Equation  4.13  . 


K{s) 


K{ts  +  1) 
ars  -b  1 


(4.131 


where 


1  -Sin^niax 
1  +  sin  (pmax 
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and 


1 


where  4>max  is  the  maximum  phase  leawi  contribution  of  the  filter.  In  this  case  a  = 
0.0718  and  r  =  1.244.  The  increase  in  gain  due  to  the  inclusion  of  the  lead  filter  at 
3  rad/sec  is  11.43  dB;  therefore,  the  filter  gain  is  set  to  -11.43  dB  or  ii'  =  0.26795. 
Theoretically  this  configuration  meets  the  specifications;  however,  it  is  experimentally 
determined  that  this  compensator  did  not  satisfy  the  steady  state  error  requirements. 
This  is  because  the  compensator  is  designed  based  on  a  linear  model  which  includes 
only  viscous  friction.  In  practice,  static  and  coulomb  friction  exist  and  these  phenom¬ 
ena  have  significant  effect  upon  steady  state  error.  A  lag  compensator  is  cascaded 
with  the  lead  compensator  designed  above  to  improve  steady  state  tracking.  The 
entire  compensator  is  given  by  Equation  4.14. 


K{s)  =  0.26795 


s  -f  .001  0.0893s  -I- 1 


|4.14| 


Figures  4.6  and  4.7  show  the  continuous  open  and  closed  loop  Bode  plots  for  the 
yaw  attitude  command  system.  One  can  see  that  the  compensated  system  meets  the 
specifications. 


•  Closed  loop  stability. 


•  Gain  Margin  =  10.43  dB 


•  Phase  Margin  =  66.93". 

•  Closed  Loop  -3  dB  cutoff  frequency  u-sdB  =  6.154  rad/sec. 

•  Zero  steady  state  error  to  a  step  input  {Kp  =  oo) 

Figure  4.8  shows  a  block  diagram  of  the  yaw  attitude  sampled  data  system.  The 
digital  compensator  K{z)  is  designed  to  emulate  the  continuous  compensator  designed 
above.  The  bilinear  transformation  with  prewarping  is  used  to  design  the  discrete 
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Open  Loop  ContinuoLS  Yaw  Attitixle  Command  System 


Figure  4.6  Bode  plot  of  open  loop  continuous  uncompensated  and  compensated 
yaw  attitude  command  system  ^'(ju>)/e»(ju;). 
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Phase  (deg) 


Closed  Loop  Continuous  Yaw  Attitude  Command  Syst 


Figure  4.7  Bode  plot  of  closed  loop  continuous  compensated  yaw  attit 

system  ^{ju>)/9com{jf*})- 


c 


8' 


fi 
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equivalent.  The  sampling  period  is  T=0.05  sec  (u;,  =  125.6  rad/sec)  which  exceeds 
the  sampling  theorem  requirement  since  the  closed  loop  bandwidth  is  u>-3jb  =  6.1504 
rad/sec.  The  warping  frequency  is  selected  to  preserve  the  frequency  response  of  the 
continuous  filter  at  cj  =  3  rad/sec.  The  discrete  equivalent  filter  is  given  below 


K{z)  = 


3.0104Z*  -  5.82762  -I-  2.8201 
2*  -  1.56192  +  0.5620 


(4.15) 


Figures  4.9  and  4.10  show  the  open  and  closed  loop  Bode  plots  of  the  discrete  system 
of  Figure  4.8.  The  properties  of  the  discrete  yaw  attitude  command  system  are: 

•  Closed  loop  stability. 

•  Gain  Margin  =  9.074  dB 

•  Phase  Margin  =  63°. 

•  Open  loop  gain  crossover  frequency  ujgc  =  2.7037  rad/sec 

•  Closed  Loop  -3  dB  cutoff  frequency  w-ajB  =  6.16  rad/sec. 

•  Zero  steady  state  error  to  a  step  input  (Kp  =  oo) 
which  meet  the  specifications. 

The  digital  filter  K(z)  is  implemented  using  the  standard  programming  technique 
discussed  in  Section  4.3.  Figure  4.11  shows  the  actual  and  simulated  step  responses 
of  the  yaw  attitude  command  system.  It  can  be  seen  that  the  desired  control  action 
is  achieved.  The  short  (<  0.5  sec)  time  delay  at  the  beginning  of  the  actual  step 
responses  is  likely  caused  by  static  friction.  The  actual  positive  step  response  does 
not  rise  or  attain  steady  state  as  quickly  as  the  actual  negative  step  response.  Recall 
that  the  viscous  yaw  axis  damping  coefficients  are  direction  dependent  and  the  linea*' 
model  used  an  average  of  the  two  coefficients.  This  explains  why  the  motion  of  the  yaw 
axis  is  more  damped  in  the  positive  direction  than  the  negative.  Both  positive  and 
negative  experimental  step  responses  indicate  that  there  is  more  damping  present  than 
predicted  by  the  linear  model.  This  is  most  likely  the  result  of  unmodeled  nonlinear 
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(0  (rad/sec) 


Figure  4.9  Bode  plot  of  open  loop  discrete  compensated  yaw  attitude  command 

system 
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Phase  (deg) 


Oosed  Loop  Discrete  Yaw  Attitude  Command  System 


(0  (rad/sec) 


Figure  4.10  Bode  plot  of  closed  loop  discrete  compensated  yaw  attitude  command 

system 


91 


Phase  (deg) 


coulomb  friction.  Despite  the  nonlinearities,  one  can  see  that  this  system  does  track 
attitude  commands  reasonably  well.  The  refined  model  of  Chapter  3  empirically 
models  the  effects  of  the  static  and  coulomb  friction  by  replacing  the  free  1/s  in  the 
denominator  of  the  yaw  tr2msfer  function  with  a  lag  term  l/(s  +  .7153).  Figure  4.12 
show  the  actual  and  simulated  results  using  the  refined  yaw  axis  model.  One  can  see 
that  the  refined  yaw  axis  model  predicts  the  behavior  rather  well. 

4.5.2  Pitch  Axis 

The  following  specifications  are  set  for  the  pitch  attitude  command  systems; 

•  Closed  loop  stability. 

•  Gain  Margin  >  6  dB 

•  Phase  Margin  ipm  >  40". 

•  Closed  Loop  -3  dB  cutoff  frequency  u7_3jb  >  3  rad/sec. 


•  Steady  state  error  to  a  step  <  5%  (Kp  >  19)  where  Kp  is  the  position  error 
constant  [16]. 


Figure  4.13  shows  a  block  diagram  of  the  closed  loop  pitch  axis  dynamics,  main  rotor 
actuator  and  cascade  compensator.  First  K{s)  is  found  to  satisfy  the  steady  state 
error  requirement. 


K„  =  lim 


K  655.2 


52  +  17.67s  +  655.2  s^  +  1.064s  +  2.722 


=  0.0184  K 


[4161 


therefore  K  =  1034.5.  The  open  loop  system  is  unstable  for  the  value  of  K  chosen 
as  indicated  by  the  negative  phase  margin  =  —4.0288".  A  lag-lead  compensator 
is  designed  to  improve  low  frequency  tracking  and  increase  the  phase  margin  at  the 
target  crossover  frequency  u>t  =  A  rad/sec.  The  continuous  compensator  which  met 
the  specification  is  given  in  Equation  4.17. 


ACs)  =  1034.5 


s  -f-  2  s  -f  0.5 
s  -1-  0.01  s  -1- 15 
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[4.171 


0  2  4  6  8  10 
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Figure  4.11  Simulated  and  2M:tual  yaw  attitude  command  system  step  responses 


93 


Yaw  Angle  (rad)  Yaw  Angle  (rad) 


Yaw  Attitude  Convnand  System  Step  Responses* 


0  2  4  6  8  10 


Time(sec) 


Tinne(sec) 

Figure  4.12  Simulated*  and  actual  yaw  attitude  command  system  step  responses.(* 

Refined  yaw  axis  model.) 
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Figure  4.13  Sampled  data  pitch  attitude  command  system. 


Figures  4.14  and  4.15  show  the  open  and  closed  loop  Bode  plots  for  the  continu¬ 
ous  pitch  attitude  command  system.  The  continuous  compensated  system  has  the 
following  characteristics: 

•  Gain  margin  =  13.98  dB 

•  Phase  margin  =  52.15" 

•  Open  loop  gain  crossover  frequency  =  4.27  rad/sec 


•  Closed  loop  -3  dB  cutoff  frequency  =  7.127  rad/sec 

•  Steady  state  error  to  a  step  0.8%  or  Kp  =  124 

Orie  can  see  that  the  design  meets  the  specifications. 

Figure  4.16  shows  a  block  diagram  of  the  closed  loop  sampled  data  pitch  axis 
command  system.  The  digital  compensator  (filter)  K{z)  is  designed  to  emulate  the 
continuous  compensator  K{s),  The  bilinear  transformation  with  frequency  prewarp¬ 
ing  is  used  to  accomplish  this  task.  The  warping  frequency  is  selected  to  be  wj  =  4 
rad/sec  and  the  sampling  period  is  chosen  to  be  T  =  0.05  sec  {u>,  =  125.6  rad/sec). 
The  sampling  frequency  exceeds  the  closed  loop  -3  dB  cutoff  frequency  by  a  factor  of 
17.64.  The  discrete  equivalent  compensator  K(z)  is  given  by  Ekjuation  4.18. 

798.372^  -  1500.30Z  -h  703.83 
2*  -  1.452  -1-  0.45 

Figures  4.17  and  4.18  show  the  open  and  closed  loop  Bode  plots  of  the  discrete  system 
of  Figure  4.16.  The  properties  of  the  discrete  pitch  attitude  command  system  are: 


95 


Open  Loop  Continuous  Pitch  Attitude  Command  System 


a 

■o 


c 

CD 

O 


o  (rad/sec) 


~  Magnitude  (Coniperisaled) 
--  PTtsse  (Corrpensaled) 

—  Magnitude  (Uncompensated) 
--  Phase  (Uncompensated) 


Figure  4.14  Bode  plot  of  open  loop  continuous  compensated  and  uncompensated 
pitch  attitude  command  system  ©(ju;)/ee{ju;). 
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Closed  Loop  Continuous  Pitch  Attitude  Command  System 


CO  (rad/sec) 


Figure  4.15  Bode  plot  of  closed  loop  continuous  compensated  pitch  attitude 
command  system  0(jui)/0com(j‘*’). 


Figure  4.16  Sampled  data  pitch  attitude  command  system. 


•  Gain  margin  =  11.87  dB 

•  Phase  margin  ^pm  =  46.03" 

•  Open  loop  gain  crossover  frequency  =  4.26  rad/sec 

•  Closed  loop  -3  dB  cutoff  frequency  w-ms  =  7.78  rad/sec 

•  Steady  state  error  to  a  step  0.8%  or  Kp  =  124 

which  meet  the  specifications.  The  digital  compensator  is  implemented  using  the 
standard  programing  technique  described  in  Section  4.3.  Figure  4.19  shows  actual  and 
simulated  step  responses  of  the  pitch  attitude  command  system.  One  can  see  that  the 
desired  control  action  is  present.  The  actual  response  to  a  positive  step  ©cm  is  lower 
than  predicted  for  the  first  four  seconds  but  approaches  the  simulated  steady  state 
value.  The  actual  response  to  a  negative  step  ©cm  overshoots  the  simulated  response 
during  the  first  second  and  reaches  steady  state  faster  than  the  simulated  response. 
This  is  likely  caused  by  the  unmodeled  main  rotor  speed  rate  limit.  Recall  that  the 
linear  model  of  the  main  rotor  accurately  described  the  actual  frequency  response  out 
to  25  rad/sec.  The  step  inputs  consist  of  an  infinite  sum  of  sines  of  different  amplitudes 
and  frequencies;  therefore,  one  can  expect  the  unmodeled  high  frequency  nonlinearity 
to  have  some  effect  on  the  step  response.  The  forcing  functions  in  the  human  pilot 
tests  consist  of  a  sum  of  5  low  frequency  (<4  rad/sec)  sine  waves;  therefore,  the 
linear  model  adequately  describes  the  main  rotor  dynamics  in  the  frequency  range  of 
interest. 


4.6  Rate  Command  Systems 


In  order  to  design  rate  command  systems,  rate  measurements  or  estimates  must 
be  available  for  feedback.  In  this  case,  estimates  of  the  angular  rates  are  used  because 
only  attitude  measurements  axe  available.  The  simplest  estimate  of  pitch  and  yaw 
rate  given  sampled  attitude  measurements  is  given  by  Equation  4.19. 


e(*)  = 


[4.19) 
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open  LcxDp  Discrete  Pitch  Attitude  Command  System 


<»  (rad/sec) 

—  Magnitude  (ConfjeraMed) 
*■  (Compensated) 


Figure  4.17  Bode  plot  of  open  loop  discrete  compensated  pitch  attitude  command 

system 
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Closed  Loop  Discrete  Pitch  Attitude  Command  System 


CO  (rad/sec) 


Figure  4.18  Bode  plot  of  closed  loop  discrete  compensated  pitch  attitude  command 

system  0(e'"^)/0com(e^'*'^). 
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Phase  (deg) 


Figure  4.19  Simulated  and  actual  step  responses  for  pitch  attitude  command 

system. 
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Unfortunately  the  attitude  measurements  are  contaminated  with  sensor  noise  and 
Equation  4.19  is  very  sensitive  to  noise.  Sensor  noise  in  the  attitude  measurements 
can  make  rate  estimates  calculated  from  Equation  4.19  have  a  large  noise-to-signal 
ratio,  especially  for  small  T.  Smoothed  rate  estimates  are  needed  and  obtained  from 
an  optimal  state  estimator  or  Kalman  filter. 

4.6.1  Kalman  Filter 

A  digital  steady  state  Kalman  filter  is  designed  to  obtain  reliable  rate  estimates 
in  the  presence  of  sensor  and  process  noise.  The  process  noise  crudely  accounts  for 
unknown  disturbances  and  model  imperfections. 

The  continuous  state  space  model  of  the  TROFM  is  discretized  for  T  =  0.05  sec. 

x(fc+ 1)  =  $x(fc)  +  ru(lf)  +  Gw{it)  [4.20] 
y(fc)  =  Cx(A:)  +  v(Ar)  [4.21] 
=  [  Hk)  Q{k)  '^{k)  0{k)  ^{k)  ^k)  eik)  ^{k)  ] 

u(fe)  = 

where  $  =  F  =  e^^drB,  w(fc)  is  the  process  noise  and  v(fc)  is  the  sensor 
noise  which  are  assumed  to  be  zero  mean,  Gaussian  “white”  random  sequences  with 
covariances  Q  and  R  respectively.  The  discrete  sensor  noise  covariance  is  assumed 
to  be  diagonal  in  form  and  the  elements  of  the  covariance  matrix  are  computed  from 
experimental  data  using  Equation  4.22. 

(‘<•221 

*=0 

where  v,{k)  is  a  sequence  of  N  actual  measurements  taken  with  the  TROFM  in  its 
nominal  configuration  (0  =  0  rad,^'  =  0  rad,$com  =  0  rad,/?  =  300  rad/sec).  The 
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experimentally  determined  discrete  sensor  noise  covariance  matrix  R  is  given  by 

1.7596  000 

0  5.060  X  lO--*  0  0 

0  0  1.825  X  10-®  0 

0  0  0  1.480  X  10-“ 

The  Kalman  filter  takes  the  form  of  the  prediction  estimator  of  Equation  4.24. 

x(fc  +  1 )  =  $x(fc)  +  ru(fc)  +  Lp[y (it)  -  Cx(ib)]  [4.24] 

where  Lp  is  the  Kalman  gain  matrix  and  the  over-bar  denotes  prediction  estimates. 
A  prediction  estimator  is  used  instead  of  a  current  estimator  to  compensate  for  the 
time  delay  between  the  measurement  and  D/A  conversion  caused  by  the  computer 
processing  time  required  to  generate  estimates  and  compute  control  outputs.  In  other 
words,  by  the  time  the  control  action  is  taken,  time  has  advanced  about  one  sample 
time.  If  a  current  estimator  is  used  the  estimates  are  about  one  sampling  period  old 
by  the  time  the  control  signal  is  sent  to  the  actuator. 

The  software  used  to  compute  the  feedback  matrix  solves  for  the  matrix  Le  asso¬ 
ciated  with  the  current  estimator  of  Ek]uation  4.25  (5). 

x{k)  =  x(jfc)  -I-  Lc(y()b)  -  Cx(fc))  [4.25] 

where  x(A:)  is  the  predicted  estimate  baised  on  a  model  prediction  from  the  previous 
estimate. 

x(fc)  =  $x(fc  -  1) -f  ru(A:  -  1)  [4-26] 

The  relationship  between  Lp  and  Lc  may  be  deduced  as  follows: 

•  substitute  Ek^uation  4.25  into  4.26  to  obtain 

x(fc)  =  $  [x(A:  -  1)  -I-  Lc  [y(fc  -  1)  -  Cx(jt  -  1)]]  -f  ru(jb  -  1)  [4.27] 

•  shifting  the  index  by  -1-1  one  obtains 

x(jt  -I- 1)  =  $x(it)  +  $Lc[y(Jfc)  -  Cx(ib)]  -t-  ru(ifc)  [4.28] 
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Comparing  the  structure  of  Equation  4.28  to  that  of  Equation  4.24,  one  may  deduce 
that  the  prediction  estimator  feedback  gain  matrix  Lp  is  related  to  Lc  by  Equa¬ 


tion  4.29. 


Lp  =  ♦Le 


14.29) 


The  design  software  determines  the  gain  matrix  Lc  such  that  the  discrete  current 
estimator  of  Eiquation  4.25  produces  an  optimal  estimate  of  the  state  vector  x  [5]. 
The  steady  state  Kalman  gain  matrix  is  computed  from  Ex}uation  4.30. 


L.  =  PC’'(CPC’'  +  R)- 


|4.30| 


where  P  is  the  unique  positive  definite  solution  to  the  algebraic  Riccati  equation 


P  =  4^  [P  -  PC^(CPC^  +  R)-‘CPj  +  GQG^ 


[4.31] 


The  matrix  P  is  the  a  priori  (before  the  measurement  update)  estimation  error  co¬ 


variance. 


A  pragmatic  design  rule  states  that  Lp  should  be  adjusted  such  that  the  residuals 
y  =s  {y(A;)  —  Cx(k))  lie  within  plus  or  minus  two  standard  deviations  of  the  variance 
of  the  residuals.  The  residual  covariance  is  given  by  Equation  4.32. 


Pi  =  CPC’^  +  R 


14.32) 


The  free  design  parameter  Q  allows  one  to  adjust  the  process  noise  and  thereby  adjust 
the  feedback  matrix  Lp.  Since  the  Kalman  filter  is  based  on  a  linear  model  of  the 
system  which  assumes  that  the  process  noise  has  the  same  power  at  all  frequencies 
(“white”)  and  it  is  known  that  the  actual  plant  has  several  low  frequency  nonlinear 
:haracteristics,  the  final  adjustment  of  Q  is  performed  by  analyzing  time  history 
lata.  The  PRBS  time  history  responses  of  Chapter  3  are  used  to  determine  how 
he  filter  performs  on-line.  The  process  noise  is  first  adjusted  to  meet  the  design 
ule  requirements.  Next  the  Kalman  filtered  rate  estimates  are  compared  to  the 
ilferentiated  attitude  measurements  and  Q  is  further  adjusted  to  produce  acceptably 
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smooth  rate  estimates.  The  process  noise  covsiriance  which  yielded  acceptable  results 
is  given  in  E^quation  4.33 


10 
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0 

[4.33] 

Q  = 

0 

0 

0 

0 

0 
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Figure  4.20  shows  the  Kalman  filtered  state  estimates  and  the  residuals  with  ±2<r(y) 
lines  where  applicable. 

The  TROFM  is  occasionally  required  to  operate  outside  of  the  range  where  the 
linear  model  is  valid  (i.e.  during  start-up  or  when  the  pilot  finds  it  necessary  in  order 
to  complete  a  task  (see  Chapter  5)).  The  pitch  axis  is  linearized  about  0  =  0  and 
the  linear  model  is  only  valid  for  small  deviations  about  0  =  0.  Recall  that  ground 
effect  is  modeled  as  a  linear  spring.  In  reality,  it  is  nonlinear  since  the  effective  spring 
“constant”  is  a  function  of  0.  For  large  positive  pitch  attitudes  Kge{Q)  =  0  and  for 
large  negative  pitch  attitudes  Kge{^)  ^  K''7e{0).  This  means  that  the  rotor  speed 
$  required  to  trim  the  TROFM  is  a  function  of  pitch  angle.  The  relationship  between 
$  and  0  is  experimentally  determined  and  is  modeled  by  the  third  order  polynomial 
fit  of  Equation  4.34. 

y9(0)  =  301 .2  -I-  22.83  0  -  1 13.8  0*  -I-  120.57  0^  [4.34] 

This  equation  modeled  the  trim  behavior  over  the  range  —0.6  <  0  <  0.6  rad/sec.  The 
estimated  trim  rotor  speed  is  subtracted  from  the  measured  absolute  ^  to  eliminate 
biased  0  and  0  estimates.  Without  accounting  for  the  biased  /?,  the  Kalman  filter 
overestimates  0  and  0  for  high  values  of  0  since  the  estimate  of  the  perturbation 
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Retidual  (nAhtt.)  Beiadat  (nd/acc)  Phi  Residuil  (lad)  Ai  (rad) 


0l2 


Time  (wc) 


Figure  4.20,  continue<i. 
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Figure  4.21  Block  diagram  of  hypothetical  continuous  pitch  rate  command  system 

A/?  is  significantly  higher  than  the  true  rotor  speed  perturbation. 

A/3  =  4  -  300  rad/sec  >  A/3  = /3  -  ^(0  »  0)  [4.35] 

Similarly  at  low  values  of  0  the  Kalman  filter  underestimates  0  and  0  because  the 
perturbation  A/3  is  significantly  lower  than  the  true  rotor  speed  perturbation. 

A/3  =  /?-  300  rad/sec  <  =  0  -}{Q  [4.36] 

Equation  4.34  is  used  to  preprocess  absolute  /3  measurements  to  yield  accurate  per¬ 
turbation  A0  estimates  which  are  required  for  proper  Kalman  filter  operation. 

4.6.2  Pitch  Rate  Command  System 

The  following  specifications  are  set  for  the  pitch  rate  command  system 

•  Closed  Loop  Stability 

•  Gain  margin  >  6  dB 

•  Phase  margin  >  40° 

•  Closed  loop  -3  dB  cutoff  frequency  >  6  rad/sec 

•  Steady  state  error  to  a  step  <  25 

Figure  4.21  shows  a  block  diagram  of  tne  hypothetical  continuous  pitch  rate  command 
system.  Figure  4.22  shows  a  Bode  plot  of  the  uncompensated  open  loop  system.  One 
can  see  that  the  uncompensated  system  has  the  characteristics  of  a  differentiator  at 
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low  frequencies  which  results  in  very  poor  steady  state  tracking  behavior.  In  order  to 
meet  the  steady  state  requirement,  an  integral  compensator  of  the  form  K{3)  =  K/s 
is  designed. 


655.2 


0.05  s 


K 


''V  =  +  17"6"77+655-2.-+ I.5S7+-2.7227  = 


The  minimum  gain  K  required  to  meet  the  steady  state  requirement  is  K  =163.32. 
The  gain  is  increased  to  K  =170  to  insure  that  the  requirement  is  met.  With  this 
compensation,  the  system  has  an  open  loop  crossover  frequency  of  u}gc  =3.2663  with  a 
gain  margin  of  12.66  dB  and  0"  of  phase  margin  at  6.366  rad/sec.  The  phase  deficiency 
is  43°  at  7.25  rad/sec;  therefore,  a  leawl  compensator  is  designed  to  contribute  65°  of 
phase  lead  at  this  target  crossover  frequency.  The  total  compensation  is  given  by 

Integral  Lead 


170  0.306s  H-  1 
s  0.0306s  +  1 


14.38) 


Figures  4.22  and  4.23  show  the  open  and  closed  loop  Bode  plots  of  the  compensated 
system.  The  compensated  system  has  the  following  characteristics: 


•  Gain  margin  =  9.38  dB 


•  Phase  margin  v?m  =  66° 

•  Open  loop  gain  crossover  frequency  =  5.98  rad/sec 

•  Closed  loop  -3  dB  cutoff  frequency  =  H-ll  rad/sec 


•  Steady  state  error  to  a  step  =  24.5 


which  meet  the  specifications. 

A  digital  filter  is  designed  to  emulate  the  continuous  compensator  of  Eiquation  4.38. 
It  is  assumed  that  continuous  pitch  rate  me<isurements  are  available  for  sampling.  Fig¬ 
ure  4.24  shows  a  block  diagram  of  the  hypothetical  sampled  data  pitch  rate  command 
system.  The  digital  compensator  is  computed  using  the  bilinear  transformation  with 


no 


(0  (rad/sec) 

—  Mignitude  (CcmperaatBd) 
-•  Phase  (Conpensaled) 

—  MmnAude  (UnoamperaaM) 
--  Phase (Unoonpensaled) 


Figure  4.22  Bode  plot  of  hypothetical  open  loop  continuous  compensated  and 
uncompensated  pitch  rate  command  system  0(ju;)/e@(jfa;). 


Ill 


Phase  (deg) 


(0  (rad/sec) 


4 

Figure  4.23  Bode  plot  of  hypothetical  closed  loop  continuous  compensated  pitch 
rate  command  system  0(ju;)/0coTn(j‘*’)- 
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Figure  4.24  Block  diagram  of  hypothetical  sampled  data  pitch  rate  command 

system 


prewarping  for  a  sampling  period  of  T  =  0.05  sec  and  warping  frequency  wj  =  6 
rad/sec. 

49.71z2  +  3.868z  -  45.839 
1.096Z  + 0.967 

Figures  4.25  and  4.26  show  the  open  and  closed  loop  Bode  plots  for  the  discrete 
system  of  Figure  4.24.  The  characteristics  of  the  discrete  compensated  system  are: 

•  Gain  margin  =  8.33  dB 

•  Phase  margin  ^pm  =  57.52® 

•  Open  loop  gzun  crossover  frequency  u)cg  =  5.96  rad/sec 

•  Closed  loop  -3  dB  cutoflf  frequency  u;_3jb  =  14.49  rad/sec 

•  Steady  state  error  to  a  step  =  24.45 

which  meet  the  specifications.  In  reality  continuous  pitch  rate  measurements  are 

•  E 

not  available  for  sampling  but  are  estimated  by  the  Kalman  filter.  Discussion  of  the 
effect  of  using  rate  estimates  instead  of  sampling  continuous  rate  measurements  is 
deferred  until  Section  4.6.4. 
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open  Loop  Discrete  Pitch  Rate  Command  System 


©  (rad/sec) 

Magnitude  (Conyxraaled) 
-•  Phase  (Con^pensMed) 


Figure  4.25  Bode  plot  of  hypothetical  open  loop  discrete  compensated  pitch  rate 

command  system  0(e'‘*'^)/eQ(e''"^). 


Phase  (deg) 


Figure  4.27  Block  diagram  of  hypothetical  continuous  yaw  rate  command  system 

4.6.3  Yaw  Rate  Command  System 

The  following  specifications  are  set  for  the  yaw  rate  command  system: 

•  Closed  Loop  Stability 

•  Gain  margin  >  6  dB 

•  Phase  margin  >  40° 

•  Closed  loop  -3  dB  cutoff  frequency  >  6  rad/sec 


•  Zero  steady  state  error  to  a  step;  Kp  =  oo. 

Figure  4.27  shows  a  block  diagram  of  the  hypothetical  continuous  yaw  rate  command 
system.  A  Bode  plot  of  the  uncompensated  system  is  shown  in  Figure  4.28.  In  order 
to  meet  the  steady  state  error  requirement  an  integral  controller  is  implemented 
K{s)  =  1/s. 


7.16 


153.42 


1 


_ _ ^  ^ 

s  +  0.7153  s*  + 14.18s  +  153.42  s 


[4.40] 


This  compensated  system  has  an  open  loop  gain  crossover  frequency  of  uigc  =  2.509 
rad/sec  with  a  1.312  dB  gain  margin  and  2.35°  of  phase  margin.  A  lead  compensator 
is  designed  to  increase  the  open  loop  crossover  frequency  to  Wjc  =  4  rad/sec  with  a 
phase  margin  of  47.33°.  The  total  compensation  is  given  in  Equation  4.41  . 


,  1  0.5074s +  0.9 

K{s)  =  - 


[4.41] 


s  0.0277s  +  1 

Figures  4.28  and  4.29  show  the  Bode  plots  of  the  open  and  closed  loop  compensated 
system.  The  system  properties  are  as  follows: 
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•  Gain  margin  =  8.80  dB 


•  Phase  margin  ipm.  =  47.34® 

•  Open  loop  gain  crossover  frequency  u>c,  =  4.01  rad/sec 

•  Closed  loop  -3  dB  cutoff  frequency  ui-uB  =  8.97  rad/sec 


•  Steady  state  error  to  a  step  =  0 

A  digital  filter  is  designed  to  emulate  the  continuous  compensator.  Emulation  is 
accomplished  using  the  bilinear  transformation  with  prewarping  for  a  sampling  period 
of  T  =  0.05  sec  and  a  warping  frequency  of  4  rad/sec.  The  discrete  equivalent 
compensator  is  given  by 

0.25182*  +  0.02152  -  0.2304 


K{z)^ 


14.42] 


z‘‘  -  1.0497r  + 0.04978 
Figures  4.31  and  4.32  show  the  open  and  closed  loop  Bode  plots  of  the  hypothetical 
discrete  system  of  Figure  4.30.  The  properties  of  the  discrete  yaw  rate  command 
system  are  as  follows: 

•  Gain  margin  =  7.235  dB 


•  Phase  margin  (pm  =  41.62® 

•  Open  loop  gain  crossover  frequency  u)cg  =  4  rad/sec 

•  Closed  loop  -3  dB  cutoff  frequency  u-sdB  —  8.97  rad/sec 

•  Steady  state  error  to  a  step  =  0 
which  meet  the  specifications. 

4.6.4  Effect  of  Kalman  Filter  in  Feedback  Loop 

Figure  4.33  shows  a  block  diagram  of  the  actual  rate  command  system  implemen¬ 
tation.  Note  the  presence  of  the  Kalman  filter  in  the  feedback  loop.  Obviously  this 
is  not  the  same  as  having  rate  measurements  available  for  feedback.  Figures  4.34 
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Open  Loop  Continuous  Yaw  Rate  Command  System 


--  Phase (Conpensaled) 

—  Magnitude  (Umonpersaied) 

-  -  F^dse  (Uncompensated) 


Figure  4.28  Bode  plot  of  hyptothetical  open  loop  continuous  compensated  and 
uncompensated  yaw  rate  command  system  ^'(ju;)/e*(ju;). 
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Phase  (deg) 


Closed  Loop  Continuous  Yaw  Rate  Command  System 


Figure  4.29  Bode  plot  of  hypothetical  closed  loop  continuous  compensated  yaw  rate 

command  system 

c 
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Phase  (deg) 


OR  EQUIVALENTLY: 
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Figure  4.32 


(0  (rad/sec) 


Bode  plot  of  hypothetical 


command  system  4’(e-'"^)/^com(e^‘^^)- 
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Figure  4.33  Block  diagram  of  actual  rate  command  system  implementation. 

and  4.35  compare  the  Bode  plots  of  the  single-input  single-output  (SISO)  closed  loop 
rate  command  systems  developed  above  to  the  rmlti-input  multi  output  (MIMO) 
model  using  Kalman  filtered  rate  estimates  for  feedbcick.  One  can  see  that  the  low 
frequency  behavior  of  the  yaw  rate  command  system  is  unchanged  but  the  resonant 
peak  is  approximately  2.3  dB  higher  and  the  -3  dB  cutoff  frequency  is  reduced  from 
8.97  rad/sec  to  about  8  rad/sec.  The  low  frequency  behavior  of  the  pitch  rate  com¬ 
mand  system  is  unchanged  but  the  resonant  peak  is  about  1  dB  higher  and  the 
magnitude  curve  drops  below  -3  dB  at  6.79  rad/sec  to  a  minimum  of  -4.2  dB  at  9.10 
rad/sec.  The  magnitude  curve  crosses  the  -3dB  line  again  at  13.64  rad/sec  and  finally 
drops  below  -3dB  at  14.05  rad/sec.  The  high  frequeii''y  behavior  of  the  system  is  vir¬ 
tually  unchanged.  Figures  4.36  shows  the  step  responses  of  the  p'*ch  and  yaw  rate 
command  systems.  The  SISO  simulated  rate  response,  online  Kalman  filter  rate  es¬ 
timate  and  differentiated  attitude  measurements  are  shown  for  comparison.  One  can 
see  that  the  Kalman  filter  rate  estimates  are  a  compromise  between  the  differentiated 
attitude  measurement  and  the  simulated  response. 
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Gain  (dB) 


Figure  4.36  Step  responses  of  SISO  simulated,  online  Kalman  filtered  rate  estimate, 
and  differentiated  attitude  measurements  for  a)  pitch  rate  b)  yaw  rate  command 

systems. 
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5.  MAN  IN  THE  LOOP  EXPERIMENTS 


A  study  of  the  effects  of  command  systems,  displays  and  tasks  on  human  pilot 
behavior  is  conducted.  Loren  P.  Dean,  a  general  aviation  pilot  uses  the  Cooper- 
Harper  rating  scale  to  evaluate  the  flying  qualities  of  the  TROFM  using  unaugmented, 
rate,  and  attitude  command  systems  for  two  types  of  tasks  and  two  different  display 
types.  Pilot  and  pilot-vehicle  describing  functions  are  obtained  from  time  histories 
of  dynamic  tracking  tasks  where  the  pilot  attempts  to  track  a  random  appearing 
input.  Compensatory  instrument  displays  and  a  visual  pursuit  scene  are  used  for  the 
dynamic  tracking  tasks.  The  basic  structure  of  a  continuous  time  single  loop  pilot- 
vehicle  system  where  the  pilot  attempts  to  track  a  random  appearing  input  (forcing 
function)  is  given  in  Figure  5.1.  Another  task  is  devised  which  requires  the  pilot 
to  move  a  stationary  object  from  one  platform  to  another  in  minimum  time.  The 
pilot  evaluates  the  suitability  of  the  three  command  systems  for  this  task  using  the 
Cooper-Harper  rating  scale.  Correlations  between  the  Cooper-Harper  ratings  and 
quantitative  descriptions  of  pilot  behavior  are  made.  The  experimental  results  are 


Figure  5.1  Basic  structure  of  a  continuous  pilot  in  the  loop  tracking  system. 


compared  to  the  classical  theory  of  man-machine  system  dynamics  which  is  based  on 
McRuer’s  Crossover  Model. 

5.1  Apparatus  and  Procedure 

The  command  systems,  special  equipment  and  experimental  procedures  used  in 
the  pilot  in  the  loop  tests  are  described  below. 

5.1.1  Control  Systems 

The  TROFM  is  configured  with  the  attitude  command,  rate  command,  and  open 
loop  systems  designed  in  Chapter  4.  Figure  5.2  shows  simplified  block  diagrams  of 
these  command  systems.  The  transfer  function  associated  with  the  rate  command 
systems  Q/Qcom  and  ^ /'^com  have  nearly  flat  0  dB  frequency  responses  out  to  ap¬ 
proximately  6.79  and  8  rad/sec  respectively.  Since  the  pilot  is  trying  to  control  0 
and  'P  by  manipulating  0  and  4’,  the  effective  controlled  element  is  Yc  KJs  at  low 
frequencies.  The  attitude  command  systems  0/0com  and  '‘i! /'^com  have  nearly  flat  0 
dB  frequency  responses  out  to  7.78  and  6.16  rad/sec  respectively.  This  makes  the 
effective  controlled  element  K  ss  1  in  the  low  frequency  region. 

5.1.2  Joystick 

A  position  hold  potentiometric  joystick  is  used  as  the  aircraft  flight  director.  Posi¬ 
tion  hold  describes  the  characteristics  of  a  joystick  which  maintains  a  fixed  orientation 
(using  friction)  in  the  absence  of  an  input.  Longitudinal  joystick  deflection  controls 
main  rotor  speed,  pitch  rate  0  and  pitch  attitude  0  for  the  open  loop,  rate  and 
attitude  command  systems  respectively.  The  lateral  portion  of  the  joystick  controls 
rotor  tilt  angle  $,  yaw  rate  and  yaw  attitude  for  the  open  loop,  rate  and  attitude 
command  systems  respectively.  The  1  KQ  joystick  potentiometers  are  configured  as 
voltage  dividers  powered  by  a  10  volt  power  supply.  The  wipers  are  connected  to  the 
A/D  board  on  the  control  computer. 
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Figure  5.2  Simplified  block  diagrams  of,  open  loop,  rate  and  attitude  command 

systems. 
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Joystick 


I _ I  Arm  Angle  with  Horizontal 

Flat  Frequency  Response  out  to  7.78  rad/sec 


_  Arm  Heading  Angle 

Flat  Frequency  Response  out  to  6.16  rad/sec 


Figure  5.2,  continued. 
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5.1.3  Pilot 


Loren  P.  Dean,  a  general  aviation  pilot,  evaluated  the  flying  qualities  of  the 
TROFM.  The  pilot  was  a  masters  student  in  dynamics  and  control  in  the  Purdue 
University  School  of  Aeronautics  and  Astronautics.  He  had  logged  230  hours  in  sin¬ 
gle  engine  land  general  aviation  aircraft  and  maintained  a  current  instrument  rating. 
He  had  no  prior  experience  evaluating  aircraft  handling  qualiti^’:.  The  pilot  was 
briefed  on  the  Cooper-Harper  rating  scale  and  the  scale  flowchart  (Figure  5.3)  was 
used  for  all  flying  qualities  evaluations.  Pilot  comments  were  recorded  on  audio  tape 
for  later  transcription. 

5.1.4  Experiment  #1 

The  pilot  is  assigned  the  task  of  maneuvering  the  TROFM  to  a  stationary  object 
and  moving  the  object  to  a  specified  position  on  an  elevated  platform  in  minimum 
time.  Figure  5.4  shows  a  diagram  of  the  apparatus  used  in  this  test.  A  tail  hook 
is  rigidly  attached  to  the  ballast  of  the  TROFM  for  the  purpose  of  picking  up  a 
small  light  weight  spool.  In  each  test  the  TROFM  started  from  a  position  halfway 
between  the  spool  and  the  target.  A  timer  is  started  by  flight  test  director  and  the 
pilot  maneuvers  the  TROF.M  to  the  spool  and  picks  it  up  using  the  tail  hook.  The 
pilot  proceeds  to  move  the  spool  to  the  elevated  platform  and  attempts  to  bring  the 
spool  to  rest  in  the  desired  performance  region  marked  on  the  elevated  platform. 
When  the  spool  is  at  rest  within  the  desired  or  adequate  performance  region  and  the 
tail  hook  removed,  the  timer  is  stopped  and  the  pilot  uses  the  Cooper-Harper  rating 
scale  to  evaluate  the  TROFM  command  system  being  tested.  The  performance  and 
suitability  of  each  command  system  for  the  task  is  measured  by  the  amount  of  time 
taken  to  complete  the  task  and  the  Cooper-Harper  handling  qualities  rating. 

5.1.5  Experiments  #2  and  #3 

The  pilot  is  assigned  the  task  of  tracking  random  appearing  inputs  and  keeping 
the  tracking  errors  within  desired  performance  bounds.  The  desired  performance 
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Figure  5.3  Cooper- liar  per  rating  chart 


Figure  5.4  Diagram  of  the  apparatus  used  in  Experiment 
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region  is  define*!  as  0e  <  0.07  amd  4fe  ^  0.10  rad.  Minor  excursions  out  of  the 
desired  performance  region  are  allowed  provided  they  are  not  excessive  or  prolonged 
and  minimal  pilot  compensation  is  required  to  get  back  within  the  desired  bounds. 

Low  (slow),  medium  and  high  (fast)  bandwidth  pseudo-random  forcing  functions 
are  used  in  the  dynamic  tracking  tests.  Each  forcing  function  consists  of  a  sum  of 
five  sine  waves. 


5 

sin(a;jfcr  +  <Pj)  (5.1) 

Each  forcing  function  contains  N  =  1024  discrete  points  which  are  sampled  at 
T  a  0.05  sec.  The  component  sinusoids  are  carefully  designed  to  have  an  integer 
number  of  cycles  in  1024  samples  and  no  frequency  component  is  an  integer  multiple  of 
another.  The  frequencies  are  roughly  equidistant  on  a  logarithmic  scale.  The  phase 
angles  Oj  are  randomly  selected  for  each  run  to  maintain  a  random  forcing  function 
appearance  while  preserving  RMS  magnitudes  and  spectra.  The  forcing  function 
signal  composition  for  sampling  frequency  of  F,  —  20  Hz  is  shown  in  Figure  5.5 
and  tabulated  in  Table  5.1.5.  Typical  forcing  function  time  histories  are  shown  in 
Figure  5.6.  In  multi-axis  tasks  there  are  no  common  frequencies  in  the  pitch  and  yaw 
disturbance  functions.  In  single  axis  tasks,  the  axis  not  in  use  is  physically  restrained 
and  the  joystick  potentiometer  controlling  that  axis  is  taken  off  line.  The  pilot  takes 
no  less  than  two  practice  runs  on  each  command  system  before  data  is  taken.  During 
this  practice  time  the  pilot  is  able  to  adjust  the  sensitivity  of  the  joystick  to  his  liking 
so  that  the  joystick  does  not  degrade  his  performance  or  affect  his  ability  to  assign 
an  appropriate  HQR. 

5.1.6  Instrument  Display 

In  Experiment  ^2  the  pilot  uses  a  compensatory  instrument  display  to  determine 
the  tracking  error.  A  compensatory  display  is  one  in  which  the  forcing  function 
minus  the  modified  control  response  is  indicated.  When  using  a  compensatory  display 
the  pilot  can  determine  the  effects  of  the  joystick  deflection  alone  only  under  zero 
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Table  5.1  Signal  composition  of  a)  Low  b)  Medium  s  c)  High  bandwidth  forcing 

functions. 


#  of  Cycles 

in  1024  samples 

Frequency 

(rad /sec) 

- - 

Amplitude 

(rad) 

2 

0.2453 

0.045 

3 

0.3680 

0.045 

5 

0.6133 

0.045 

11 

1.3492 

0.0225 

17 

2.0S52 

0.0225 

Frequency 

•Amplitude 

(rad /sec) 

(rad) 

0.4906 

0.045 

0.7359 

0.045 

1.1039 

0.045 

1.7172 

0.0225 

2.5758 

0.0225 

#  of  Cycles 

in  1024  samples 

Frequency 

(rad /sec) 

Amplitude 

(rad) 

5 

0.6138 

0.045 

8 

0.9813 

0.045 

12 

1.4719 

0.045 

18 

2.2078 

0.0225 

26 

3.1891 

0.0225 
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Figure  5.5  Forcing  function  signal  composition 
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input  conditions.  Figure  5.7  shows  the  instrument  panel  used  in  Experiment  ^2. 
This  panel  is  displayed  on  the  CRT  of  the  control  computer  and  is  generated  by 
National  Instruments  LabView  program.  The  compensatory  displays  are  made  using 
LabView’s  virtual  oscilloscope  elements  since  these  are  the  fastest  displays  available. 
The  scope  element  indicates  pitch  error  through  the  vertical  position  of  a  horizontal 
line.  This  error  indication  is  aligned  with  the  longitudinal  joystick  axis  which  controls 
the  pitch  axis.  Unfortunately  yaw  error  is  indicated  the  same  way  as  pitch  error  which 
is  counter-intuitive.  Yaw  error  would  be  better  represented  by  the  horizontal  position 
of  a  vertical  line.  This  would  be  aligned  with  the  lateral  joystick  axis  which  is  the 
primary  yaw  axis  controller.  Every  combination  of  available  elements  was  tried  in 
an  effort  to  obtain  an  intuitive  yaw  error  display  but  each  combination  caused  the 
sampling  rate  to  drop  to  unacceptably  low  levels;  therefore,  two  scope  elements  were 
the  only  choice  in  this  case. 

The  pilot  adjusted  the  joystick  sensitivity  to  his  liking  during  the  practice  runs 
using  the  sliders  shown  in  Figure  5.7.  After  the  practice  runs  the  pilot  is  given  20 
seconds  to  stabilize  the  aircraft  and  obtain  near  0  error.  The  “Taking  Data”  light 
on  the  panel  is  then  illuminated  and  the  forcing  function  is  sampled.  The  difference 
between  the  desired  attitude  and  the  actual  aircraft  attitude  is  displayed  on  the  virtual 
oscilloscopes.  Single  axis  tests  are  performed  first  to  allow  the  pilot  to  get  the  feel  of 
the  individual  axes  before  attempting  to  control  both  axes.  Discrete  time  histories 
of  the  forcing  function  i(kT),  vehicle  output  o(kT)  and  joystick  deflection  c(kT)  are 
recorded  on  diskette  for  analysis.  The  pilot  uses  the  Cooper- Harper  rating  scale  to 
evaluate  the  flying  qualities  of  each  configuration. 

5.1.7  Visual  Target 

In  Experiment  #3  the  pilot  is  assigned  the  task  of  tracking  a  moving  target. 
The  target  Is  mounted  on  a  custom  built  vertical  x-y  plotter  shown  in  Figure  5.8. 
A  high  intejibity  LED  is  enclosed  in  an  open  ended  tube  which  is  mounted  on  the 
ballast  of  the  TROEM.  The  pilot  is  asked  to  keep  the  light  beam  within  the  desired 
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LabView  Virtual  Oscilloscopea 


Longitudinal  Joystick  Sensitivity  Lateral  Joystick  Sensitivity 


Adjusment  Sliders 


Figure  5.7  Diagram  of  the  instrument  panel 
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Figure  5.8  Diagram  of  the  x-y  plotter  with  target 
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performance  bounds  defined  by  a  rectangular  card  with  cross-hairs  denoting  ideal 
performance.  The  motion  of  the  target  is  imparted  by  two  independent  computer 
controlled  stepper  motors  which  were  salvaged  from  two  Citizen  120D  dot  matrix 
printers.  Stepper  motors  are  chosen  as  motion  drivers  because  of  their  ability  to 
operate  as  open  loop  position  controllers  using  only  digital  outputs  from  the  control 
computer.  Both  analog  outputs  and  7  out  of  the  8  analog  input  channels  are  already 
in  use  by  the  TROFM;  thus,  open  loop  digital  control  is  the  only  choice  for  additional 
motion  driver  mechanisms.  The  external  circuitry  required  to  drive  the  steppers  is 
also  custom  built  (refer  to  Appendix  A  for  schematics).  The  target  is  able  to  move  in 
0.047”  increments  in  both  axes.  The  steppers  are  driven  by  the  sum  of  sines  forcing 
functions  previously  described,  however  in  in  the  multi-axis  tests  the  component 
amplitudes  are  decreased  by  10%,  to  compensate  for  a  decrease  in  sample  rate  that 
is  caused  by  excessive  computer  processing  time  required  to  drive  the  two  steppers. 

The  target  position  in  steps  is  correlated  with  the  attitude  of  the  TROFM  by 
starting  from  a  reference  point  and  moving  the  target  a  specific  number  of  steps  and 
manually  adjusting  the  TROFM  attitude  such  that  the  light  beam  is  centered  on  the 
target  cross-hairs.  The  pitch  and  yaw  axis  calibration  procedures  yield  the  following 
relations; 

rad 

Qcom  =  0.00138 -  ^Pitch  Steps  +  0.1825  rad  [5.2] 

step 

=  0.00129-^^^  aw  Steps  —  0.1571  rad 
step 

The  number  of  steps  to  move  in  each  sampling  period  T  is  calculated  using  Equa¬ 
tions  5.3  and  5.4. 


'AOjtep 

AT 


int  [  727. 2 {Qcomik  +  1)  -  -  0.1825)1 


[5.3] 


AT 


int  [  781.3 +  1)  -  'l'com(fc)  +  0.1571)] 


(5.4) 
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The  true  forcing  function  time  history  is  recorded  in  steps.  When  the  test  is  over 
the  true  forcing  function  time  history  in  radians  is  extracted  from  the  step  record  by 
Equations  5.5  and  5.6. 


0com(^  ) 


0.0138 


Yi  A0„,p(*:)  -  200steps 

U=o 


+  0.1825 


(5.5] 


'i'com(n) 


0.0129 


Y  +  200s<eps 


0.1571 


(5.6] 


Lk=o  J 

The  200  step  offsets  are  used  to  center  the  target  between  the  stops  that  are  used  as 
reference  positions. 


5.2  Discussion 


Analysis  of  Experiment  #1  is  straight  forward.  There  are  two  measures  of  perfor¬ 
mance;  handling  qualities  rating  and  time  taken  to  complete  the  task.  Experiments 
#2  and  #3  are  designed  to  contain  much  more  information. 

5.2.1  Experiments  #2  and  #3 

Pilot  describing  functions  and  open  loop  pilot-vehicle  describing  functions  are  de¬ 
termined  from  the  time  histories  recorded  in  Experiments  #2  and  #3.  Describing 
functions  are  extracted  from  the  discrete  time  histories  using  a  Fast  Fourier  trans¬ 
form  algorithm.  The  FFT  yields  r'^sults  equivalent  to  the  Discrete  Fourier  Transform 
(DFT)  but  the  number  of  data  points  (N)  in  a  set  must  be  a  power  of  two.  The 
data  sets  collected  in  Experiments  #2  and  ^3  contain  ;V  =  2’°  =  1024  points  to 
facilitate  fast  computation  of  frequency  responses.  Equation  .5.7  gives  the  definition 
of  the  DFT  of  a  general  discrete  time  sequence. 

N-l 

X(n)  =  =  DFT[x(kT)]  =  Y  (5.7] 

fc=0 

It  can  be  shown  that  the  DFT  is  defined  only  at  discrete  frequencies  a;„  =  2nFtn/N 
where  F,  is  the  sampling  frequency  in  Hz.  Consider  the  problem  of  determining 
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Figure  5.9  Block  diagram  of  a  sampled  data  system 


the  frequency  response  of  the  sampled  data  system  in  Figure  5.9  from  experimen¬ 
tal  time  history  data,  where  the  input  is  composed  of  a  finite  number  of  sinusoids 
and  the  sampling  frequency  is  much  greater  (in  this  case  40x  greater)  than  the  high¬ 
est  frequency  component  of  the  input.  Equation  5.8  mathematically  expresses  the 
relationship  between  the  transfer  function  G(z)  and  the  input-output  sequences. 


G{z) 


Ojz) 

I{2) 


[5.8] 


The  frequency  response  of  the  the  discrete  transfer  function  G(r)  can  be  determined 
from  Equation  5.9. 


)  =  Pfr|o(fer)i 
'  '  DFT[HkT)] 

where  Ui  is  a  component  frequency  of  the  sum  of  sines  input 


15.9] 


5.2.2  Open  Loop  Pilot-Vehicle  Describing  Functions 

Figure  5.10  shows  a  block  diagram  of  a  single  loop  sampled  data  pilot-vehicle 
tracking  system.  The  open  loop  pilot-vehicle  describing  functions  are  assumed  to  be 
of  the  form  of  the  extended  crossover  model  (Equation  5^10)  described  in  References 
[9]  and  [13]. 

L''„-ji(TW  +  o/w) 

YpYcij^)  = - ^ -  [5-10] 

ju: 
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Figure  5.10  Block  diagram  of  single  loop  sampled  data  pilot-vehicle  tracking  system. 

where  r  is  the  open  loop  time  delay  and  a  is  the  open  loop  phase  droop  parameter. 
The  experimental  open  loop  describing  function  {YpYc)^^^  is  determined  by  computing 
the  ratio  of  the  DFT  of  the  vehicle  output  o(kT)  and  the  DFT  of  the  error  e{kT) 
perceived  by  the  pilot  from  the  display. 


...7)  ^  DFT[o{kT)] 


[5.11] 


DFT[e(kT)]\ 

where  w,  are  the  n  component  frequencies  of  the  forcing  function  i{kT).  The  continu¬ 
ous  time  transfer  function  of  Equation  5.10  is  fitted  to  the  experimentally  determined 
open  loop  describing  function  function  of  Equation  5.11.  The  parameters  A',r,  and 
Q  which  best  fit  the  experimental  open  loop  describing  functions  in  a  leaist  squares 
sense  on  a  logarithmic  scale  are  calculated  by  the  procedure  outlined  below.  Note 
that 


[5.12] 


arg(VpVcOw))  =  +  o/-’) 


First  find  K  which  best  fits  the  experimental  magnitude  data  (measured  in  dB).  Since 
the  magnitude  of  A'/ju.’  system  has  a  20dB/dec  slope  on  a  logarithmic  frequency  scale, 
one  may  write 

=  |VpV;(e^-^)|„p  =  -201og,o(u.’.)  +  b  [5.13] 
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Gain  (db) 


where  b  is  the  point  where  the  magnitude  line  intersects  the  logjQ(lrad/sec)  =  0  axis 
(Figure  5.11)  . 

Compute  the  value  of  6  which  best  fits  the  experimental  magnitude  data  in  a  least 
squares  sense  by  solving  the  linear  algebra  problem  of  the  form  A6  =  Y  given  below. 

(5.141 

+  20/ogriou.v 

The  least  squares  estimate  of  b  is  given  by 

6  =  (A^A)-^A^Y  =  -  f^(y;  +  201og,ou.-.)  [5.15] 

The  least  squares  estimate  of  K  may  now  be  computed  from  Equation  5.16. 


r 

1 

Y\  +  20  log  10  U.'! 

; 

6  = 

1 

Yn  +  20/ogriou.v 

K  =  1050 


[5.16] 


Recall  that  a  K/ju  system  has  90°  or  1x12  (rad)  of  phase  lag  at  all  frequencies. 
To  obtain  le£ist  squares  estimate  of  r  and  a  the  residual  phcise  or  the  difference 
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between  the  experimentally  determined  phase  and  the  A'/jw  phase  must  be  computed. 


Or,,.  =  arg[y;j;(r'“'’’')l„,  -  (- j) 


15.171 


The  parameters  r  and  o  which  best  fit  the  residual  phase  data  are  computed  by 
solving  the  linear  algebra  problem 


‘  JWl 


r 

a 


X  =  (pres 


(5.181 


The  least  squares  estimate  of  the  parameter  vector  x  is  given  by 

X  =  [r  a]^  =  (A^ A)' ’Andres-  [5.19] 

Thus,  least  squares  estimates  of  the  extended  crossover  model  can  be  obtained  ana¬ 
lytically. 

5.2.3  One  Third  Law 


The  “one  third  law”  has  been  developed  from  the  crossover  model  to  describe 
the  performance  of  a  human  operator  in  a  closed  loop  compensatory  task  [9|.  This 
rule  <issumes  that  the  input  bandwidth  is  much  less  than  the  pilot-vehicle  crossover 
frequency,  relatively  small  remnant  or  noise  is  present  in  the  control  loop  and  the 
input  spectrum  is  rectangular  and  the  open  loop  pilot-vehicle  describing  function  is 
of  the  form  h  The  one  third  law  is  given  in  Equation  5.20. 

2 


tr?  3V.-J 


[5.20] 


where  is  the  variance  of  the  input,  el  is  the  variance  of  the  error,  u;*,  is  the 
effective  input  bandwidth  and  is  the  pilot- vehicle  crossover  frequency.  The  effective 
bandwidth  of  the  forcing  function  is  estimated  from  Equation  5.21  found  in  Reference 
[9]. 

[5.21] 
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where  is  the  power  spectral  density  of  the  input.  The  effective  bandwidths  of  the 
low,  medium  and  fast  bandwidth  forcing  functions  are  calculated  using  a  discretized 
Equation  5.21.  The  low,  medium  and  fast  effective  bandwidths  are  1.3684,1.8170  and 
2.0619  rad/sec  respectively.  The  variance  of  each  forcing  function  used  in  Experiment 
#2  and  the  single  axis  tests  of  Experiment  #3  is  =  0.595  rad.  The  variance  of  the 
multi-axis  forcing  functions  used  in  Experiment  ^3  is  cr  =  0.536  rad.  One  third  laws 
for  the  slow,  medium  and  fast  bandwidth  forcing  functions  are  given  below  for  input 
variances  of  0.595  rad. 


[5.22] 

[5.23] 


The  crossover  frequency  is  taken  to  be  equal  to  the  least  squares  estimate  of  K 
computed  by  the  method  outlined  in  the  previous  section. 

In  practice  the  one  third  law  must  be  modified  to  include  a  non-zero  intercept. 
The  existence  of  an  intercept  implies  that  zero  RMS  error  is  unattainable.  This  is 
because  the  pilot  is  not  a  linear  element  since  he  produces  a  remnant  or  pilot  induced 
noise  at  frequencies  other  than  those  which  compose  the  forcing  functions. 


5.2.4  Pilot  Describing  Functions 


Reference  [9]  gives  the  general  form  of  a  precision  pilot  model  for  single  loop 
systems,  and  the  form  for  simplification  of  two  of  the  dynamic  elements. 
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1 


where 

"""  TV  ~ 

and 

Ts  ^Tn,  +  ~  [5.25] 

Ufs 

The  major  mid-frequency  elements  are  the  gain  A'p,  time  delay  r  and  the  lead-lag 
or  lag-lead  series  equalization.  The  pilot  model  (Equation  5.26)  in  this  study  took 
advantage  of  the  very  low  and  high  frequency  approximations  shown  in  Equation  5.24. 


Yp  ^  Kp  { [5,26] 
\  YlagJ^'  +  1  / 

The  effective  time  delay  parameter  is  defined  by  =  (r  -f  TV).  Reference  [10] 
suggests  that  the  series  equalization  parameters  are  limited  to  the  following  ranges: 


0  <  Tuad  <  5sec  [5.27] 

0  <  Tug  <  20sec 

The  gain  Kp  is  adjusted  for  overall  system  stability  and  good  low  frequency  response. 
The  pilot  adjusts  Tuad  and  Tcag  such  that  the  system  can  be  stabilized  by  proper 
gain  selection  over  a  broad  frequency  range  and  VpYc  has  the  characteristics  of  a  K/s 
plant  over  a  considerable  frequency  range  centered  on  the  crossover  frequency. 

To  determine  appropriate  values  of  KpJ^iead,  Ttag^Te  and  Q  from  experimentally 
determined  frequency  response  data,  an  iterative  constrained  optimization  routine  is 
used  to  minimize  the  cost  function  defined  in  Equation  5.28. 


V'=  E 


W,=U/l 


0.02 


arg 


[DFT{e{kT))] 

( [gFr(c(A:7’))] 

I  [DFT{e{kT))] 


-gain{Yp„t) 


+ 


[5.28] 


t2 


~  Mg(yj,  eat) 


Where  the  gain  is  in  decibels  and  the  phase  is  in  degrees. 
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Minimizing  the  weighted  sum  of  squares  cost  function  has  been  shown  to  be  ef¬ 
fective  in  matching  frequency  responses  of  low  order  equivalent  systems  to  frequency 
responses  extracted  from  flight  time  histories  using  FFTs  [14].  A  constrained  opti¬ 
mization  routine  from  the  Matlab  Optimization  Toolbox  is  used  to  solve  the  opti¬ 
mization  problem. 


min  Tuad-,  Tug,  t*,  a)  Subject  to  0  <  Tuad  <  5se<{5.29] 

(Ap.  Ti,  ead>  T' Lag* 

0  <  Tiag  <  lOsec 

The  10  second  constraint  on  Tiag  is  chosen  because  the  lag  pole  IfT^ag  can  not 
be  accurately  determined  in  all  cases  since  the  lowest  frequency  component  sinusoid 
is  approximately  0.21  rad/sec  which  is  greater  than  the  minimum  lag  pole  location 
1/10  rad/sec;  thus,  estimates  of  T^ag  >  10  sec  can  only  be  interpreted  as  significant 
lag  generation  in  the  region  where  experimental  frequency  response  data  is  available. 

5.2.5  Product  Rule 

The  product  rule  is  an  empirical  formula  which  combines  single  axis  HQRs  to 
form  an  estimate  of  the  overall  multi-axis  HQR.  A  pilot  can  control  one  axis  better 
than  he  can  control  multiple  axes;  therefore,  one  can  expect  multi-axis  HQRs  to  be 
worse  than  any  single-axis  HQR.  The  claissical  product  rule  [9]  for  a  two  axis  system 
takes  the  form  of  a  general  hyperbola. 

Rm  —  o,R\  /?2  "I"  h/?2  "h  cR\  d  [5.30] 

The  regression  coefficients  (a,  b,  c,  d)  are  computed  from  a  least  squares  fit  of 
Equation  5.31  to  experimental  data. 

Rm  —  qRqR^  "h  bR<f  -1-  cRq  -|-  d  [5.31] 

where  Rq  and  are  single  axis  pitch  and  yaw  HQRs  respectively. 

The  regression  coefficients  are  computed  by  solving  the  linear  algebra  problem  of 
Equation  5.32 
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[5.32] 


•  - 

/2e,  ft*,  ft*,  fte,  1 

^e„  I 

A 

The  least  squares  solution  to  Equation  5.32  is  given  by 

x  =  [ahc^i]^  =  (A^A)-'A^Rm  [5.33] 

It  should  be  noted  that  hyperbolic  product  rules  have  some  serious  shortcomings. 
In  most  cases  if  a  rating  in  one  axis  is  extremely  good,  the  product  rule  will  predict 
the  two  axis  rating  to  be  better  than  the  single  axis  ratings.  These  counter-intuitive 
predictions  arise  from  the  product  rule’s  lack  of  recognition  of  the  HQR  scale.  The 
HQR  scale  is  limited  to  values  between  1  and  10  while  the  general  hyperbolic  equation 
is  not  limited.  One  must  therefore  use  judgement  when  interpreting  the  results  of 
any  product  rule  regression. 

5.3  Results 

5.3.1  Experiment  #1 

The  amount  of  time  taken  to  move  a  stationary  object  between  the  two  platforms 
and  the  handling  qualities  rating  assigned  to  a  configuration  by  the  pilot  are  the  two 
measures  of  performance  used  in  this  test.  Table  5.2  shows  the  results  of  Experiment 
#1- 

The  single  letters  (A,B,C)  denote  the  first  runs  performed  at  the  beginning  of  this 
study.  The  double  letter  runs  are  performed  after  all  of  the  experiments  have  been 
run  to  see  how  experience  affects  both  HQR  and  time  performance.  One  can  see  that 
experience  significantly  reduces  the  amount  of  time  taken  to  complete  the  task  in  all 
categories  but  has  little  effect  on  the  HQR.  The  HQR  margin  of  error  is  found  to  be 


a 

b 

c 

d 

X 


ftm, 

ftmn 
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Table  5.2  Results  of  Experiment  #1. 


Run 

Command  System 

HQR 

Time(sec) 

D 

Open  Loop 

50.67 

AA 

Open  Loop 

5.5 

36.15 

B 

Attitude 

36.70 

BB 

Attitude 

ra 

19.70 

C 

Rate 

mm 

37.00 

CC 

Rate 

4 

27.00 
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+/-  .79  (HQR)  for  all  experiments.  This  shows  that  the  pilot  HQR  is  a  reliable  and 
consistent  measure  of  performance.  Table  5.2  shows  that  the  pilot  achieved  the  best 
performance  using  the  attitude  command  system.  The  attitude  command  system  also 
receives  the  best  HQRs  .  The  pilot  assigns  better  HQRs  to  the  rate  than  the  op)en 
loop  command  system.  He  also  achieves  better  performance  using  rate  rather  than 
the  open  loop  command  system. 

5.3.2  Experiments  #2  and  #3 
5.3.2. 1  Comparison  of  Displays 

The  visual  scene  is  clearly  preferred  to  the  instrument  display  as  can  be  deduced 
from  Figure  5.12  and  Table  5.3. 

It  is  assumed  that  if  a  configuration  receives  a  HQR  of  1  in  the  instrument  task, 
then  the  visual  task  will  also  receive  a  HQR  of  1.  The  best  fit  lines  computed  from 
Ic'ast  squares  fits  are  therefore  constrained  to  pass  through  the  point  (1,1).  A  line  with 
a  slope  of  1  indicates  that  the  pilot's  HQR  is  unaffected  by  display  type.  Lines  with 
slopes  less  than  1  indicate  that  it  easier  for  the  pilot  to  perform  dynamic  tracking 
tasks  using  the  visual  scene.  The  slopes  of  the  best  fit  lines  for  single  axis  pitch  and 
yaw  cases  are  0.4793  and  0.6932  respectively.  A  plot  for  the  multi-axis  raise  is  not 
shown  because  Table  5.3  clearly  shows  that  the  visual  scene  is  preferred.  Many  of  the 
multi-axis  instrument  tasks  receive  HQRs  of  10  indicating  that  control  is  lost  during 
some  part  of  the  operation. 

It  is  not  surprising  that  the  pilot  prefers  the  visual  scene  to  the  instrument  display 
because  the  pilot  can  perceive  more  information  from  the  visual  scene.  When  using 
the  visual  scene,  the  pilot  can  see  and  hear  what  the  individual  components  of  the 
TROFM  are  doing  at  any  given  time.  This  generates  the  possibility  of  the  pilot  open¬ 
ing  more  feedback  channels  and  allows  the  pilot  to  act  as  a  multi-variable  feedback 
control  system  to  obtain  desired  performance.  The  pilot  is  also  able  to  determine 
whether  an  error  increase  is  caused  by  joystick  deflection  or  forcing  function  which 
is  not  possible  with  the  compensatory  instrument  display.  Pilot  comments  suggest 
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Table  5.3  Handling  qualities  ratings  for  instrument  display  and  visual  scene  with 

identical  tasks. 


.Axis/Command  System 

Instrument  HQR 

Visual  HQR 

Pitch/Attitude 

2 

1 

Pitch/ Attitude 

6.5 

2 

Pitch/At  litude 

7.5 

3.5 

Pitch/Rate 

6 

4 

Pitch/Rate 

7.5 

4.5 

Pitch/ Rate 

6.5 

5.5 

Pitch/Open  Loop 

3.5-^4 

3 

Pitch/Open  Loop 

4 

3 

Pitch/Open  Loop 

8 

3.5 

A’aw/ Attitude 

4 

3 

Yaw/Attitude 

3.5 

3^3.5 

Yaw/Attitude 

3.5 

3.5^4 

Yaw  /  Rate 

3.5 

3 

Yaw/ Rate 

5.5 

3.5 

Yaw/ Rate 

6.5 

3.5 

Yaw/Open  Loop 

4.5 

4 

Yaw /Open  Loop 

8 

5.5 

Yaw/Open  Loop 

9 

8 

Multi/Attitude 

8 

3 

Multi/Attitude 

8 

5.5 
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Table  5.3,  continued. 


Axis/Command  System 

Instrument  HQR 

X’isual  HQR 

Multi/.\ttitude 

9.5 

4.5 

Multi/Rate 

9 

6.5 

Multi/Rate 

10 

7 

Multi/Rate 

10 

9 

Multi/Open  Loop 

10 

9 

Multi/Open  Loop 

10 

8 

Multi/Open  Loop 

10 

7 
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a) 


10 
9 
8 
7 
6 
5 
4 
3 
2 
1 

1  23456789  10 

Instrument  HQR 


b) 


r  i  :  1 _ i _ I _ : _ 1 _ ^ _ I 

1  23456789  10 


Instnjmenl  HQR 


Figure  5.12  Visual  versus  instrument  HQRs  for  a)  pitch  axis 


b)  yaw  axis  tasks 
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Command  Type 

*  HQR  Ol  Enor  .79  d  « 


Figure  5.13  HQR  cis  a  funct  ion  of  command  system  and  forcing  function  bandwidth 
using  instrument  display  for  single  axis  pitch  task. 

that  the  multi-axis  instrument  tasks  would  have  been  much  easier  if  the  yaw  error 
indicator  was  aligned  with  the  lateral  joystick  axis. 

5. 3. 2. 2  Instrument  Display 

Figures  5.13,  5.14  ’.15  and  Table  B.l  show  the  handling  qualities  ratings  for  the 
instrument  test  as  a  function  of  command  system  and  forcing  function  bandwidth. 

From  Figure  5.13  one  can  see  that  for  the  single  axis  pitch  task,  HQRs  increase 
as  task  bandwidth  increases  for  all  command  systems  except  rate,  where  the  order 
goes  slow,  fast,  medium.  The  attitude  command  system  receives  the  best  rating  (2) 
for  the  slow  task  while  the  open  loop  receives  the  best  rating  for  the  medium  task  (4) 
and  the  rate  command  system  receives  the  best  rating  (6.5)  for  the  fast  task. 

From  Figure  5.14  one  can  see  that  the  yaw  attitude  command  system  receives  good 
HQRs  for  all  forcing  function  bandwidths.  The  open  loop  command  system  receives 
better  ratings  than  the  rate  command  system  for  all  bandwidth  forcing  functions 
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Instnunent  Dirolay 

SrigkAMY^ 


Op^  Loop  Rate  Attitude 

Command  Type 


Figure  5.14  HQR  as  a  function  of  command  system  and  forcing  function  bandwidth 
using  instrument  display  for  single  a.xis  yaw  task. 


Instrument  Display 

Mufli-Axs  ' 


0  I - ^ ^ - ) 

Open  Loop  Rate  Attitude 

Command  Type 


*K3W>i^yiclEffcr**^.78 


Figure  5.15  HQR  as  a  function  of  command  system  and  forcing  function  bandwidth 
using  instrument  display  for  multi-axis  tasks. 
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Visual  Display 


CommandType 


Figure  5.16  HQR  cis  a  function  of  command  system  and  forcing  function  bandwidth 
using  visual  scene  for  single  axis  pitch  task. 

and  is  rated  better  than  the  attitude  command  system  for  tracking  the  slow  forcing 
function. 

Figure  5.15  shows  that  the  open  loop  command  system  is  uncontrollable  in  the 
multi-axis  task.  The  rate  command  system  is  uncontrollable  in  all  cases  except  for  the 
combination  of  low  bandwidth  pitch  and  medium  bandwidth  yaw  forcing  functions. 
The  attitude  command  system  is  controllable  in  all  cases  but  receives  very  poor 
ratings.  The  pilot  stated  that  the  confusing  yaw  axis  instrument  display  made  the 
multi-axis  tasks  more  difficult  than  they  should  have  been. 

5. 3. 2. 3  Visual  Scene 

Figures  5.16,  5.17  and  5.18  show  HQR  as  a  function  of  command  system  and 
forcing  function  bandwidth  for  Experiment  #3. 

Figure  5.16  shows  that  for  the  single  axis  pitch  cases,  high  bandwidth  tasks  receive 
higher  (worse)  HQRs  than  low  bandwidth  tasks  in  all  cases.  The  attitude  command 
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Command  Type 


Figure  5.17  HQR  as  a  function  of  command  system  and  forcing  function  bandwidth 
using  visual  scene  for  single  axis  yaw  task. 


Visual  Display 

krWti-AxiS  ' 


Open  Loop  Rate  Attitude 

Command  Type  - 


Figure  5.18  HQR  as  a  function  of  command  system  and  forcing  function  bandwidth 

using  visual  scene  for  multi-axis  tasks. 
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system  receive  the  best  overall  HQRs  followed  by  the  open  loop  and  rate  command 
systems. 

Figure  5.17  indicates  that  the  rate  and  attitude  command  systems  essentially  tie 
(i.e.  within  the  +/'.79  HQR  margin  of  error)  for  the  best  HQRs  while  the  open  loop 
command  system  is  clearly  not  the  command  system  of  choice. 

Figure  5.18  illustrates  that  the  open  loop  command  system  receives  poor  ratings, 
but  the  ratings  improve  as  the  bandwidth  of  the  tasks  increase.  The  rate  and  attitude 
command  system  HQRs  improve  as  the  bandwidth  of  the  tasks  decrease.  The  attitude 
command  system  receives  the  best  overall  HQRs. 

5. 3. 2. 4  Product  Rule 

A  product  rule  regression  is  performed  on  the  HQR  data  (Table  5.4)  from  the 
visual  dynamic  tracking  test  (Experiment  #3).  An  estimate  of  the  multi-axis  HQR 
given  the  two  single  axis  HQRs  can  be  obtained  from  Equation  5.34. 

=  -0.6779 Re R^  +  5.8322R4-  +  3.5767/?e  -  17.71  (5.34) 

One  can  see  that  the  product  rule  estimates  predict  the  actual  multi-axis  HQR  to 
within  -I-/-  1.02  rating  for  7  out  of  the  9  cases.  It  should  be  noted  that  only  the  7 
“good”  points  in  Table  5.4  are  used  in  the  calculation  of  the  regression  coefficients. 

The  iso- rating  curves  in  Figure  5.19  indicate  that  the  multi-axis  rating  degrades 
£is  either  single  axis  rating  degrades.  Upon  close  examination  one  can  also  deduce 
that  a  small  degradation  in  the  single  axis  yaw  HQR  will  result  in  large  multi-axis 
HQR  degradation.  One  may  also  see  that  HQR  degradations  in  the  pitch  axis  have 
a  less  pronounced  effect  upon  the  multi-axis  HQR.  This  indicates  that  the  yaw  axis 
is  the  most  difficult  to  control  in  multi-axis  tasks  and  one  can  expect  the  the  pilot 
to  concentrate  more  on  the  yaw  axis  than  on  pitch  axis.  Note  that  the  product  rule 
computes  ratings  which  are  beyond  the  bounds  of  the  Cooper-Harper  scale  since  it 
does  not  recognize  these  bounds.  As  stated  previously  product  rules  do  not  work  if 
either  or  both  single  axis  ratings  are  extremely  good.  In  this  case  experimental  data 
is  available  for  a  limited  number  of  cases  so  the  product  rule  should  not  be  considered 
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Yaw  HQR 


Pitch  HQR 


Figure  5.19  Product  rule  for  visual  scene. 
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Table  5.4  Actual  single  and  multi-axis  HQRs  with  product  rule  estimates  of 
multi-axis  HQR  (*  Not  used  in  regression). 


Pitch  HQR 

Yaw  HQR 

Multi-axis  HQR 

Est.  .Multi  HQR 

Error 

2 

3 

3 

2.87 

0.13 

3.5 

3 

4.5 

5.19 

-0.69 

2 

3.5 

5.5 

5.11 

0.39 

4.5 

3 

7 

6.73 

0.27 

4 

3.5 

6.5 

7.51 

o 

* 

4.5 

3.5 

9 

8.12 

0.88 

3 

5.5 

9* 

13.9 

4.91 

3 

4 

8 

8.21 

-0.21 

3.5 

5.5 

7* 

13.83 

6.84 

accurate  far  outside  of  the  range  covered  by  the  experimental  data  (2  <  ^  <  4  and 
3  <  /?♦  <  4.5). 

5. 3. 2. 5  Assessment  of  adherence  to  crossover  law 

Recall  that  the  crossover  law  states  that  in  a  closed  loop  compensatory  task,  the 
pilot  will  compensate  so  that  the  combined  pilot- vehicle  open  loop  dynamics  is  similar 
in  appearance  to  a  Kjs  plant  with  a  pure  time  delay  in  the  region  of  crossover.  The 
Bode  plots  of  the  experimentally  determined  describing  functions,  as  well  as  best 
fit  extended  crossover  models  to  the  experimental  data  are  shown  in  Appendix  B. 
The  extended  crossover  model  fits  the  experimental  data  well  in  caises  where  the  pilot 
assigns  HQRs  <  7.  This  is  as  expected  because  HQRs  >  7  indicate  that  controllability 
is  in  question. 

The  one  third  law  condenses  the  information  shown  in  the  open  loop  Bode  plots. 
Figures  5.20,  5.21  and  5.22  show  the  one  third  laws  for  the  slow,  medium  and  fast 
runs  for  the  visual  scene  test.  Figures  5.23,  5.24  and  5.25  shows  the  one  third  laws 
for  the  slow,  medium  and  fast  runs  for  the  instrument  test.  One  can  see  that  the 
modified  one  third  laws  describe  the  trends  in  the  all  of  the  visual  data  well  and  most 
data  scatter  can  be  described  by  adjusting  the  intercept  of  the  one  third  law  line. 
The  modified  one  third  laws  describe  the  medium  and  high  bandwidth  instrument 
cases  fairly  well  with  the. exception  of  the  cases  affected  by  the  crossover  regression 
phenomenon  [13].  Crossover  regression  occurs  when  the  pilot  reduces  his  gain  to 
reduce  the  pilot-vehicle  crossover  frequency  to  maintain  a  reasonable  tracking  error. 
This  can  also  occur  when  the  pilot  focuses  his  attention  on  one  axis  in  a  multi-axis 
task  which  opens  the  loop  on  the  ignored  axis.  One  can  see  that  most  instances  of 
crossover  regression  occur  when  the  pilot  is  involved  in  a  medium  or  fast  multi-axis 
tracking  task.  These  plots  also  indicate  how  effective  the  pilot  is  in  following  the 
forcing  function.  The  ratio  egjai  can  be  thought  of  as  normalized  performance.  A 
ratio  of  one  indicates  that  no  reduction  of  error  is  taking  place.  Ratios  less  than  one 
indicate  that  the  pilot  is  successfully  tracking  the  forcing  function  and  is  effectively 
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One  Thifd  Law  For  Slow  Visual  Runs 


Figure  5,20  One  third  law  for  the  visual  scene  with  low  bandwidth  forcing  functions. 


reducing  tracking  errors.  Figures  5.26  and  5.27  show  the  open  loop  describing 
functions  and  best  fit  crossover  models  for  the  pilot-vehicle  using  the  instrument 
display  and  pitch  and  yaw  attitude  command  systems  respectively.  Recall  that  the 
pitch  and  yaw  attitude  command  systems  had  nearly  flat  0  dB  frequency  responses 
out  to  7.78  and  6.61  rad/sec  respectively.  The  fact  that  the  crossover  model  fits 
the  experimental  data  well  indicates  that  the  pilot  is  generating  lag  compensation 
in  the  region  of  crossover.  Upon  examination  of  the  pilot  describing  functions  for 
these  cases  (Figures  5.28  and  5.29)  one  can  see  that  the  pilot  is  indeed  generating  lag 
compensation  which  is  indicated  by  the  negative  slopes  of  the  magnitude  curves.  The 
best  fit  pilot  describing  functions  associated  with  Figures  5.28  and  5.29  are  found  to 
be 


Pitch: 


Yaw: 


V  —  n  _-j(0.0385..^-t-0.3696/w) 

"  ~  ^  2.045s  +  1 


v  —  7  co9oiLlIIrf_jlJi„-j(0.3507u/-(-0.1769/u/) 
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[5.35] 

[5.36] 
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One  Third  Law  For  Medium  Visual  Runs 


Figure  o.i 


One  third  law  for  visual  scene  with  medium  bandwidth  forcing 
functions. 


One  Third  Law  For  Fast  Visual  Runs 


1/Kbf  S  (1/(i>e) 


Figure  5.22  One  third  laws  for  the  visual  scene  with  fast  bandwidth  forcing 

functions. 


One  Third  Law  For  Slow  Instnimeot  Runs 


Figure  5.23  One  third  laws  for  the  instrument  display  with  low  bandwidth  forcing 

functions. 


One  Third  Law  For  Medium  Instrument  Runs 
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Figure  5.24  One  third  laws  for  the  instrument  display  with  medium  bandwidth 

forcing  functions. 
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One  Third  Law  P<x  Fast  Instrument  Runs 


Figure  5.25  One  third  laws  for  the  instrumeiii  display  with  fast  bandwidth  forcing 

functions. 
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Run  1  Ooen  Loon  Describing  Function 


Cl)  (rad/sec) 


Figure  5.26  Open  loop  pilot-vehicle  describing  function  for  a  low  bandwidth  single 
axis  pitch  forcing  function  with  instrument  display  and  attitude  command  system. 


[b]  Figures  5.30  and  5.31  show  open  loop  pilot- vehicle  describing  functions  and  best 
fit  extended  crossover  models  for  pitch  and  yaw  rate  command  systems  used  with 
instrument  displays.  One  can  see  that  the  crossover  model  describes  the  trends  quite 
well. 

Recall  that  the  pitch  and  yaw  rate  command  systems  have  nearly  flat  0  dB  fre¬ 
quency  responses  out  to  approximately  6.79  rad/sec  and  8  rad/sec.  Since  the  con¬ 
trolled  elements  K  are  simple  integrators  in  the  forcing  function  frequency  range, 
one  can  expect  that  the  pilot  will  act  eis  a  pure  gain  element  to  make  the  open 
loop  pilot- vehicle  describing  function  appear  eis  a  k/s  plant.  Figures  5.32  and  5.33 
show  pilot  describing  functions  associated  with  the  open  loop  describing  functions 
of  Figures  5.30  and  5.31.  One  can  see  that  the  p:  ot  is  primarily  acting  as  a  pure 
gain  compensator  with  time  delay  as  indicated  by  the  relatively  flat  magnitude  and 
slightly  negative  phase  curves.  Estimates  of  the  pitch  and  yaw  describing  functions 
are  given  in  Equations  5.37  and  5.38. 


V'  -  7  316/'^'^^“^'^  ^  -j(0.1554u;-0.0339/^) 

tp-  1.3164  e 


[5.37] 
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Run  4  Open  Loop  Describing  Function 


(u  (rad/sec) 


Figure  5.27  Open  loop  pilot-vehicle  describing  function  for  a  low  bandwidth  single 
axis  yaw  forcing  function  with  instrument  display  and  attitude  command  system. 
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Figure  5.28  Pilot  describing  function  for  a  low  bandwidth  single  axis  pitch  forcing 
function  with  instrument  display  and  attitude  command  system. 
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Run  4  Pilot  Describing  Function 
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Figure  5.29  Pilot  describing  function  for  a  low  bandwidth  single  axis  yaw  forcing 
function  with  instrument  display  and  attitude  command  system. 


Vp  =  1.4513  -;(o.49i5w+o.i:69/.)  3g| 

The  open  loop  yaw  command  system  shown  in  Figure  4.2  has  the  frequency  re¬ 
sponse  characteristics  shown  in  Figure  4.4.  One  can  see  that  it  is  primarily  a  rate 
command  system  with  a  crossover  frequency  at  2.43  rad/sec  and  has  very  little  phase 
margin.  This  suggests  that  the  pilot  will  have  to  generate  lead  compensation  to  in¬ 
crease  the  phase  margin  when  the  forcing  function  bandwidth  approaches  the  yaw 
axis  crossover  frequency.  Figure  5.34  shows  the  open  loop  pilot-vehicle  describing 
functions  for  a  high  bandwidth  =  2.06  rad/sec)  forcing  function.  One  can  see 
that  the  crossover  model  fits  the  magnitude  and  phase  data  well.  Figure  5.35  shows 
the  pilot  describing  function  for  this  case.  Note  the  positive  slope  of  the  magnitude 
curve  around  2  rad/sec  indicates  the  pilot  is  generating  lead  compensation.  The  phase 
curve  is  slightly  above  0°  between  1  and  2  rad/sec  and  then  sinks  below  0°  above  2 
rad/sec  which  indicates  that  the  pilots  phase  lead  is  counteracted  by  his  reaction  time 
delay.  The  best  fit  pilot  describing  function  for  this  case  is  given  in  Equation  5.39. 


Y  _  n  .5qi4iM?^?ijti]p-r(0  1275w+1.2351/u,) 

0.1157s  4-1 


[5.39] 
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Figure  5.30  Open  loop  pilot-vehicle  describing  function  for  a  high  bandwidth  single 
axis  pitch  forcing  function  with  instrument  display  and  rate  command  system. 
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Figure  5.31  Open  loop  pilot-vehicle  describing  function  for  a  low  bandwidth  single 
axis  yaw  forcing  function  with  instrument  display  and  rate  command  system. 
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Figure  5.32  Pilot  describing  function  for  a  high  bandwidth  single  axis  pitch  forcing 
function  with  instrument  display  and  rate  command  system. 
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Figure  5.33  Pilot  describing  function  for  a  low  bandwidth  single  axis  yaw  forcing 
function  with  instrument  'display  and  rate  command  system. 
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Figure  5.34  Open  loop  pilot-vehicle  describing  function  for  a  high  bandwidth  single 
axis  yaw  forcing  function  with  instrument  display  and  open  loop  command  system. 


5.3.3  Pilot  Compensation  and  Handling  Qualities 

The  slope  of  the  pilot  magnitude  curve  at  the  best  fit  open  loop  pilot-vehicle 
crossover  frequency  is  used  as  a  measure  of  pilot  compensation.  This  is  chosen  because 
it  unambiguously  indicates  the  type  of  compensation  generated  by  the  pilot  and  is 
not  affected  by  time  delay. 

Figures  5.36  and  5.37  show  how  HQRs  are  affected  by  the  amount  of  lead  com¬ 
pensation  for  the  single  and  multi-axis  visual  cases  respectively.  Second  and  first 
order  polynomials  are  fitted  to  the  experimental  data  and  are  given  in  Equations  5.40 
and  5.41  for  the  single  and  multi-axis  visual  cases  respectively. 

HQ  Rest  =  0.0017  +  0122  +  3.7925  [5.40] 

[dlog(u;)J  dlogiu) 

+  6.7776  [5.41] 

One  can  see  that  lead  compensation  degrades  HQRs  while  lag  compensation  improves 
them.  This  confirms  the  trends  observed  in  References  [11]  and  [9].  No  conclusions 
as  to  how  compensation  affects  HQRs  can  be  made  for  the  instrument  cases  because 
of  data  scatter. 


HQRest  =  0.1201 


dlog(u;) 
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Figure  5.35  Pilot  describing  function  for  a  high  bandwidth  single  axis  yaw  forcing 
function  with  instrument  display  and  open  loop  command  system. 
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Normalized  performance  {calcr,)  as  a  function  of  slope  of  the  pilot  magnitude 
curve  at  the  open  loop  pilot-vehicle  crossover  frequency  is  shown  for  the  single  axis 
instrument,  single  axis  visual  and  multi-axis  visual  tasks  in  Figures  5.38,  5.39  and  5.40 
respectively.  One  can  see  that  pilot  performance  improves  as  lag  generation  increases 
and  degrades  as  lead  compensation  increases.  Figures  5.38  and  5.39  indicate  that 
the  pilot  achieves  the  best  performance  when  generating  significant  amounts  of  lag 
{Slope  <  —IbdB/dec).  The  experimental  data  shown  in  Figure  5.40  suggest  that 
there  is  a  lower  limit  to  how  much  lag  the  pilot  may  generate  and  still  improve 
performance.  The  pilot  achieves  the  best  performance  when  generating  about  -10 
dB/dec  of  lag  compensation  and  performance  degrades  as  lag  compensation  increases. 
Figures  5.38,  5.39  and  5.40  all  indicate  that  pilot  lead  compensation  degrades  overall 
performance. 

From  the  above  analysis  one  may  conclude  that  the  preferred  pilot  model  is  of  the 
form 

"  "■  TugS  -h  1 

This  means  that  pilots  prefer  to  generate  lag  compensation  which  implies  that  attitude 
command  systems  are  desirable.  More  generally  a  desirable  controlled  element  Yc  in 
a  man-machine  system  should  have  the  characteristics  of  a  simple  tracker  (unity 
gain  element).  This  modifies  the  assumption  of  Reference  [l]  that  “human  control 
behavior  reaches  optimum  when  the  man  becomes  the  analogue  of  a  simple  amplifier” 
and  tasks  of  integration  (lag)  and  differentiation  (lead)  degrade  overall  performance. 
The  results  of  Experiments  ^2  and  ^3  validate  the  modification  to  the  £issumptions 
of  Reference  [1]  made  by  Reference  [11]  which  states  controlled  elements  which  require 
the  pilot  to  generate  low  frequency  lag  compensation  to  achieve  VpYc  ^  Ke~'^* js  are 
desirable. 
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Figure  5.37  HQR  as  a  function  of  slope  of  pilot  Bode  magnitude  curve  at  the  open 
loop  pilot-vehicle  crossover  frequency  for  the  single  axis  visual  cases. 
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Single  Axis  Instrument  Tasks 


Figure  5.38  Normalized  performance  as  a  function  of  slope  of  pilot  Bode  magnitude 
curve  at  the  open  loop  pilot-vehicle  crossover  frequency  for  the  single  axis 

instrument  cases. 
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Single  Axis  Visual  Tasks 
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Figure  5.39  Normalized  performance  as  a  function  of  slope  of  pilot  Bode  magnitude 
curve  at  the  open  loop  pilot-vehicle  crossover  frequency  for  the  single  axis  visual 

Ccises. 
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Muld  Axis  Visual  Tasks 


Figure  5.40  Normalized  performance  as  a  function  of  slope  of  pilot  Bode  magnitude 
curve  at  the  open  loop  pilot-vehicle  crossover  frequency  for  the  multi-axis  visual 

cases. 
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6.  CONCLUSIONS 


A  study  of  the  effects  of  command  systems,  displays  and  tasks  on  human  pilot 
behavior  is  conducted.  A  general  aviation  pilot  uses  the  Cooper-Harper  rating  scale 
to  evaluate  the  flying  qualities  of  the  TROFM  using  unaugmented,  rate,  and  attitude 
command  systems  for  two  types  of  tasks  and  two  different  display  types.  Pilot  and 
pilot- vehicle  describing  functions  are  obtained  from  time  histories  of  dynamic  tracking 
tasks  where  the  pilot  attempts  to  track  a  random  appearing  input.  Compensatory 
instrument  displays  and  a  visual  pursuit  scene  are  used  for  the  dynamic  tracking  tasks. 
Correlations  between  the  Cooper-Harper  ratings  and  quantitative  descriptions  of  pilot 
behavior  are  made.  The  experimental  results  are  compared  to  the  classical  theory  of 
man-machine  system  dynamics  which  is  based  on  '  IcRuer’s  Crossover  Model. 

The  pilot  prefers  attitude  command  systems  when  maneuvering  to  stationary  ob¬ 
jects  and  when  tracking  moving  targets.  This  preference  is  independent  of  the  display 
type  and  in  most  cases  the  forcing  function  bandwidth.  It  is  found  that  lag  compen¬ 
sation  is  ecisier  for  a  pilot  to  generate  than  either  lead  or  pure  gain  compensation. 
Handling  qualities  ratings  and  over-all  performance  degrade  as  the  pilot  is  required  to 
produce  more  lead  compensation.  This  result  implies  that  attitude  command  systems 
have  desirable  dynamic  characteristics  for  controlled  elements  in  pilot-vehicle  systems 
where  the  pilot  is  attempting  to  track  an  attitude  forcing  function.  More  generally, 
this  implies  that  desirable  controlled  elements  in  man-machine  systems  have  unity 
gain  (flat  0  dB  frequency  responses)  over  the  frequency  range  of  interest. 

The  pilot  achieves  the  best  performance  in  dynamic  tracking  tasks  when  using  a 
visual  pursuit  scene  as  opposed  to  a  compensatory  instrument  display.  It  is  found 
that  a  product  rule  describes  the  handling  qualities  degradation  associated  with  a 
multi-axis  visual  task  within  ±  1  HQR  for  78  %  of  the  cases. 
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The  extended  crossover  model  describes  most  of  the  experimentally  determined 
open  loop  describing  functions  well  in  cases  where  the  HQR  <  7.  It  is  also  found 
that  the  type  of  pilot  compensation  can  be  predicted  from  knowledge  of  the  machine 
dynamics,  forcing  function  bandwidth  and  application  of  the  theoretical  crossover 
model.  Pilot  ratings  are  found  to  be  repeatable  since  the  HQRs  for  a  given  configu¬ 
ration  are  consistent  (±  0.79  HQR). 
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7.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


Future  researchers  can  use  this  equipment  to  conduct  a  study  of  how  different  types 
of  manipulators  (joysticks,  etc.)  effect  performance  and  handling  qualities  ratings.  It 
is  recommended  that  a  spring  loaded  center  stick  be  one  of  the  candidates  because 
of  its  widespread  use  in  other  handling  qualities  studies.  Improvements  in  the  multi¬ 
axis  compensatory  instrument  display  are  also  in  order.  Studies  similar  to  the  one 
performed  here  can  make  use  of  command  systems  designed  by  modern  multi- variable 
methods.  The  response  and  adaptation  of  the  human  operator  to  step  changes  in  the 
TROFM  dynamics  can  also  be  studied.  This  can  simulate  the  in-flight  failure  of  a 
stability  augmentation  system.  Comparison  of  the  pilot’s  response  when  maneuvering 
from  one  point  to  another  to  nonlinear  time  optimal  “bang-bang”  models  can  also  be 
made. 

In  practice  tilt  rotor  aircraft  use  constant  speed  rotors  with  variable  pitch  blades 
to  govern  the  amount  of  lift  produced.  This  configuration  allows  variations  in  lift 
to  occur  much  faster  since  these  variations  are  not  limited  by  the  b^Lndwidth  of  the 
motor.  A  variable  pitch  rotor  could  enhance  the  TROFM  discussed  in  this  report. 
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Appendix  A:  Schematics  of  Electronic  Circuits 
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Figure  A. 2 


Schematic  of  connections  to  Copley  Model  201  Amplifier. 
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Figure  A. 3  PWM  circuit  for  tilt  servo. 


See  Reference  [19]  for  details. 
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Figure  A. 4  Stepper  motor  logic  circuit. 


See  Reference  [19)  for  details. 
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Appendix  B:  Runlog  for  Piloted  Experiments 


This  appendix  contains  the  run  log  for  the  man  in  the  loop  experiments.  The  lists 
and  tables  use  keys  to  indicate  the  configuration  of  the  TROFM.  The  general  form 
of  the  configuration  key  for  the  pilot  comment  list  is  “A/B:  C/D  C/D;  /E".  In  single 
axis  cases  and  individual  table  listings  the  form  is  “A/B/C/D/E”.  The  following  is  a 
break  down  of  the  configuration  key: 

.\.  Display  configuration 

I.  Instrument  Display 
V'isual  Scene 

B.  Number  of  axes  to  control 

S.  Single  Axis 
M.  .Multi  Axis  (2) 

C.  Axis  to  control 

P.  Pitch 
Y.  Yaw 

D.  Bandwidth  of  task 

L.  Low  (Slow) 

M.  Medium 
H.  High  (Fast) 

E.  Command  System 


A.  Attitude  (position) 

R.  Rate 

0.  Open  loop  (un-augmented) 
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B.l  Pilot  Comments 


The  following  are  pilot  comments  transcribed  from  audio  tape.  The  comments 
w'ere  recorded  immediately  after  each  run.  Because  of  a  recorder  malfunction  during 
the  first  twelve  runs,  no  audio  record  was  available  for  transcription.  In  those  cases, 
key  words  or  phrases  recorded  on  paper  by  the  flight  test  director  were  used. 

Run-1  I/S/P/L/A:  No  comments  recorded. 

Run-2  I/S/P/M/A:  Controllability  in  question.  Jerking  the  stick  makes  it  uncontrol¬ 
lable. 

Run-3  I/S/P/H/A;  If  I  jerk  the  stick  at  all  it  is  uncontrollable.  Smooth  stick. 

Run-4  I/S/V/L/.A:  Jerking  stick  still  a  factor.  Bad  display  (up-down  yaw  error  dis¬ 
play). 

Run-5  I/S/Y/M/A:  Low  workload  but  same  problems  as  above  (4) 

Run-6  I/S/Y/H/A:  No  comments  recorded 

Run-7  I/M;  Y/M  P/L;  /A;  Alot  of  concentration.  More  pilot  compensation. 

Run-8  I/M;  Y/H  P/L;  /A:  Alot  of  concentration.  Learning  curve. 

Run-9  I/M;  Y/L  P/M;  /A:  Poor  in  pitch. 

Run-10  I/M;  Y/H  P/M;  /A:  Uncontrollable,  both  (axes)  bad.  Let  error  slide. 

Run-11  I/M;  Y/L  P/H;  /A;  Adequate. 

Run-12  1/M;  Y/M  P/H;  /A:  Intense  concentration.  Very  high  workload. 

Run-13  I/S/P/L/R:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  I  don’t  like  the  control  characteristics  of  this.  It  is  slow  to  react  at  first 
and  too  quick  later.  .Aircraft  characteristics  were  objectionable,  but  you  could 
compensate  for  them.  I’ll  give  it  a  6. 
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Run- 14  I/S/P/M/R;  Controllable?  Yes.  Performance  not  adequate  but  felt  like  I  had 
more  instantaneous  control,  very  jerky  and  hard  to  keep  between  the  bounds. 
Alot  of  compensation.  I’ll  give  it  a  7  1/2. 

Run-15  I/S/P/H/R:  I’m  going  out  of  bounds  alot  but  not  because  it  is  uncontrollable, 
but  because  it  reacts  so  quickly  to  an  input.  You  can’t  keep  very  narrow  margins 
of  error  on  this  task.  It  was  definitely  controllable.  I  like  this  better  than  the 
other  ones  (13,14)  but  didn’t  do  as  well.  I’d  give  it  a  6  1/2. 

Run-16  I/S/Y/L/R:  Controllable?  Yes.  Performance  adequate?  Yes.  Satisfactory 
without  improvement?  I  don’t  like  this  type  of  system  that  well,  but  it  was 
an  easy  test.  Required  a  fair  amount  of  compensation.  The  deficiencies  were 
mildly  unpleasant,  but  moderate  compensation,  so  3  1/2. 

Run-17  I/S/Y/.M/R;  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  No.  Workload  was  a  little  bit 
higher  on  this  one.  If  the  rate  of  change  of  an  error  was  relatively  fast,  if  I 
gave  a  large  input,  1  almost  got  into  a  situation  similar  to  a  PIO  (pilot  induced 
oscillation).  That  is  objectionable.  Considerable  compensation.  I’d  give  it  a  5 
1/2. 

Run-18  I/S/Y/H/R:  Controllable?  Yes.  Performance  adequate?  marginal.  Big  prob¬ 
lem  with  overshoot.  I  liked  the  fact  that  the  responsiveness  wcis  better  on  this, 
but  I  kept  going  out  of  bounds  on  my  control  inputs.  I  would  consider  the 
overshoot  very  objectionable.  Considerable  pilot  workload  involved.  I’d  give  it 
a  6  1/2. 

Run-19  I/M/;Y/M;P/L; /R;  Controllable?  Almost  no.  Performance  adequate?  Defi¬ 
nitely  not.  At  best  it’s  a  9.  Two  tasks  at  once,  the  pilot  concentration  is  very 
intense.  It’s  hard  in  the  single  axis  C2ise  but  almost  impossible  to  do  here.  I’ll 
give  it  a  9  but  that's  being  generous. 
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Run-20  I/M/;  Y/H  P/L;  /R;  Controllable?  I  think  it  is  controllable,  actually  I  like  it 
better  than  the  last  test  we  just  did.  The  yaw  is  more  responsive,  almost  too 
responsive.  The  same  problems  we  had  before,  the  overshoot.  Pilot  concentra¬ 
tion  is  intense.  Deficiencies:  big  problems  with  overshoot,  problems  with  two 
things  at  once.  It’s  a  difficult  task.  I’ll  go  down  to  a  9  1/2. 

Run-21  I/M/;  Y/L  P/M;  /R:  Controllable?  No.  It’s  a  10.  I  had  to  let  one  axis  go  to 
correct  for  the  other.  Pitch  was  responsive  when  I  needed  it.  Trying  to  control 
for  yaw  was  tough  to  do. 

Run-22  I/M;  Y/H  P/M;  /R:  It  only  gets  worse,  another  10.  I  can  control  one  axis,  but 
keeping  the  other  one  within  bounds  is  next  to  impossible.  1  think  the  display 
has  alot  to  do  with  this.  If  the  yaw  was  left-right  and  the  pitch  was  up-down 
the  way  I  think  it  should  be,  this  would  most  likely  be  controllable. 

Run-23  I/M:  Y/L  P/H;  /R:  I  w’ouldn’t  say  uncontrollable,  but  pretty  close.  The  yaw 
error  was  kept  within  bounds  pretty  well,  but  the  pitch  error  kept  going  out. 
I’d  give  it  a  9  1/2. 

Run-24  I/M;  Y/M  P/F;  /R:  No  doubt  about  this  rating,  it’s  a  10. 

Run  25  I/S/P/L/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Not  much  required  in  the  way 
of  pilot  compensation,  alot  of  control  input  needed  in  order  to  control  this.  I 
would  say  the  deficiencies  are  minor  but  annoying.  I’ll  give  it  a  3  1/2  to  4. 

Run-26  I/S/P/M/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Alot  of  control  input  is  required 
for  the  task  given.  Have  to  have  an  input  for  every  error  shown  on  the  display. 
It  was  a  little  bit  harder  to  keep  the  medium  bandwidth  task  in  bounds  than 
it  was  the  slow  bandwidth.  Pilot  compensations  was  a  little  bit  more,  but  not 
overly  burdensome.  I'd  give  it  a  4. 
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Run-27  I/S/P/H/0:  Controllable?  Yes.  I  originally  said  no  because  of  the  excessive 
overshoots  involved  with  this.  The  major  deficiency  is  the  overshoot.  It’s  ex¬ 
tremely  hard  to  control.  The  pilot  compensation  is  not  intense  but  there  is  alot 
of  pilot  compensation  involved.  The  biggest  problem  is  the  overshoot  which 
would  cause  me  to  give  it  an  8. 

Run-28  I/S/Y/L/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  If  you  give  it  a  large  input, 
you  have  to  counteract  that  input  for  a  short  period  of  time  and  then  center 
the  joystick  in  order  to  keep  from  oscillation  back  and  forth.  Moderate  pilot 
compensation,  not  alot  going  on.  Minor  to  moderately  objectionable.  4  1/2. 

Run-29  I/S/Y/M/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  No.  Overshoot  a  very  big  problem  and  again  the  problem  of  overcoming 
your  previous  inputs.  Not  alot  of  pilot  compensation  involved.  You  have  alot  of 
inputs  going  in  but  it  doesn't  take  alot  of  thinking,  but  the  deficiency  is  major. 
I’d  give  this  one  an  8. 

Run-30  I/S/Y/H/0:  Same  problems  as  last  time,  workload  is  alot  higher.  I’ll  give  it 
a  9. 

Run-31  I/M;  Y/M  P/L;  /O:  Definitely  a  10,  it’s  uncontrollable. 

Run-32  I/M;  Y/H  P/L;  /O:  Controllable?  Barely.  Performance  adequate?  No.  I’ll 
give  it  a  9. 

Run-33  I/M;  Y/L  P/M;  /O:  Controllable?  No.  Response  was  too  sluggish  in  the  yaw 
axis.  I  was  concentrating  so  much  on  yaw,  I  forgot  about  pitch.  Definitely  a  10. 

Run-34  I/M;  Y/H  P/M;  /O;  Another  10.  I  had  trouble  with  yaw  again.  Trouble  with 
yaw  seems  to  translate  into  pitch. 

Run-35  I/M:  Y/L  P/H:  /O;  .Another  10.  Yaw  control  is  very  sluggish.  Pitch  seems 
OK  but  I  have  to  concentrate  on  yaw  too  much. 
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Run-36  I/M;  Y/M  P/H;  /O:  The  worst  of  them  all,  a  10.  1  was  pegged  at  the  limits 
for  most  of  the  test. 

Run-37  V/S/Y/L/R:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Yes.  Deficiencies:  a  little  hard 
to  keep  light  centered  on  line  during  initialization  phase.  I  think  this  is  because 
it’s  a  rate  controller,  at  least  I  think  it  is.  Deficiencies  are  not  a  big  deal,  e2ay 
to  control,  minimal  compensation.  I’ll  give  it  a  3. 

Run-38  V/S/Y/M/R:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Yes.  Little  bit  more  overshoot, 
harder  to  control,  I  don’t  like  it  as  much  as  the  last  one.  Borderline  satisfactory 
with  out  improvement  so  I'll  give  it  a  3  1/2. 

Run-39  V/S/Y/H/R;  Pretty  much  the  same  comments  as  the  last  one.  I  liked  the 
slower  bandwidth  tasks  better  though.  I’ll  go  with  a  3  1/2. 

Run-40  V'/S/P/L/.A:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Yes.  Pilot  compensation  is 
minimal,  very  easy  to  keep  the  light  centered  on  the  target  at  all  times.  I’ll  give 
it  a  1. 

Run-41  V/S/P/M/A:  Same  results  «is  last  time  but  quick  movements  cause  overshoot. 
Deficiencies  are  negligible  though.  I’ll  give  it  a  2. 

Run-42  V/S/P/H/A:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  No.  Overshoot  is  much  more 
pronounced.  Large  movements  cause  a  large  overshoot  which  is  hard  to  correct 
for.  Minor  deficiencies,  minimal  compensation.  I’ll  give  it  a  4. 

Run-43  V/S/P/L/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Yes.  Seems  slow  to  react  to 
commands  but  the  tracking  is  relatively  good.  Not  a  whole  lot  of  compensation 
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involved  provided  you  don’t  get  a  large  acceleration  going  where  you  have  to 
recover  from  it.  I'd  give  it  a  3. 

Run-44  V/S/P/M/0:  Pretty  much  the  same  comments,  overshoot  becoming  more 
pronounced,  but  the  reaction  time  is  better.  I  almost  prefer  this  task  more.  I 
give  it  a  3. 

Run-15  V/S/P/H/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Borderline,  overshoot  is  more 
pronounced.  I’m  catching  myself  trying  to  predict  where  the  target  will  stop 
moving  in  order  to  prevent  overshoot.  Compensation  is  a  little  bit  higher  but 
not  demanding  at  all,  but  the  overshoot  is  an  annoying  deficiency.  I’d  give  it  a 
3  1/2  to  a  4. 

Run-46  Y/S/Y/L/O;  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  ’ies.  Satisfactory  without  improvement?  There’s  a  trade  off  between  stick 
sensitivity  and  how  quick  you  react.  For  a  smaller  stick  sensitivity,  the  aircraft 
doesn’t  react  as  quickly  to  your  inputs  but  you  have  less  overshoot.  For  a  large 
sensitivity,  it  reacts  quickly  but  there’s  alot  more  overshoot.  It’s  not  in  the 
aircraft  itself  but  it  causes  a  problem.  There  is  a  little  bit  more  compensation 
involved  with  this  and  this  deficiency  is  a  little  annoying.  I’d  give  it  a  4. 

Run-47  V/S/Y/M/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Definitely  not.  It’s  extremely 
hard  to  keep  the  light  centered  on  the  target.  Alot  of  overshoot  involved  the 
whole  time.  There’s  more  pilot  compensation  involved  but  it’s  not  extensive 
but  would  say  it’s  considerable.  The  deficiencies  are  going  from  moderate  to 
very  objectionable.  I’d  give  it  a  5  1/2. 

Run-48  V/S/Y/H/0:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Debatable.  I  had  to  drop  the  sensitivity  on  the  stick  to  even  attempt 
to  stay  on  the  card.  Any  input  was  either  a  full  deflection  to  the  left  or  right 
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in  order  to  get  the  aircraft  moving.  Alot  more  compensation  involved.  Major 
deficiencies.  I’d  give  it  an  8. 

Rup-49  V/S/Y/L/A:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Yes.  In  order  to  get  maximum 
deflection  to  the  left  or  right,  the  joystick  sensitivity  has  to  be  set  a  little  bit 
higher  than  I  preferred.  It  seems  a  little  bit  jerky.  Other  than  that  it’s  a  good 
controller.  It  works  really  well.  I’d  give  it  a  3. 

Run-50  V/S/Y/M/A:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Yes.  As  the  bandwidth  in¬ 
creases,  overshoot  is  becoming  more  of  a  problem.  A  little  bit  harder  to  control 
,but  it’s  not  that  big  of  a  deal,  not  alot  of  compensation  involved.  I’d  give  it  a 
3  to  a  3  1/2. 

Run-51  V/S/Y/H/A:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Overshoot  again  is  a  problem. 
Not  alot  of  compensation  involved.  I’d  give  it  a  3  1/2  to  a  4. 

Run-52  V/S/P/L/R:  Controllable?  Yes.  Performance  adequate  with  tolerable  work¬ 
load?  Yes.  Satisfactory  without  improvement?  Tendency  to  overshoot  alot. 
It’s  a  little  bit  annoying.  I  came  off  the  card  a  little  bit  a  few  times.  It’s  not 
too  difficult  of  a  task  to  track  the  card  with  this  type  of  controller  though. 
Compensation  is  moderate.  But  satisfactory  without  improvement  no,  because 
of  minor  but  annoying  deficiencies  with  moderate  pilot  compensation.  I  give  it 
a  4. 

Run-53  V/S/P/M/R;  Pretty  much  the  same  comments  as  last  time.  Overshoot  is 
getting  to  be  more  of  a  problem,  cant  use  any  jerky  inputs.  Jerky  inputs  give 
you  a  large  amount  of  overshoot  and  you  spend  too  much  time  correcting  for 
them.  Again  the  deficiencies  are  annoying  to  moderately  objectionable.  Pilot 
compensation  is  moderate.  I’ll  give  it  a  4  1/2. 
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Run--54  V/S/P/H/R:  Same  problem  as  last  time  but  overshoot  is  getting  to  be  a  much 
bigger  problem.  I'm  having  trouble  keeping  it  on  the  card  at  some  points.  Pilot 
compensation  is  definitely  higher,  but  not  extensive.  The  overshoot  is  getting 
very  objectionable.  I’ll  give  it  a  5  1/2. 

Run'55  V/M;  Y/L  P/M;  A;  Controllable?  Yes.  Performance  adequate  with  tolerable 
workload?  Yes.  Satisfactory  without  improvement?  Yes,  I  am  having  a  little 
bit  more  problems  with  overshoot  when  we  get  into  this  multi-axis  case.  Not 
alot  of  compensation.  Overshoot  is  mildly  unpleasant  so  I’ll  give  it  a  3. 

»Run-56  V/M;  Y/M  P/H;  A:  Controllable?  Yes.  Performance  adequate  with  tolerable 
workload?  Yes,  I  only  went  off  the  card  a  couple  of  times.  Satisfactory  without 
improvement?  Little  bit  more  compensation  required  since  the  task  is  faster. 
Overshoot  is  a  problem  if  you  make  large  inputs.  I  don’t  have  trouble  doing  the 
task  but  there  is  something  about  it  I  don’t  like.  I’ll  give  it  a  4  1/2. 

Run-57  V/M;  Y/H  P/M;  A;  Controllable?  Yes.  Performance  adequate  with  tolerable 
workload?  Debatable.  Requires  extensive  pilot  compensation.  Deficiencies, 
overshoot  is  a  problem.  I’ll  give  it  a  5  1/2  to  a  6. 

Run-58  V/M;  Y/M  P/L;  R:  Controllable?  Yes.  Performance  adequate  with  tolerable 
workload?  Debatable.  Having  trouble  with  picking  joystick  gains.  Each  axis  is 
very  controllable  but  putting  them  together  makes  it  alot  more  difficult.  I’m 
having  problems  with  overshoot.  I’m  having  a  hard  time  keeping  centered  on 

^  the  card  but  I  can  keep  it  on  the  card  pretty  well,  but  my  tendency  is  not  to 
always  be  heading  toward  the  center  of  the  card,  I  may  be  heading  away.  I’ll 
give  it  a  7. 

• 

Run-59  V/M;  Y/M  P/L;  R;  Controllable?  Yes.  Performance  adequate  with  tolerable 
workload?  I  think  my  performance  w^ls  adequate.  It  was  definitely  more  diffi¬ 
cult,  alot  of  compensation  involved.  I  think  my  performance  was  better  on  this 
one  because  my  workload  was  higher,  I  went  with  what  was  natural  to  me  in 
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the  way  of  controlling  an  aircraft.  I  didn’t  think  as  much.  I  didn’t  get  confused 
as  much  with  the  (joystick)  sensitivities.  I’d  give  it  a  6  1/2. 

Run-60  V/M;  Y/H  P/M;  R:  Controllable?  Yes.  Performance  adequate  with  tolerable 
workload?  No.  Very  hard  to  keep  on  the  card  with  a  tendency  to  be  off  the 
card  for  a  very  long  time.  Pilot  compensation  was  maximum.  I  did  not  like  this 
at  all.  I’d  give  it  a  9. 

Run-61  V/M;  Y/M  P/L;  O:  Controllable?  Barely.  Performance  adequate?  Definitely 
not.  Almost  impossible  to  keep  on  card.  Very  slow  to  react  and  extremely  hard 
task  to  perform.  Being  generous,  I’ll  give  it  a  9. 

Run-62  V/M;  Y/L  P/M;  O:  Controllable?  Yes.  Performance  adequate?  No.  It  is  a 
difficult  controller  to  work  with,  alot  of  pilot  compensation,  slow  to  react  and 
once  it  finally  does  react,  I  need  to  counteract  that  input.  I  don’t  like  this  kind 
of  controller.  I  felt  more  comfortable  with  this  controller  than  the  previous  task 
(61).  I’ll  give  it  an  8. 

Run-63  V/M;  Y/M  P/H;  O:  Controllable?  Yes.  Performance  adequate?  I  don’t 
believe  it  was  adequate  but  definitely  better  than  the  past  (runs  61  and  62).  I 
have  no  idea  why  but  I  felt  more  comfortable  with  this  than  I  did  the  slower 
bandwidth  tasks.  I  could  keep  it  on  the  card  pretty  well  and  if  I  wasn’t  I  was 
always  heading  back  toward  it.  Deficiencies:  alot  of  overshoot  and  slow  to 
react  and  a  little  slow  to  recover.  I  don’t  thing  adequate  performance  could  be 
achieved  with  maximum  pilot  compensation.  I’d  give  it  a  7. 

Run-65  I/S/Y/L/0:  Controllable?  Yes.  Performance  adequate?  Yes.  Satisfactory 
without  improvement?  Yes.  Task  is  almost  hands  off,  a  “  no  brainer”,  no 
compensation.  Keeping  within  the  bounds  is  extremely  easy  to  do.  Again  there 
is  a  problem  with  the  overshoot  but  for  this  slow  bandwidth  it  doesn’t  matter. 
I’d  say  a  2  1/2. 
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Run-66  I/S/Y/M/0:  Controllable?  Yes.  Performance  adequate?  Yes.  Satisfactory 
without  improvement?  Very  little  compensation  involved  at  all.  Minimal  com¬ 
pensation.  Very  easy  task  to  perform.  Deficiencies  were:  slow  to  react  amd 
overshoot.  I’d  rate  this  around  a  4. 

Run-67  I/S/Y/H/0:  Controllable?  Yes.  Performance  adequate?  Yes.  Satisfactory 
without  improvement?  No.  Alot  of  overshoot.  It’s  real  tough  to  keep  it  on  zero 
so  I  just  aimed  to  keep  it  within  the  bounds.  For  every  input,  I  need  a  counter 
input.  Moderately  objectionable.  I’d  give  it  a  5. 

B.2  Describing  Functions  and  Performance  Measures. 

The  graphs  and  tables  of  this  section  are  the  results  of  the  man  in  the  loop  exper¬ 
iments.  Describing  function  data  and  the  corresponding  best  fit  describing  functions 
are  plotted  on  the  same  graph.  Key  performance  measures  and  best  fit  describing 
function  parameters  are  tabulated.  The  following  is  a  list  of  symbols  for  the  tables 
of  this  appendix.  (*  Indicates  definitions  are  taken  from  Reference  [13].) 

tg  One  standard  deviation  rms  value  for  tracking  error  during  the  run  (radians  of 
pitch  or  yaw  tracking  error  as  appropriate)  . 

e<,/(T,  Normalized  performance.  If  less  than  one,  the  pilot  is  effectively  reducing 
errors.  If  greater  than  one  the  pilot  is  making  matters  worse. 

Cg  One  standard  deviation  rms  value  for  joystick  deflection  during  run  (inches  of 
longitudinal  or  lateral  stick  as  appropriate). 

•jJc  Crossover  frequency.  Frequency  at  which  the  open-loop  pilot-vehicle  0  dB  line 
intersects  the  Bode  amplitude  asymptote  calculated  from  a  linear  interpolation 
between  the  two  describing  function  data  points  immediately  above  and  below 
crossover  (rad/sec).  See  Figure  B.2. 
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*V?m  Bode  open-loop  phase  margin  at  frequency  of  closed-loop  gain  crossover,  wd 
computed  from  a  straight  line  interpolation  between  the  two  describing  function 
data  points  immediately  above  and  below  Uc  (deg).  See  Figure  B.2. 

‘SlopCo/  Slope  of  Bode  open-loop  amplitude  asymptote  between  two  data  points  im¬ 
mediately  above  and  below  gain  crossover  frequency  (dB/dec). 

Tj  oi  Plant  open  loop  high  frequency  time  delay  parameter  from  the  exponential  r  ui 
(sec).  Computed  from  a  least  squares  fit  of  r  u;  -f  q/u?  to  the  residual  phase 
(experimental  phase  -  (-90”)). 

Oo/  Plant  open-loop  low  frequency  phase  droop  parameter  from  the  exponential  ajuj. 
Computed  from  a  least  squares  fit  of  ru,'  +  a/u:  to  the  residual  phase  (experi¬ 
mental  phase  -  (-90°)). 

Koi  —^ctat  Least  squares  estimate  of  open-loop  pilot-vehicle  gain  or  equivalently 
open  loop  crossover  frequency. 

HQR  Cooper-Harper  handling  qualities  rating 

Kp  Pilot  gain  computed  from  low  order  equivalent  system  fitting  procedure. 

Tiead  Low  order  equivalent  system  estimate  of  pilot  generated  lead  time  constant 
(sec). 

Tlag  Low  order  equivalent  system  estimate  of  pilot  generated  lag  time  constant  (sec). 

Tp  Low  order  equivalent  system  estimate  of  high  frequency  pilot  time  delay  parameter 
from  the  exponential  Tpuj  (sec). 

Qp  Low  order  equivalent  system  estimate  of  low  frequency  phcise  droop  parameter 
from  the  exponential  Qp/w. 

Slopep  Slope  of  low  order  equivalent  pilot  Bode  amplitude  curve  at  the  least  squares 
estimate  of  the  open-loop  pilot-vehicle  crossover  frequency. 


202 


Figure  B.l 


Definition  of  linearly  interpolated  crossover  frequency  and  phase  margin. 
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Pitch  Bandwidth  Yaw  Bandwidth  Command  System  Instrument  HQR  Visual  BQR 
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Table  B.3,  continued. 
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Figure  B.2  Open  loop  pilot-vehicle  and  pilot  describing  functions. 


217 


3''9(Yp)  (deg)  arg(YpYc)  (deg) 


Run  2  Open  Loop  Describing  Function 


o>  (rad/sec) 


Run  2  Pilot  Describing  Function 


0)  (rad/sec) 


Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 


266 


arg(Yp)  (deg)  arg(YpYc)  (deg) 


0.1  0.2  12  10 
(a  (rad/sec) 


Rim  59  Pilot  Describing  Function 


(0  (rad/sec) 

Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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Figure  B.2,  continued. 
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